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ABSTRACT 


Whistler mode waves that propagate through the magnetosphere 
exchange energy with energetic electrons by wave-particle interaction 
mechanisms. Using linear theory, a detailed investigation is presented 
of the resulting amplitude variations of the wave as it propagates. 
Arbitrary wave frequency and direction of propagation are considered . A 
general class of electron distributions that are nonseparable in particle 
energy and pitch-angle is proposed . Comparison with data is obtained 
by computing the total amplitude variation between two locations along 
the wave ray path. It is found that the proposed distribution model is 
consistent with available whistler and particle observations. In 
particular, this model yields insignificant amplitude variation over a 
large frequency band, a feature commonly observed in whistler data. 

This feature of the data implies a certain equilibrium between waves and 
particles in the magnetosphere over a wide spread of particle energy, 
at least during certain (magnetically quiet) times, and is relevant to 
plasma injection experiments. Application of our analysis for monitoring 
the distribution of energetic electrons in the magnetosphere is discussed. 
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DEFINITIONS OF IMPORTANT TERMS USED THROUGHOUT THE TEXT 
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I- INTRCDUCTIGN 


The magnetosphere extends many earth radii above the ionosphere. It 
is the region in space where the earth’s magnetic field has dominant con- 
trol over the motions of charged particles. 

The presence of these charged particles in turn determines the 
characteristics of the different electromagnetic wave modes which can 
propagate in the magnetosphere- Our interest will be focused on the 
"’whistler"' mode. The frequency range of the whistler mode extends approx- 
imately between the proton and electron gyrof requencies . Its sense of 
polarization corresponds to the sense of gyration of electrons compelled 
to gyrate and drift along the earth’s magnetic field lines. 

The charged particles can be divided into two classes according to 
their energy: the ’cold’ or ’thermal’ particles and the ’hot* or ’ener- 

getic’ particles. The cold particles have a distribution approximately 
Maxwellian with an average energy of the order of 10 ^ eV [Angerami, 1966] 
and determine the ray path of a wave between two locations in the magneto- 
sphere. The hot particles can generate waves [Kennel and Petscheck, 1966; 
Helliwell, 1967] or cause amplitude variations of a preexisting wave along 
its ray path. Such an interaction involves a possible exchange of energy 
between waves and particles. This transfer of energy into the wave can 
be either positive or negative, depending upon the precise values of wave, 
cold, and hot plasma parameters. 

Since Storey’s [1953] investigation, the study of whistler wave prop- 
agation characteristics has provided an invaluable diagnostic tool to 
determine the distribution of thermal particles in the magnetosphere. 

In contrast, the study of whistler-electron interaction characteristics 
has not led to the development of a successful diagnostic tool for the 
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determination of the energetic particle population- The importance of 
this interaction study is considerable though- Such a study may provide 
a better knowledge of the energetic particle population- It may provide 
also a better understanding of the magnetospheric medium. The interaction 
mechanism in itself may furnish an explanation of features of the magneto- 
sphere or the generation of electromagnetic waves. 

Knowledge and understanding of particle distribution characteristics 
are of fundamental importance in the whole field of plasma physics, from 
laboratoiy plasma experiments to astrophysics. In the laboratory, for 
example, knowledge of particle distributions is primordial in the field 
of plasma ion sources [Bernard, 1967; Benoit-Cattin and Bernard, 1968]. 
Understanding of wave-particle interaction mechanisms is most important 
in fusion research- These mechanisms certainly play a role in astro- 
physical phenomena, like pulsars. 

The purpose of our work is to compute the whistler amplitude growth/ 
damping rates due to whistler-electron interaction, to discuss its varia- 
tions with regard to all parameters of concern, and, using a number of 
whistler observations, to infer possible models of the energetic electrons 
of the magnetosphere. 

Wave amplitude variations are due both to "instability*^^ and "accessi- 

t»2 

bility and it is often difficult to separate the two effects. When there 
is a frequency cutoff of a broadband signal, the question arises whether 
the frequencies above (or below) the frequency cutoff were damped by 

^"instability, as defined here, means wave amplitude variations due to 
wave-particle interaction . 

2 

A location in space is "accessible** if, according to geometrical optics, 
it can be reached by a wave starting from another point. 
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wave-particle interaction, or whether these frequencies could not propa- 
gate from the location of generation to the location of reception. This 
frequency cutoff phenomenon is illustrated in Figure 1.1 which represents 
an example of a ’’nose whistler" spectrogram recorded at a ground station. 
Typically, such a signal presents an upper cutoff at some frequency f^. 

The spectral shape of the signal is explained schematically in Figure 1.2 . 

A broadband wave signal is generated by a lightning discharge close to 
the earth’s surface. The generated frequencies are traiismitted in part in 
the earth-ionosphere waveguide to the receiver in the other hemisphere 
where they appear on the spectrogram practically without dispersion as a 

S 

frequency impulse or ’’atmospheric.” The generated frequencies travel also 
in the magnetosphere and are guided or ”ducted" along a magnetic field 
line. The magnetosphere is highly dispersive, that is, different fre- 
quencies travel with different velocities. As a result, the whistler has 
the characteristic shape shown on Figure l.l^ Notice from the figure that 
there is a characteristic frequency, namely the frequency of minimum time 

delay t or "nose" frequency f . Measurements of t and f have 
n n n n 

provided an invaluable technique for magnetospheric diagnostics [Smith, 

1960; Carpenter, 1962, 1963, 1966, 1970; Carpenter and Smith, 1964; 

Carpenter et al., 1972; Helliwell, 1965; Angerami, 1966; Angerami and 
Carpenter, 1966; Park, 1970, 1973; Park and Carpenter, 1970; Bernard, 1973; 
Ho and Bernard, 1973]. 

Smith [1961], proposed an elegant explanation of the upper frequency 
cutoff f^* As the magnetosphere is practically a collision free medium, 
density irregularities can exist for a long time. He showed that small 
field -aligned enhancements of ionization or "ducts" can trap low frequencies 
below (path 2 on Figure 1.2), whereas high frequencies above f^ aie 
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FIGURE 1.1. TYPICAL NOSE WHISTLER SPECTROGRAM SHOWING AN UPPER FREQUENCY CUTOFF. 
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FIGURE 1.2. RAY PATHS INVOLVED IN PRODUCTICffl OF A TYPICAL NOSE WHISTLER. Frequencies generated 
at point G by a lightning discharge travel in part close to the earth along path 1 
and in part along a magnetic field line along path 2 up to an upper frequency f 
above which they propagate along path 3. Wave components of frequency below f can be recorded 
on the ground to point R. Frequencies above f can only be recorded on satellites. 


untrapped at some point A (path 3) and cannot reach the receiver location 
R. 

Scarf [1962] proposed an alternative explanation for the upper cut- 

off f based on hot plasma effects* Taking an isotropic Maxwellian 
u 

distribution for the hot electrons he found that whistlers could be 

attenuated due to wave-particle interaction, the attenuation increasing 

rapidly with frequency- Postulating that Hie high frequency cutoff was 

indeed caused by such a mechanism, he determined in turn what should be 

5 o 

the temperature. He found an order of 10 K. This analysis was subse- 
quently refined by Liemohn and Scarf [1964]- 

5 o 

Guthart [1964, 1965] then showed that a temperature 10 K would 
produce a slight, but a measurable change in the dispersion curve f(t) 
of Figure 1.1. Because he could not observe this change, he proposed 
another type of distribution which would produce the observed cutoff but 
no change in the curve f(t). Liemohn [1967] made a quantitative analysis 
where he integrated the variation of amplitude along the field line path 
(2). He chose a particle distribution function f(W,a) which was a product 
of a function of W alone and a function of a alone where W and a 
are respectively the particle energy and pitch-angle: 

f(w,a.) oe w""^ sin\ (1.1) 

and found that a value of the parameter q = 2 would produce the right 
value of frequency cutoff. He found also an increasing in amplitude just 
below f^, compatible with observation. The question of whether this 
cutoff is caused by accessibility or wave-particle instability was reopened 
when Carpenter [1968] made a statistical study of over 500 nose whistlers 
propagating along a wide variety of locations. He showed that the ratio 
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f u/f eg (where fe^ is the equatorial gyrof requency of the path) was 
statistically very constant: ~ - 0*03- This value is 

almost exactly the value predicted by Smithes [1961] theory of trapping 
(f^/feg = 0-50). These facts favored the accessibility explanation 
rather than the instability explanation because the predicted value of 
the accessibility ratio practically independent of any 

parameter whereas the instability ratio depends upon many parameters 
and should vary with different propagation paths* This is contrary to 
observation . 

Evidence of the accessibility explanation was put forth by Angerami 
[1970], who observed some unducted frequency components of whistlers on 
a satellite. Briefly, referring to Figure 1.2, high frequency components 
are trapped up to a certain point A on the field line and become untrapped 
after, that is, no longer follow a field line but follow a path such as 
(3). A satellite intersecting the ray path (3) can record these components 
(see Figure 4*1) . 

Since the validity of the accessibility cutoff appears to be clearly 
established, the distribution proposed by Liemohn [1967] seems somewhat 
doubtful. One of our aims has been to find more realistic distributions. 

Furthermore, Liemohn *s [1967] computations include only propagation 

parallel to the magnetic field. Now that a body of very interesting 

satellite data, such as Angerami [1970], has been obtained showing 

evidence of nonparallel propagation, these computations should include 

nonparallel propagation as well. A qualitative study of wave instability 

for nonparallel propagation limited to low frequency (f « f ), has been 

c 

done by Kennel [1966], and Kennel and Thorne [1967]. Brinca [1972] extended 
this study to waves of frequency ~ Thorne [1968] made a quantitative 
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instability study for *'mr whistlers*' (see Section 4E) which are again 
low frequency whistlers. The observations of whistler signals at high 
normalized frequencies [Angerami, 1970; Dunckel and Helliwell, 1972], 
show the necessity to make a quantitative study of the high frequency 
oblique whistler instability. This study has been our second goal* 

The contributions of the present work can be briefly stated as 
follows: 


1. We propose a new model for the energetic particle distribution 
function. This model has some important characteristics and 
among them, it yields very small amplitude variations over a 
large band of wave frequencies, and we believe it represents 
the energetic electron distribution more realistically than 
previous models. 

2. We derive a new general expression for the variation of wave 
amplitude caused by wave-particle interactions for an arbitrary 
angle of propagation. This expression is algebraically simpler 
than previously derived expressions and permits easier numerical 
computations. A Fortran program has been developed to compute 
gain rates of a whistler integrated along its ray path for a 
certain class of hot plasma distributions. 

3 . We have integrated gain ra tes along whistler ray paths deduced 
from some whistler observations. We have determined some 
bounds on the values of the parameters of our proposed distri- 
bution. We make some suggestions concerning the explanation of 
certain features of the observed particle spectrum. 

4. We discuss the application of our analysis to the development 
of a diagnostic tool for monitoring the distribution of hot 
electrons, and to plasma injection experiments. 

In Chapter 2, the theory of whistler-mode wave particle interactions 
is presented. First, we describe briefly the physics of the interaction, 
followed by the formulation of solving the coupled system of Maxwell 
equations and linearized Vlasov equation. Then we derive a general 
expression for the rate of variation of wave amplitude for an arbitrary 
angle between the wave normal and the earth’s magnetic field and for an 
arbitrary frequency below the gyrof requency . 
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In Chapter 3, we study in detail the variation of the wave gain for 
parallel propagation and for both a distribution separable and nonseparable 
in energy and pitch angle* We present a detailed study of the influence 
of the angle between wave normal and earth*s magnetic field* 

In Chapter 4, we compute gain rates along ray paths deduced from 
some whistler observations - 

Finally, in Chapter 5, conclusions are drawn and recommendations are 
made for future work. The details of numerical computations of gain 
rates, and various derivations and auxiliary material are presented in 
the appendices - 
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II- THEORY 


A. INTRODUCTION 

This work is concerned with the amplitude variations of an electro- 
magnetic wave caused by the presence of energetic electrons* 

The physics of such a mechanism can already be understood from the 
interaction of a single particle with an electromagnetic wave- We describe 
the interaction from this point of view in Section B- We ^'define there 
such notions as ’’resonance” and "trapping.” It is then possible to 
establish a qualitative relation between the interaction of one particle 
to the interaction of a distribution of particles. 

More rigorously, the interaction between wave and particle is a 
solution of a wave -plasma system. As such, it involves the solution of 
Maxwell equations coupled with a kinematics equation. Since the magneto- 
spheric medium can generally be treated as collision-free, the evolution 
in time of the particle distribution function can be described by the 
Vlasov equation. 

Section C describes the formulation for solving the coupled system 
of equations after linearization of Vlasov's equation. 

In Section D, a general expression is derived for wave amplitude 
variation, including relativistic effects, arbitrary wave frequency, and 
arbitrary direction of propagation. 

B. PHYSICS OF THE INTERACTION DESCRIBED FROM A TEST PARTICLE MOTION 

1 . Interaction of One Particle with an Electromagnetic Wave 

The motion of a particle interacting with an electromagnetic wave 
has been described in great detail by a number of authors (e.g., Roberts 
and Buchsbaum [1964]; Laird and Knox [1965]; Bell [1965]; Lutomirski and 
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Sudan [196G]; Laird [1968]; Dungey [I960]; Roux and Solomon [1970]; 
Palmadesso and Schmidt [1971]; Dysthe [1971]; Palmadesso [1972]). In 
this section we will give a brief description of the interaction. 

Because of its motion a particle can experience the electric field 
of an electromagnetic wave, which is time-varying in the laboratory frame, 
as a constant electric field- In general this condition is realized only 
when the momentum or a component of the momentum of the particle takes a 
particular set of values. This condition is referred to as the resonance 
condition. Because the particle sees a constant electric field at 
resonance, it can experience a strong acceleration (deceleration) which 
means there is an energy exchange between the particle and the wave as 
the particle gains (loses) energy at the expense (benefit) of the wave 
energy. This mechanism is responsible for damping (growth) of the wave. 

Off resonance, the particle sees a time-varying field and therefore 
it alternatively experiences both acceleration and deceleration; the 
exchange of energy between the wave and the particle averages to zero as 
time elapses. A particle initially at resonance will not stay in this 
state for a very long time, its momentum changing its value due to its 
acceleration (deceleration). It is shown that, provided the value of the 
particle momentum is close enough to one of the resonance values, a stable 
situation called trapping can develop: a particle initially in phase with 

the wave field will be accelerated, will eventually see a phase reversal 
of the field, be decelerated and will again be in phase with the wave field, 
the process repeating itself in a stable way. 

Not too close to resonance, the particle will not be trapped, that 
is it will drift along successive peaks and valleys of the wave. 
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In any case, the interaction is most important at resonance and its 
vicinity because the further from resonance the less the particle is 
perturbed by the wave . 

For this reason, we rederive the resonance condition [Bell, 1964] 
which is deduced in the linear treatment. We also rederive the expression 
for the minimum energy of resonance. 

Linear behavior : The electric field of the wave is taken of the 

form E = E Within the linear approximation, we evaluate 

the phase of the field at the unperturbed (E = 0) location of the particle 
at the time t [Quemada, 1968]. According to what we have said previously, 
the resonance condition will be obtained when the phase of the field becomes 
constant. The unperturbed motion of the particle is defined by: 

dp p 

5T - “e bx 2o - X E 


where the relativistic mass 


/ 


M-Mq «/■■•+ 22 

MqC 


( 2 . 2 ) 


and 


% M 


(2.3) 


The motion of the particle (neglecting radiation) is a helix whose 

axis is parallel to therefore: 

“0 


= PzO 




p^(t) = cosio^t , ^ x(t) = ^ 

c 


p (t) = p sinoi t 
y 1 c 


= - SET '“““o' ' ^0 

C 


(2.4) 
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M/ith k = (k ,0,k ), the wave phase becomes: 

— X z 


Vi Pqz 

kT - cot = k*r^ + sinm t + (k — — -cS)t 

OMo) c zM 


(2.5) 


where r^=rr(t=rO). 

Using the identity: 


ixsincp ^ T / \ 
e ^ = S J (x)e ^ 


m=-cn 


m 


(2.6) 


ik a 

l(kT'-03t) 0 

e = e S J 


m \ Mcp / 


^Oz 

exp(i(k — — + mm “•m)t) (2.7) 

z M c 


This last expression shows that the particle sees a superposition of 
waves propagating in the Oz direction and sees a constant field when 

Pr 


Qz 
^z M 


+ mm “ o) = 0 
c 


( 2 . 8 ) 


Now, m' = m - k 


Oz 


is the Doppler shifted frequency of the wave 


z M 

seen by the particle and the resonance condition is rewritten as: 


m' mm , m=:0,± 1,±2 , • . 


(2.9) 


The resonance m = 0 is excited by both the parallel and the per- 
pendicular components of the electric field (referenced to the static 
magnetic field) and includes the well known Landau resonance. (For con- 
venience we will refer to the m = 0 resonance as the Landau resonance.) 
All other resonances m ^ 0 are excited by the perpendicular component 
of the electric field and are called the cyclotron resonances. The 
infinity of harmonics is due to the spatial variation of the electric 
field in planes perpendicular to the static magnetic field. 

In the general case, there are two roots of Eq. (2.8) (m ^ 0; 
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m = 0 only one positive root): 


^Rml,2 



( 2 . 10 ) 


The minimuin energy of resonance W (p = 0) is given by: 

K X. 


Rni 0 


1 + -■ 1 


Wq = MqC 


( 2 . 11 ) 


Pq = MgC 


0=0 (k.^) 


A = cd/co where here co is the rest mass gyrof requency 


cO cO 

N = refractive index 
N = Ncos6 


In the non-relativistic limit the resonance condition is satisfied 
for one value of for each m: 


= AlE, ^ 

RmO A N,, 


( 2 . 12 ) 


In Figure 2.1 we have plotted versus N„ with A a 

parameter for the Landau (m = 0) and the fundamental cyclotron (m - 1) 

resonances. Note in Figure 2.1 that departure from a straight line 

represents departure from the non-relativistic limit and that in the 

non-relativistic range W _ is lower (higher) than W when A is 

Ro R1 

smaller (larger) than 0.5. 

Nonlinear behavior ; To see what happens after the initial (linear) 
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FIGURE 2.1, A PLOT OF THE MINIMUM ENERGY OF INTERACTION OF AN ELECTRON WITH A WHISTLER 
VERSUS THE COMPONENT OF THE WHISTLER REFRACTIVE INDEX ALONG THE STATIC 
MAGNETIC FIELD (Nj, ) * Landau and fundamental cyclotron resonances are 
shown for different values of A. 




development of the interaction, we have to integrate the equations of 
motion. Though the general case (fi / 0) is fairly involved and has only 
been tackled recently [Palmadesso, 1972], it is basically the super- 
position of Landau and cyclotron cases and it suffices to look at these 
two cases separately. 

In the Landau interaction (E field alone), the equation of motion 
is simply; 

,2 

— ^ ^ E sin(kx-cot) (2.13) 

dt^ “o 

which becomes by change of variable E = kx-ot>t; 


d E qEk . _ 

_| . sin5 

dt 0 


(2.14) 


This equation is analogous to a pendulum equation. Note that ^ 
is merely the deviation from the resonance position ^ = — • t . 

Now the solution of Eq. (2.14) can be of two types depending upon 
initial conditions. The particle can oscillate back and forth around 
the resonance position or can ^^rotate” completely. In the first case 
the particle is said to be trapped. Trapping is an important phenomenon 
because the motion of the particle is significantly altered by the 
presence of the field and energy exchanges with the wave can be quite 
important . 

In the case of trapping, the solution is periodic. 

The case of small oscillations gives an order of magnitude of the 
period (called trapping time T ): 

It 


T 


L 




(2.15) 
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Note that particles close to resonance are most likely to be 
trapped (small oscillations case) • Note also that Landau trapping means 
bunching of particles in space - 

The situation will be found to be similar for cyclotron interaction- 
Consider a particle in the field of a circularly “polarized electromagnetic 


wave propagating along 

Let B, = [B, ^cosr (kz-cot) , B, .sinr (kz-<Bt) , 0] be the magnetic field 
— 1 10 -LO S 

vector of the wave. (f. - ±1, which defines polarization for normal 

interaction with either an electron or an ion.) 

The equations of motion of the particle are [Dysthe, 1971]: 


V „ = 


r V oi ^sin\lf , 
± cl ^ 


cn - “ C a:) C ^ 


(Vn -V^) 


cl v 


cos\]; 


(2,16) 


where the particle velocity is v = [v cosco,v sinco>Vj| ] and where 

eB eB 

♦ - <) J ■ “cO = IT ' “cl = — • “S' 

the angle ijf can be written as ^ = - + ^^(kz-o>t) and two differenti- 

ations with respect to time produces the expression 


= “ CD 


+ r kv 

^s 


(2.17) 


Since O) « O) ^ in the magnetosphere, the last equation in (2.16) 
cl cO 

shows that (provided v^ is not too small) co constant and therefore 
Eq. (2.17) becomes \(f ^ Thus the first equation in (2.16) can be 

rewritten 

^ sinijf (2.18) 


where 


2 m , , .1/2 

^ ^ -.kv ) ' 

TC T Cl L 

C 
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We find again a pendulum type equation, this time with respect to 
the variable that is, in velocity space. The same conclusions apply 

as in the Landau case. There will be trapped particles and the period 
of trapping will be of the order of T^ ^ . 

From the first two equations of (2.16) a constant of motion is 
immediately found [Dysthe, 1971]: 


2 

V 


“ 2v,rV 


P 


:= constant ♦ 


(2.19a) 


2 . Relation of Particle Motion to Wave Growth 

Now we would like to know what happens when we include an 
ensemble of particles. It is possible to get a qualitative answer with- 
out going to the full Maxwell Vlasov treatment described in some detail 
in the next section. 

In the Landau interaction case, particles initially having energy 
slightly above the energy of resonance will tend to be trapped and on the 
average oscillate at the resonance energy while it is the opposite for 
particles initially having energy slightly below resonance. Therefore, 
the former particles will tend to lose energy while the latter will tend 
to gain energy. If more particles tend to lose energy than to gain 
energy, there will be wave growth and vice versa. 

In terms of the distribution function: 


> 0 GROWTH 
< 0 ^ DAMPING 

For the cyclotron interaction case, it is useful to recall 
Eq. (2.19a). It is written also as: 
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(2.19b) 


1 Mv^ + i M(v, -V = constant = W' 

2 i 2 " p 


Equation (2.19b) is very easily interpreted. In the wave frame (where 
we shall denote quantitites by prime superscript), by Lorentz transforma- 
tions [Feynman, 1964]: 


E' 

—II 


Si 


E' 
— i 


(E + V XB) 
-P - X 



Now from Maxwell's equations: 


E =5 - V X B, or, e; = 0 . 

-X -P -X -X 


Also since Ej| = 0, 0 . 

Therefore Eq. (2.19b) simply means that the particle energy W' 

in the wave frame stays constant. Following Gendrin [1968], Eq. (2.19b) 

is graphically interpreted in Figure 2.2 . 

In the (V|| ,v^) plane, particles of constant energy follow a circle 

centered at the origin. From Eq . (2.19b) particles interacting with the 

wave follow a circle centered on v (taking v as > 0; remember that 

P P 

V < 0). From Figure 2.2, it is obvious that particles with parallel 
R 

velocity initially slightly above |v^| and resonating with the wave 
will gain energy. At the same time their pitch angle increases. The situ- 
ation is reversed for particles with slightly below 

cyclotron interaction, the ^slow* particles cause wave growth in contrast 
to the Landau interaction, a condition already drawn by Bell [1964]. This 
reflects merely the fact that for (normal) cyclotron interaction, wave and 
resonating particles are moving in opposite directions. For anomalous 
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FIGURE 2.2. GRAPHICAL INTERPRETATION OF THE FIRST CONSTANT OF MOTION OF AN 
ELECTRON INTERACTING WITH A V-HISTLER WAVE PROPAGATING ALONG THE 
STATIC MAGNETIC FIELD. 


are of 


cyclotron interaction (e-L or p-R interaction) v and v 

the same sign and ’fast* particles cause wave growth. In this case 

particle energy increase is associated with pitch angle decrease, which 

was algebraically demonstrated by Brice [1964] . 

From Figure 2.2, it is also seen that resonant particles with high 

energy (|v I » v ) interact with the wave with practically no change of 

energy because the circles of W' = constant are practically coincident 

with the circles W = constant. There is ’pure pitch angle* diffusion 

[Kennel and Petschek, 1966]. For particles with |v^j comparable or 

less than v , the curves of constant W and differ significantly. 

P 

There can be energy diffusion at a rate comparable to pitch angle 
diffusion. This is confirmed quantitatively by the quasi-linear theory 
[Kennel and Engelmann, 1966]. 


Now 



( 2 . 20 ) 


Therefore high-energy resonant particles interact with low frequency 

(cu « CD ) ivaves and vice versa, 
c 

We have seen that in the Landau case, gain rates are related to the 
derivative of the distribution function measured at the resonant velocity. 

We expect a similar situation in the cyclotron interaction but in con- 
trast to the Landau case which was clearly a one -dimensional problem, 
several derivatives should be involved. The azimuthal variable plays a 
fundamental role in trapping but in the linear case this variable is not 
of primary importance (we shall find more rigorously in the next section 
that a necessary condition for the time invariance of the unperturbed system 
is that f must be cp^independent) . So two variables remain to be considered, 
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V 


and for example. However, since in experimental measurements of 

particle fluxes, the energy W and the pitch angle a of the particle 
are most readily measurable, we shall use rather those two variables. 

From Figure 2.2 we can qualitatively see how the signs of the deriva 
tives and ^ influence the wave gain. On Figure 2.2, particles 

(1) and (2) refer respectively to particles with | V|| j slightly faster 

and slower than |v I. Either particle resonating with the wave follows 
' R ' 

a curve of constant W'. The trajectories can be decomposed as paths at 
constant pitch-angle, A B- and A B , and paths at constant energy 

JL JL ^ ^ 

B C and B C . 

1 2 

At constant pitch angle : particle (1) causes damping while particle 

(2) tends to cause growth. Now particle (1) has a lower level of energy 
than particle (2). Therefore, at constant pitch angle, if f is such 
that there are more particles of type (1) than of type (2), there will 
be damping, in other words: 


Sf 


< 0 DAMPING 


At constant energy ; particle (1) causes damping while particle (2) 
causes growth. More particles of type (2) than of type (1) will cause 
growth; 

df I 

^ j > 0 -* GROWTH 

This gives a qualitative picture of wave gain by wave^particle inter 
action, a result found more rigorously in the next section. The correct 
expressions for gain rates involve such derivatives of the particle dis- 
tribution function evaluated at the resonant velocities. The gain 
expression depends upon the number of resonant particles as anticipated, 


SEL 73-043 


22 



but because it depends also upon derivatives of f not only the ampli- 
tude of the distribution function is important but also the fine details 
of the distribution function- For this reason there is possibly a great 
range of wave gain rates while still considering realistic models of 


particle distributions in the magnetosphere. 

3 , Whistler Mode Refractive Index Characteristics and the Standard 
Approx ima t ion 

To know the energies of resonance of particles, we need the 
refractive index values. In a cold collisionless plasma immersed in a 
static magnetic field electromagnetic waves of frequency below the 

electron gyrof requency can propagate in either of two modes. These modes 
are elliptically polarized and the mode whose sense of polarization 
corresponds to the sense of gyration of electrons is defined as the 
'whistler mode." The other mode corresponds to a polarization in the 
opposite sense to the sense of gyration of the electrons. As we will be 
concerned essentially with electron interaction and wave of frequency 


^ we will be concerned uniquely with the whistler mode. It can 
ci 

be shown that in fact throughout most of the inner magnetosphere the 
other mode does not propagate. 

The full expression of the refractive index is given in Stix 
[1962]. A first approximation is to neglect the ions. It is valid when 
where is the lower hybrid resonance frequency. For a 


LHR 


high density, two species plasma consisting of electrons of mass 


and ions of mass 


“lhr % 



( 2 . 21 ) 


Neglecting ions, 


can then be written in the form 



N = 1 + 


/~2 2 2 2 2 ’ 
'Tl cos 6 + C sin 0 + Csin 6-A 


( 2 . 22 ) 


where 




Equation (2*22) shows that the whistler mode is a propagating mode 
over the angular range from 0=0 to the resonance cone angle 


. -1 A A2 , 2 

= cos -/P + 1 - y\ 


(2.23) 


High-plasma frequency approximation ; Equation (2.22) can be simpli- 
2 2 2 

fied when g » 1 or 03 » oi • With this approximation cos0 ^ A 

P R 

and C can be neglected in Eq . (2.22) as well as the factor 1 in front 
of the fraction terra: 


A(cos0-A) 


(2.24) 


Quasi-parallel approximation : This approximation holds when Csin 0 

2 2 

can be neglected in Eq . (2.22), i.e., |c|sin 0 « A • It is related to 

the "quasi-longitudinal" approximation of Stix [1962] and Helliwell 

[1965] or "quasi-circular" approximation of Allis et al [1963]. It 

holds for small 0 or for |c|. « A^» that is, 2(p^-A^) » 1 or 
2 

p » 1* This shows that the quasi-parallel approximation can be con- 
sidered as a particular case of the high plasma-frequency approximation. 
Non-relativistic approximation and standard approximation : When 

2 

6 » 1 ) relativistic corrections can be neglected for Landau resonance 

energy and, except in the case of extremely low frequencies, for funda- 
mental cyclotron resonance energy. For that reason, the non-relat ivis tic 
and the high plasma frequency (and a fortiori the quasi-parallel) 
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approximations overlap almost completely • From now on, we will use the 
term '’standard approximation” to mean explicitly that the three approxi^ 
mations are satisfied. 

Ion correction for very low frequencies: For u) comparable to 

(n , the ions cannot be neglected but a simple approximation can be used 
LHR 

[Edgar, 1972] within the standard approximation: 


N 


2 


2 


B 

A(cos0-gA) 


(2-25) 


where : 



For CJD » CO , we find expression (2-24). Equation (2-25) shows that 
LHR 

propagation is possible for all angles when oi < to „ . N(6) is sketched 

LiilK 

in polar coordinates on Figure 2.3 to point out the different topologies 
of the refractive index surface according to frequency. Precise values 
of N = NCOS0 are plotted versus 0 for different values of A> 

Figure 2.4 . 

Knowledge of the polarization of the whistler mode will be important 


also . 


E . ft a E Do 

z N sin0cos0 y _ . S 

N^sin^e-P N^-S 


(2.26) 


where P, D , and S are defined later in Eq. (2.58). The square of 
S 

these ratios is represented in Figure 2.5 , using the standard approximation. 

E 


They increase monotonically from 0 = 0 to 0 . As 9 9„> — cotan9„, 

_ R R E R 

E X 

^ 0, and the wave becomes longitudinal (^/k) and linearly polarized. 

A at resonance, that is the parallel component of 


For {^ « 1, 


E 


the electric field stays small compared to the perpendicular component. 


25 


SEL 73-043 



SEL 73-043 


to 

0 > 



FIGURE 2.3. POLAR DIAGRAMS OF WHISTLER REFRACTIVE INDEX FOR THREE RANGES 
OF FREQUENCY, 
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FIGURE 2.4. VARIATIONS OF N VERSUS THE ANGLE 0 BETWEEN THE WAVE NORMAL AND THE STATIC MAGNETIC 
FIELD FOR DIFFERENT VALUES OF THE PARAMETERS A AND B* Solid lines: g » 1; dashed 
lines: p - 2- 







c, MAXWEU^“VXASOV DESCRIPTION OF THE INTERACTICN 

To deal with an ensemble of particles we introduce the particle 
distribution functions f (r,p,t) where subscript s refers to the type 

g ^ ^ 

of particle ♦ We chose a momentum variable p in order to include 
relativistic effects. Assuming no collisions between particles, Liouville’s 
theorem implies that each distribution function is conserved as time 
elapses : 

^ f^(r,p,t) = 0 (2.27a) 

Equation (2.27a) when expressed in terms of partial derivatives yields 
the Vlasov equation: 

p Sf / PXB \ 5f 

In the last equation, the force term is written explicitly (Loren tz 
force law; we assume only electromagnetic forces are important). 

The solution of the self-consistent system of Maxwell and Vlasov 
equations is the solution of the wave-plasma system. The precise form 
of the Maxwell set depends upon whether the plasma is regarded as a 
collection of particles in free space, or as a dielectric with an 
equivalent permittivity. 

The Maxwell equations written below show the relationship between 
one description and the other one: 
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Description as particles in free space: 


Dielectric description: 


= - 3t 2 

IXK = 3t £ + i * iext 

V -B = 0 


V -D 


^ Pext 


^ = SqE 

B = 


{2.28a) 


VXE = - 


^ 2 


^XH = ^ D H- 

£’B = 0 
= Pext 

B = 4qH 


(2.28b) 


In the preceding sets of equations, p and J represent external 

0Xt 6 a L 

particle and current densities. We need also the charge relations: 


p = 2 
s 


J = 2 q 


cs 


f d^p 
S 


f d^p 
S / M s 
s s 


(2.29) 


(2.30) 


In the context of the linear theory we assume that the wave is a 
small perturbation which causes a perturbation of the particle distribution 
function which is small with respect to the equilibrium distribution 
function which exists in the absence of the wave: 


E = 0 + Ej^(r,t) 


B = ^ + Bj^(£>t) 


^s = ^Os^P^ + 


(2.31) 


Strictly speaking, this expression is valid only for plane waves. If 
the medium is homogeneous and stationary, the most general linear 
relationship between D and E is a convolution whose Fourier-Laplace 
transform yields this expression (cf. Quemada [1968]). 
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Assuming the medium as infinite, we can Fourier transform fields and 
distribution functions in space* Laplace transform is used to transform 
the time variable. Use of the Laplace transform introduces the physical 
principle of casuality (the perturbation is considered as beginning at a 
given time). This condition, initiated by Landau [1946], enables one to 
give a sense to some integrals otherwise indefinite. After transformation 
the kinetic equation reads: 

f^ [k,p,co] =Q [k,p,cn] *E_ [k,0)] + i.t. (2.32) 

Is ^ ^ — S — — —i 

(i.t. stands here and hereafter for initial terms involving the spatial 
Fourier transform of quantities at t = 0) . 

Square brackets are used to distinguish spectral functions from 
functions in the space-time domain, A transform pair is defined by: 

CO CO 

f [k,p,(Jo] - f e — d^r / dt (r , p, t ) (2.33a) 

CO 

f(r,p, t) = — J d^ke^^-'- J [k,p,co] (2.33b) 

(2it) -<0 *'w 

(similar transform pairs are defined for the electromagnetic field) - 

The quantity ^ = ±1 is introduced to take care of the two most common 

conventions in the literature. The contour W is a straight line parallel 

to the axis Read and lying above (for ^ = +1) or below (for Q - -1) 

the poles of the integrand in order to conserve casuality in the system. 

Q is an operator whose expression is determined by solving the 
— s 

coupled Maxwell and Vlasov equations after linearization. 

The current relation is 

oo 

/ — f_d^p = a[k,a)]’E[k,m] + i.t. (2.34) 

—1 s ^ M Is _ ^ —1 ^ 

S -03 S — 
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where the components of the conductivity tensor are given by 


a 


ab 



a b ,3 


M 


d P 


(2.35) 


where C is a mapping of the contour W. It has been introduced by 
Landau [1946] in order to facilitate the integration of Eq. (2.35) (see 
Stix [1962] for example) . 

In both descriptions (particles in free space and dielectric 
description) the first Maxwell equation reads: 


kXE^ = + i.t. (2.36) 

while the second equation, written differently in each description, 
defines the equivalent dielectric tensors K: 

- ilXli = ^ ^ ^ -1 + (. 2 . 37 ) 

c c ” 

From this equation we deduce the equivalent dielectric tensor 


K[k,o] 


CT [k,co] 

= ~ 


(2.38) 


The elimination of vector in both Eqs . (2.36) and (2.37) yields 

the wave equation: 


A [k,co] = S[k,o)] 


(2.39) 


where S takes into account all the initial terms, and 


Aok = NN, - N6-+K. 
ab a b ab ab 


(2.40) 


where the refractive index 


N = - k . 
— to — 


(2,41) 


The electric field solution is obtained from Eq . (2.39): 
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where 


[k,o>] -S [k,co] 

Ej^Ck.m] = A [k,co] -SEk^O)] = ~ p^■^^ > 

D[k,(J 0 ] = Det(A) . (2.42) 

By inverse Fourier -Lap lace transform of Eq. (2.42), we arrive at the 
solution for in the space-time coordinate system. By the residue 
theorem (assuming the components of S to be entire), the only contribution 
to the integral of the inverse transform comes from the poles of the 
integrand, that is, the zeros of the denominator. 

The equation which gives the zeros of the denominator in Eq . (2-42), 

D[k,CD] = 0 (2.43) 

is called the dispersion equation. Its solutions define the wave modes 
of the wave-plasma system. Equation (2.43) can be satisfied for complex 
values of k or cu (in fact to be consistent with previous convention 
m has to be considered in general as complex). For instance a complex 
value of CD for real k means either wave growth or damping according 
to signs of CD^ = ImcD. If ^cd^ > 0, the perturbation will grow without 
limit with time and the distribution is unstable. However, the fact that 
the perturbation can grow with time is not enough to determine the 
spatial characteristics of the growth. The perturbations are not single 
monochromatic plane waves but a superposition of them given by the 
Fourier -Lap lace integrals (2.33a, b). A single frequency may grow exponen- 
tially with time but the amplitude of the wave packet as a whole may 
remain finite at a fixed point in space. This leads to the distinction 
of two kinds of instabilities, ’’convective’^ and ’’absolute” (or nonconvec- 
tive) instabilities. In convective growth, the wave packet is amplified 
as it moves along (that is, at each point in space it first grows and then 
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decays), whereas in the second case it grows without limit in each point 
of space as t - oo. For further discussion, see for instance Akhiezer 
et al 11967]. This distinction has been discussed Sturrock 11958]. 
Rather involved criteria have been developed to distinguish between the 
two kinds of instability [Briggs, 1964; Derfler, 1967, 1970]. Yot the 
first kind of instability, the medium acts like an amplifier, whereas for 
the second case, it works like an oscillator. 

The work of Lee [1969], and Lee and Crawford [1970], though 
restricted to strictly parallel propagation, showed that in the magneto- 
sphere absolute instabilities appear to occur only under extreme condi- 
tions . Therefore, we shall assume hereafter that we deal only with 
convective instabilities. For convective instabilities it is immaterial 
whether we consider co real and k complex or vice versa. Provided 
|k, I « |k| and |cd^| « o:>, it can be shown that k^ and are 

related by the relationship 

CD = - k “V (2.44) 

i -i -g 

where v is the group velocity of the wave packet. The physical content 
of Eq. (2.44) is that CD^ describes the wave packet amplitude variation 
in the wave group velocity frame. Convenience will dictate the choice 
of the complex variable. For present purposes it will be more convenient 
to discuss gain rates per unit distance whereas in ray tracings 
parameterized in the time variable, it will be more convenient to work 
with gain rates per unit time. 

We derive in the next section the expression giving the growth rates 
of small amplitude waves propagating in a magnetoplasma composed of a 
cold plasma permeated by a tenuous energetic particle population. 
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D. 


GROWTH RATE EXPRESSIONS 


In order to derive the dispersion relation for our system it is 
first necessary to determine the tensor K. The derivation of the 
expression of K is somewhat lengthy but does not present any particular 
difficulty- It has been given in a number of textbooks [Montgomery and 
Tidman, 1964; Bekefi, 1966] and we do not repeat it here. In general 
can be expressed in the form (where we have explicitly derived the 
relativistic form) : 


where : 



P.. 'P 


Hm 


00 00 




S r / 2it dp 

m=-<» •'O •'-oo ^ 


dP„ 


p„-p 


II Rm 




(»|. -Vi 


-'’Rn.2 


(2.45) 


V(p ,p ) = 

II 1 




+ \/p^ + 


0 ■ Pi + P|| 


/ 


m 


2„2 

A N„ 


2 1 N/ 2 2, 

Pn + ^ Hp_ + p, ) 


N, 


In Eq. (2.45) the expressions for P„ t o given by Eq. (2. 10), 

Rml } 2 

provided in the expression of N,, = kc/cD, CD is considered now as 


complex. In the non-relativistic limit V(p^ -.iP.) V(p 0, 

Rml X Rm2 X 

and Pr. T P„ with _ given by Eq . (2.12) (where again o) must 

Rml RmO RmU 

be considered as complex). 

In Eq . (2.45) 

—s — s 


and 
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Id 


tnJ J' ^ 
. mm 2 

- 3 — 

B 


mj 

m 

a ^1 ^1 
s 


J = 
— s 


raJ J' ^ 
. mm2 

- Sts‘ -;— "j. 

S 


, 2 2 
P| 

m X 


CC iJ J' P P I (2*46) 
^^B m m II X ' 


mJ 

ra 

— P||Pi 

s 


rr iJ J' p p, 

m m II X 


2 2 
J p 
m II 


where T. = ±1 and J and J' are the Bessel functions of order m 
^i,e m m 

and their derivatives of argument: 


k p 

a = ^ ^ 

B M O) 

s cs 


(2,47) 


The quantity U has the definition: 


U = 
— s 


U 0 0 

s 

0 U 

0 


s 

0 u 


Is' 


^ 0 , , ^Os 

"s = P,, ) -5T * ■'' "ST ■ 

X " 




Os 


J. = IVl OJ — T 

Is s dp 


, ^ /^Os ^OS \ 

5 Pj. \ dPi, "j- dp, '’■/ 


(2.48) 


where is the equilibrium distribution function of particles of 

species s normalized to unity: 



(2,49) 


When can be represented as a relatively high density cold plasma 

background permeated by a relatively low density energetic particle 
distribution (as is the case in the inner magnetosphere), then we can 
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(2.50) 


^0 = ^Oc ^ ^OH 


[3 3 

e: dp « dp 

^OH ^ ^Oc ^ 


and use the approximation: 
u(p„ ) 


/- 

J p„ 


dp, P 

p II «/ P -p 

^Rjn - ^ 


a 

/ 


U(P„) 


dp„ + Qin^uiT?^) 


Rm 


(2.51) 


r Im p„ > 0 

^ Rm 


and its analytical continuation in the other half -plane in the Landau sense 
P means that the principal part of the integral, in the Cauchy sense, has 
to be taken and u stands for an arbitrary function of . Inserting 
the expression of given by Eq . (2.8) into Eq . (2.48) we find: 


U 


Is 


P|| 

= — • U 


Rm 


(2.52) 


Rm 


and in this case the Onsager relations hold (K = K >K = K ; the 

xz xz yz zy 

relation K = K was already fulfilled). 

' xy yz 

Condition (2.50) enables one to develop Eq. (2.43) in the form 
[Kennel, 1966]: 


D[k,CD] D^[k,0)} + D^[k,0)] = 0 


(2.53a) 


where is the cold plasma dispersion relation (given later by Eq. 

(2.57)) and |d^| « |i^qI * Since the energetic population is dilute, 
the real part of the refractive index is determined by [k , 0 )] = 0 
(see Appendix C) . Assuming real co and complex k, Eq. (2.53a) is 
rewritten as 


—r 


ik, ,cd] ^ 

D„[k ,m] + ik. 


0 — r —X 

r 

T [k ,co] 1 

+ i Im <D. [k ,< 

1 -r ; 

1 1 — r 
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which shows that: 


h ' 'sr 


(2.54) 


where In D, [k ,o)] is rewritten as Im D for short. 
1 —T J-r 


^0 /^o 


^ " 'Sir/ ^ 


(2.55) 


that is 


is colinear with v and Eq . (2.54) enables one to evaluate 


k (component of k along v ), which is the physically meaningful 
ig -i -g 

quantity since wave amplitude gain has to be evaluated along the ray path. 
It can be shown that so that, according to Eq . (2.55) 


■^i* 'Sir 


fPM' h)i 


(2.56) 


where ; 


“o = "s"' - v' " '=s 


(2.57) 


where A , B and C are the parameters introduced by Stix [1962]: 


S' S 


2 2 

A - sin 9 + P cos 0 , 

o 

2 2 
= RL sin 9 + PS (1 + cos 9) , 

O 


Cg = PRL , 

S = i (R + L) , Dg = I (R - L) , 


R = 1 - S Y ? 

s cs 


p = 1 - S Y, : 

s 


S '='S cs 


(2.58) 
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so that: 


Im D 


k. = 


Ir 


2 2A N^- B — 
S S m 


cosa 


(2.59) 


If time gain is more appropriate, it is immediately given by Eq. (2.44), 
i .e . , 


(0. = - k. V 
X ig'-gi 


(2.60) 


where 


^ d(cuN) 


c 

Cod? 

<5od 


ik ■ I le ^"<N)^e 


-0 


(2.61) 


and e and e. are the unit vectors along and across k respectively, 
-k -9 

With condition (2.50) we have: 

l^abl l^abl (2.62) 


where and are respectively the Hermitian and antiHermitian part 


of the dielectric tensor: 


2K^ = K + k'*' , 2iK^ = K - 


(2.63) 


where t denotes the Hermitian conjugate (complex conjugate of the 
transpose). With Eq . (2.63): 

2 2 2 

Im = K (S-N )(P-N sin 0) 

ir XX 

+ r(S-N^cos^9)(P-N^sin^e)-N‘^cos^9sin^0] 

yy 

+ [2d (P-N^sin^9) ]-K^ [2(S-N^)N^cos0sin0] 

zz S xz 

- 2pg U(^(P-N^sin^0) + K^^(N^cos9sin0)] 


and with approximation (2.51): 

m=a> ^ oo 2 


K = -c S 


.. 2 
2rt 


(2,64) 


m=-oo 0 x=l 
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so that: 


^ 2 

2jr Y “ 
im D = -c E S 

s II m=-co 


j *oo 2 

0 i=l ' 




2 

2 22 / \ 22 

(S-N )(P-N sin e) I 1 + 2Q D_(P-N sin 0) 

\ a / o 

V s / 


+ [(S-N^cos^0) (P-N^sin^0) - n"^cos^0] (J')^ 

1 ni 


mJ J' 
m m 

a 

s J 


^ ^ cos0sin0 

L s 

* "Li c(sVco.^exs-N^ - d|] / jdp^ 


— (S-N^) + Q D_J J' 
a S m m 


(2.66) 


Apart from the relativistic correction and the simplification (2,52) 
this expression has been derived by Kennel [1966]. Equation (2.66) has 
been further simplified by Brinca [1972] by making use of the cold plasma 
dispersion relation rewritten as: 


[(S-N^cos^6)(sV)-Dg](P-N^sin^9) = (S-N^)/sin\os^9 . (2.67) 


Using Eq. (2.52) brings still further simplification compared to 
Brinca’s expression (enabling one to use once more Eq. (2.67)) and Eq. 
(2.66) can be cast into the compact and relatively simple form; 

i 


''^PRmi'V®ms<PRmi' V^s^PRmi’Pi^Px ^P^ 


00 

Im D, = r s 

Ir 


where : 
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r = -c s ^ 


ms 


Ys 

N^-P 

IJ 


^ " 2 

m ]k „ |h 

N -S 

1 m2 m+1 

m ^ 0 , 


%S = h ’ 

Pj^ 



R 


e,i 


R } L j 


L 

e,i 


4 

L,R - 


( 2 . 68 ) 


The expression of ^*g easily obtained by inserting Ec| . (2*68) 


into Eq . (2*56) or Eq * (2.59): 

2 


2 cosa 

k = -r — 

ig ^4 COS0 


03 CD ^ ^ 2 

s o> N FH m=^“CD •'O i=l 

1/2 


F = 


2 4 2 2 2 

(RL-PS) sin 9 + 4P D cos 6 

D 


(2.69) 


It can be shown that ^^P^jni'Pj, ^ positive-definite 

quantities, therefore the integral in Eq. (2.69) is positive-definite. 
Furthermore the quantity H is positive. Therefore Eq. (2.68) facilitates 
a straightforward discussion of stability. The distribution function g^, 
being a solution of the Vlasov equation, can be expressed according to 
Liouville’s theorem in terms of invariants of the particle motion. In 
the magnetosphere the simplest choice of invariants is p (we assume no 
dc electric field) and the first adiabatic invariant (where 

^ is the dc magnetic field normalized to equatorial value). 

“^Notice the inversion between R and L according to the sign of the 
particle charge. We have defined as m = +1 the resonance at parallel 
propagation for either sign. In that case only C ^0 as N -R^ = 0 
either for normal electron interaction (N = R) or for normal ion interaction 

(N^= L) . 
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Therefore it will be more convenient to express U given by Eq . 
(2.48) in terms of derivatives with respect to p and sina: 


U = M m 
s s 


sina 


^Os 

<^P 




Os 


+ — 


p d(sina) 


I -2 m ,.| 
(-sin a + 


(2.70) 


If we now assume that the energetic component of the plasma consists only 


of electrons whose distribution function decreases monotonically with 

see Section 3B) and increases monotonically with pitch 


energy 
ang le 
s tability : 


hi ") <i 

I ^0 \ 

^ (3(sina) ^ / ' 


we reach the following conclusion concerning 


1* Landau interaction: m = 0, always damping (U < 0 always) 

2. Anomalous gyroresonance interaction: m < 0, always damping 

(U < 0 always) 

3- Normal gyroresonance interaction: m > 0, U can be > 0 or 

< 0 and there will be accordingly growth or damping; because of 
the positive term m/A in Eq. (2.70), small values of A 
will favor growth [Kennel, 1966; Liemohn, 1967]. Isotropic 
distribution (^^)/S(sina) = 0 cause only damping. 

For parallel propagation (k =0), a =0 and all Bessel functions 

s 

are null except J^, that is, only the term m = 1 is to be considered. 

Kennel [1966] found that the gain is maximum for parallel propagation for 


low frequency waves though a thorough investigation not restricted to 
only very small frequencies may yield an opposite conclusion [Brinca, 1972]. 

Within the standard approximation, Eq. (2.69) can be written in a 
simpler fonii (neglecting ions): 
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k . 
ig 




2 cosa 


2 COS0 


/ 


0. 


N m~a:> 0 C^cos 0 


m tt 

2 — r WPi' 


% “ = 0> = 2|p| , 


'=ml,2/‘^0'l<V *2^' V *2 = f 


2 • 2 
sin 9 

cos9-A)m ^ 


m 





(2,71) 


Expression (2,71) is valid in the context of the standard approxi- 
mation and for arbitrary frequency (provided the wave is not in the 
vicinity of both the LHR and the resonance cone). It is practically as 
simple as the low frequency approximation (A « 1) of Kennel [1966], 
There are two differences between Eq, (2,71) and Kennel's expression: 


1, For the cyclotron harmonics (m 0) , the coefficients ^ 

are now the sum of two terms: t^, given already by Kennel and 

a new term t^ which increases with 0. 

2, For the Landau interaction (m = 0) , is now the sum of a term 

involving J^, given by Kennel, and a new term involving J^* 

Let us test the validity of Kennel's [1966] approximations* At the 


Gendrin angle 9 , for low frequencies: 
G 


S = 


A(1-4A ) 


m 


m 


(2.72) 


and still can be considered as small compared to the first term 

t^ = 1 ± 2A ^ 1 . (2.73) 

But this is no longer true when 0 0 , where on the contrary t 

R 2 

is the dominant term and t^ can be neglected. Therefore Kennel's 
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approximation for the cyclotron harmonics is valid for 0 < 9^ but not 
for 9 > 9^- 

u 

Neglecting the first term in is more drastic because by 

expansion of J^, both terms in are of the same order (in small 

/\ or small 9). This point requires further investigation with respect 
to Kennel's [1966] conclusions (see Section 3E) . 

It is possible to give a physical interpretation of each terra in 
0^. The first terra corresponds to the "classical" Landau Interaction. 

It is a result of the action of the parallel component of the electric 
field upon the motion of the particle. The second terra corresponds to 
"transit time" Landau interaction. It is a result of the presence of a 
non-zero parallel component of the magnetic field for non-parallel 
propagation [Stix, 1962] . This component acts on the particle according 
to the equation 

,, z z 

“ “ = ■ 3 ^ 


2 y 

where |j. = Mv^/2B^ is the magnetic moment of the particle. Note that 

both effects tend to disappear when A -* 0- 

When e 0 : 

R 


Ct-C„=C =n^-, 
ml m2 m j. m 


0 


m 


I "o f 


for all m . 


(2.74) 


In the next chapter, we will obtain quantitative values of by 


choosing realistic models of distribution functions in the magnetosphere. 
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III. VARIATIONS OF GAIN RATES WITH MAGNETQSPHERIC PARAMETERS 


A. INTRODUCTION 

The variations of whistler wave amplitudes caused by wave-particle 
interaction depend upon both the cold plasma and hot plasma distributions* 
Although models of the cold plasma are fairly extensively known, mainly 
from, nose whistler measurements (cf * references of Chapter 1), energetic 
particle distributions are much less known. Both for lack of information 
and for computational convenience only simple models of hot plasma will be 
considered here. Limits on these models will be set by whistler observations. 

The purpose of the chapter is to study in detail the variation of 
wave amplitude gain with respect to various magnetospheric parameters. 

In Section B we review models of the cold plasma as deduced from 
nose whistler measurements, and a few energetic particle measurements. 

The assumption of parallel propagation has played a considerable role in 
whistler analysis, in part because the first whistler observations were 
made on the ground where only ducted whistlers could be recorded, and in 
part because of substantial simplifications in the analysis. In Section C, 
we study in detail the variation of wave gain for parallel propagation, 
assuming that the class of distributions is similar to that chosen by 
Liemohn [1967]: 

f(W,a) « sin\ (3.1) 

This model has several advantages. It is mathematically simple. The 
energy variation is reasonable, compared to particle data. As of now, 
no experimental data have confirmed the pitch angle variation, but a 
general pitch angle distribution can be decomposed as a superposition of 
such functions. 
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The model (3-1) is most likely oversimplified- It assumes separ- 
ability in energy and pitch-angle- This simplification is relaxed in 
Section D where we study a more general class of distributions that are 
nonseparable in energy and pitch-angle. This new class of distributions 
can have radically different characteristics from the previous class. 

In Section E, we make a detailed study of the variation of the 
gain with respect to arbitrary angles of propagation and frequency, and 
in Section F, we review the important contributions of the chapter. 

B. PARTICLE DISTRIBUTIONS IN THE MAGNETOSPHERE 
1 - Cold Plasma Distribution 

Inside the plasmapause, the cold plasma distribution has been 
very successfully described by the diffusive equilibrium model of 
Angerami [1966] (see also Park [1973]). 

Apart from the vicinity of the top of the ionosphere and beyond 
L > 6, the full model with three types of ions and a gravity term does 
not differ significantly from the following idealized model: 



Subscripts E and 1 refer respectively to equator and base level 
(1000 km altitude), 

n = density 

r = geocentric distance 

R^ = geocentric distance at 1000 km altitude 
R = dipole magnetic field equatorial geocentric distance 

Ej 

Hjj = hydrogen scale height = 

K = Boltzmann's constant 
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T = temperature 


g = acceleration of gravity - 

To complete the description of the model, we need the equatorial 

density profile. C. G. Park [private communication] has established an 

average equatorial profile based on whistler measurements in the month 

of June, 1965. It is given in Figure 3.1 (solid line). An equatorial 

-2 “4 

density decrease as L for L < 4 and as L for L > 4 seems to 

fit closely the actual variation. Contours of constant 3 = o> /m are 

p c 

plotted in Figure 3 . 2 , assuming an average variation ^ 

“3 

n (4) = 250 cm , to which corresponds 3„(4) 10. Using values of 

p given by the figure, one can deduce from Figures 2.1 and 2.4 the 
minimum energy of resonance for cyclotron and Landau interactions at 
each point of the plasmasphere . This is a very useful quantity to know, 
since the number of particles available at a given resonance provides an 
estimate of the importance of the resonance. Since ^ ^ it is 

a simple matter to deduce from a density profile different from 

Figure 3.2. It can be seen from Figure 3.2 that the standard approxi- 
mation is valid practically everywhere inside the plasmasphere. A few 

values of W (A =0.5, 0=0) are given in the figure. 

R1 

Outside the plasmapause , a completely satisfactory model is not 

yet available though the collisionless model of Angerami [1966] seems 

reasonable. This model is close to an r model. With such a model 

the contours of constant p are plotted on Figure 3.3 . assuming 
-4 -3 

n (L) oc L , and n (4) = 10 cm to which corresponds 6^(4) ^ 2. 

E F 

2 . Hot Plasma Distribution 

Most of the early particle observations in the magnetosphere 
were obtained with instruments (such as Geiger tubes or scintillation 
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FIGURE 3-2- 


CONTOURS OF CONSTANT g FOR A SIMPLIFIED 
ELECTRON DENSITY VARIATION oc L“^ . Energy 
(A = 0.5; e = 0). 
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numbers correspond to 
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FIGURE 3.3. A PLOT SIMILAR TO FIGURE 3.2 BUT FOR AN r MODEL ASSUMING oc L . 



counters) capable of observing electrons only above an energy of ~ 40 keV 
(see review in Hess 11968] and also Russell and Thorne 11969] and 
Vasyliunas [1969]). Now the bulk of whistler-electron resonance inter- 
action occurs within the energy range 100 eV - 40 keV , at least inside 
the plasmapause and knowledge of particles in that energy range is 
essential* However, information on the particle distribution in this 
energy range is still meager. 

As of now, to the author's knowledge, the only published data of 
differential fluxes of 100 eV - 40 keV particles are from Schield and 
Frank [1970]. We have reprinted on Figure 3.4 two of the figures from 
this paper which show fluxes observed during times following a long 
period of low magnetic activity. The features of the distribution 
functions seem fairly repeatable. 

Let us focus our attention on the plasmasphere spectrum. Above 

1 keV, the distribution function falls off smoothly like — (see 

v 

Figure 3.4). There is a flattening of the distribution in the 500 eV - 

5 

1 keV range. In the lower energy range (W < 500 eV), the spectrum 


^Confusion may arise in the term ” spectrum. Let us be clear about our 
definitions. The number of particles in a given velocity range is, by 

convention: o 

dn = f(v)d v 

where dn is the number density and f(v) is the "particle distri- 
bution function." For comparison with data, it is convenient to define 
the "particle distribution function in energy" F: 


dn F(W,0)dWdn 

where Q is solid angle. Experimentally, the current density is 

measured in terms of the quantity: 


3 

dj = evf(v)d v 

from which the differential flux d$ is obtained: 


d$ 


dJ 

edWdn 


vr(W,Q) 
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FIGURE 3 .4 . DIFFERENTIAL ELECTRON FLUXES IN DIFFERENT 
REGIONS OF THE MAGNETOSPHERE, FRCM SCHIELD 
AND FRANK 11970], 
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becomes softer with f(v) oc ~ * For further details, the reader is 

V 

referred to Schield and Frank 11970]. 


C. PARAUJSL PROPAGATION FOR A DISTRIBUTION SEPARABLE IN ENERGY AND 
PITCH ANGLE 

1 . Introduction 

Parallel propagation has played a considerable role in whistler 
stability analysis (e.g., Scarf [1962]; Tidman and Jaggi [1962]; Liemohn 
and Scarf [1962a, b, 1964]; Guthart [1964, 1965]; Liemohn [1967, 1969], 
to limit ourselves to work explicitly related to whistler gain in the 
magnetosphere). One reason for that choice is that the first observa- 
tions were made before the satellite era and in that case only ducted 
signals could be observed. Moreover the abundance of ground data proves 
that it is a very important case. Therefore the parallel propagation 
approximation is a valid approximation for all ground whistler data, 
at least as long as cold plasma propagation is concerned. A second 
reason is that the gain expression (2.69) which is rather complex 
simplifies considerably for parallel propagation. It reduces to: 


li 


2 oo 

2 03 /• 


In the standard approximation and neglecting ions Eq . (3.3) is 
rewritten following Kennel and Petschek [1966]: 


and the change of units is the following: 

d$ [(cm^XsecxsrxeV) ^] = 1.76xlO^(vF) [C.G.S.] 


Note that: 


f oc vF , 
d$ oc v^f . 


When we discuss a spectrum oc — , there may arise some confusion about 

whether we mean f, F or dl. Unless otherwise specified, we will 
always implicitly define y as the energy parameter of the ”spectrum^ ^ 
associated with the distribution function f(v). 
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< A - » V 


oo 

/ 


■'o 3(H5^ 


2 / So'’l‘‘'’l 


(3.4) 


Equation (3.4) shows that 


Aj^p(A) > ^ GROWTH , 


Arp (A) < 3^ “* damping, 


A„„(A . ) = rr-, MARGINAL STABILITY ( k. 

KP St II ig 

'^st ^ \p^^st^^ * 


(3.5) 


Therefore A must be at least positive for amplif ication# A 


is related to the sign of 


through Eq 4 (3.4) which shows that 


I 2feo \ 

isotropic distributions ^ ^^(sina) ^ yield damping (A^p = 0). 

The same conclusion is reached a fortiori for distributions with 
Ws±^) ^ ^ everywhere (A^p < 0). Only distributions with > 0 

somewhere may cause amplification • For such distributions and because 
is an increasing function of A, generally low (high) frequencies 
will be unstable (stable). 

With a simple choice of distribution g oc p "^sin^Q, expression 

Url 

(3.4) can be expressed in a simple analytic form. We are then able to 
discuss relatively easily the variations of the whistler gain rate with 
various magnetospheric parameters. This discussion is presented in 
Section 2. 
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2. The Variation of the Gain with Various Magnetospheric Parameters 


. “V . Q 

For a Distribution cc p sin g 

We represent the distribution function by (see Eq . (2.50)) 


^0 ®0c ''' ^OH’ 

, A . q 

1 xc vq v-3 sin a 

s ^„ = — h(p-p„)6^, p_ 


’OH 4rt 


®0c 4n 


^q/2 

/ ^EH , \ 6(p) 


.V 




(3.6) 


where b = magnetic field normalized to equatorial field value, 
corresponds to some convenient normalization energy (we will use 

later 100 eV unless otherwise specified); ~ > W )/n (W > 0). 

H(p) is the Heaviside (step) function introduced to limit the total 

1/2 

number of particles* The quantity p is chosen » (M KT) ^ to insure 

H U 

both that 6 « 1 (where 6„„ = n (W > W )/n (W >0)), and that the 

EH EH E — HE"" 

temperature correction need not be included in the cold plasma refractive 

index [Montgomery and Tidman, 1964]* The condition « 1 is necessary 

to validate the treatment of Chapter 2 used to find the complex root of 

the dispersion function. In these conditions it does not matter what is 

the precise functional form of , which we represent conveniently 

Oc 

“3 

by a Dirac distribution 6(p) plus a term 6^p„ to represent particles 

O H 

from thermal energies ( .1 eV) to 10 eV * This assumes a constant 

distribution in this last energy range, an assumption which may be far 
from the true physical representation. However since our results are 
insensitive to the exact form of g^^ in this energy range we feel this 
representation is as adequate as any. 

A general function of pitch angle can be expanded in Fourier 
series in terms of sin(na) and cos (na) > where n is an integer- It 
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is reasonable to assume symmetry with respect to the plane a = ?t/2 
(particle mirroring back and forth along the field line from each side 
of the equator) . Thus the distribution function is expandable in terms 
of sin(2na-a) and cos(2na) which can be developed in turn in powers 
of sina« Therefore we will consider q as an integer, while v is 
arbitrary. (In the following for mathematical convenience, we will 
derive expressions explicitly when q is even.) The normalizing condition 
requires : ^ ^^2 


A 

va 


v-3 r p dp r 

^ \ Jo 


- q + 1 -■ n 
sxn ada = 1 


or : 

A = 2(v-3) n (a+l/2)/(q/2) *. , q even. (3.7) 

vq 


The coefficient increases with both the energy parameter v 

and the pitch-angle parameter q. For a given energy spectrum, the co- 
efficient A must compensate for the increasing loss of particles at 

vq 

low pitch angle for increasing anisotropy. A similar conclusion is 
reached for the case in which the anisotropy is given and the spectrum 
softness parameter v is increased. 

From Eq - (3.4), it is immediately seen that 


KP 


(3*8) 


Arp is in this case independent of A and there is correspondingly 
one single frequency of marginal stability: 


KP 


'st 


1 + A. 


KP 


q 

2+q 


(3.9) 


the low frequencies A < A are unstable and the high frequencies 

s t 
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The coefficient B decreases with parameters v This 

vq 

expresses the fact that 

/ <» 

g p dp (3.12) 

0 0 i i 

R 

is 3 niGSsuro of the number of psrticles avsilsble at resonance. 

Clearly, for increasing anisotropy the average energy of a particle 
at resonance increases for a given minimum energy of resonance, that is 
the number of available resonant particles decreases. The same conclusion 
holds for increasing spectrum softness. Now, as opposite conclusions were 
reached regarding the coefficients A , there is a cancellation effect 

vq 

between A and B and the coefficient C stays fairly constant 
vQ vq vq 

for different values of the parameters v q- 
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The expression of somewhat complex with all parameters 

intricately mixed. There are two geometrical (space) parameters, L and 

one wave parameter A, and three hot plasma parameters # V 

£1 

q. The first three parameters influence the cold plasma propagation 
characteristics- Therefore the gain rate is an intricate combination of 
wave frequency and both cold and hot plasma parameters. 

As our main interest lies in comparing gains with data, it will be 


more useful to express gains in terms of decibels along a certain distance 

S (1000 km for instance). Let us call this quantity G : 
u u 


G = -C 8.686 k. S 160C ^ (-^ 

I u ^ II ig u ^ -3 II \ k 

L 


(3.13) 


From Eqs . (3.13) and (3.10) we can study the variations of ,.G with all 

II u 

the magnetospheric parameters . 

First of all, is directly proportional to the density of hot 

particles as expressed by the parameter 6„, • This is a characteristic of 

£1 

linear theory in contrast to nonlinear theory [Sudan and Ott, 1971; 

Helliwell and Crystal, 1973]. It is very important to emphasize the 

normalization of the density of hot electrons we have adopted and we will 

use throughout. We have chosen to keep 6 constant (ratio of number 

£ 1 

of hot electrons above 100 eV to number of cold electrons) for a given L 
shell value. This is in contrast with previous normalizations (e.g., 
Liemohn [1967]; Kennel and Thorne [1967]; Thorne [1968]) where the 
choice was to fix a normalization energy of the order of 10 keV or more. 
The reason for normalizing at such high energies is that at that time 
only fluxes above a few keV were known. With the more recent measure- 


ments of Schield and Frank [1970], fluxes down to 100 eV are known. There 


is an essential difference between these two types of normalizations which 
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is as follows: the bulk of resonant cyclotron energy is between 100 eV 

and 10 keV and therefore for the previous choice of normalization, a 
softer spectrum means usually more particles at resonance in contrast to 
our normalization. This difference is clearly shown in Figure 3.5. 

For a given value of equatorial gain is 

-3 

With a variation n oc ^ L , 

E 

.“uE**-’ “ «E1<‘'> 

For a hard spectrum (v = 4), G ^ varies like 6^-. (L) and for softer 

II uE El 

spectra, G increases with L more rapidly than 6„_ (L) . There is no 

II UJii El 

detailed information as how (L) varies but the data of Schield and 

El 

Frank [1970] suggest that the number of hot particles stays fairly constant 

across different L shells and it is reasonable to assume that either 

6„t(L) increases or at least stays constant with L. For a ducted ray 
£1 

path over a complete field line we have to multiply by a factor 

roughly proportional to L. All these factors add up to indicate that 

the gain rates are increasing with L. This hypothesis appears to be 

supported by observation. For instance Dunckel and Helliwell [1969] 

observed that inside the plasraasphere emissions are more frequent close 

to the plasmapause. Conclusions would be roughly the same for a variation 
~2 -4 

Hg cc L , L , and therefore outside the plasmapause as well as everywhere 
inside . 

Now we look at the variation of the gain with latitude keeping 

A constant: 


59 


SEX, 73-043 



SEL 73-043 



ENERGY (eV) 

FIGURE 3,5, COMPARISON BETWEEN EXPERIMENTAX AND MODEL FLUXES. The solid curve represents the 

experimental flux, and is taken from Schield and Frank [1970] (see Figure 3.4a). 2 

The solid straight lines represent fluxes for a model of distribution f(v,a) ~ v sin a 
and are normalized to n(W > 100 eV) = 1 cm“^ . The broken curves show Liemohn’s [1967] normalization 
and the cross shows Thorne’s [1968] normalization. 




(3.16) 


G (X) 


II u 




jj3-q/2^ 



G (X) depends now on q. Away from the equator the gain rate wil 

— Q / 2 

decrease because of the factor b in Eq. (3.16) for an anisotropic 

distribution. This fact is simply a consequence of the first adiabatic 
invariant law. The parameter v affects the gain variation quite 
differently inside and outside the plasmapause. Inside the plasmapause , 
n(X)/n constant and the gain decreases with increasing \ because the 
resonant energy is increasing away from the equator. Therefore the number 
of resonant particles decreases away from the equator. This effect, 
combined with the particle decrease due to the first adiabatic invariant 
law, yields a fast decrease in I G I with \ as illustrated in Figure 

II ^ 

3.6. Therefore the major contribution to gain is concentrated in a 
relatively narrow latitude range around the equator. 

Outside the plasmapause n(\)/n b and from Eq. (3.16), G (X) 

h II u 

is about v“ independent . The decrease of the gain is therefore much 
slower outside the plasmapause but notice we have not yet taken into 
account the decreasing of A away from the equator and that |„G^| 
decreases with A reasonably soft spectra (v > 5; see Figure 3.7 ). 

Because of this supplementary factor, the gain is still relatively 
concentrated around the equator. 

In order to have a complete picture of the gain variation, it 
remains to study its variation with v ^nd q as illustrated in Figure 

3.7. The values of the gain are high for a hard spectrum v = 4 but 
decrease for softer spectra. Notice it is not too high for \} - 5 and 
reasonable anisotropy (q = 2). As pointed out before the softer the spectrum 


61 


SEL 73-043 






SEL 73-043 


FIGURE 3.7. A PLOT OF PARALLEL WHISTLER GAIN VERSUS NORMALIZED FREQUENCY FOR A FEW 

VALUES OF ENERGY PARAMETER v AND PITCH-ANGDE PARAMETER q. is the 

density of particles of energy > 100 eV normalized to the total density 
in the magnetic equatorial plane. 







the less particles in the mid -frequencies (A ^ 0,5) because of the 
normalization of the distribution function we have adopted (same number 
of particles above 100 eV) . This effect clearly overtakes the effects 
of increasing: slope. The qualitative behavior of the gain remains the 
same for a given value of q (same number of marginal stability frequency) 
and different values of v because of the factor ^ ^ 

Eq . (3 .10) . 

D, PARALLEL PROPAGATION FOR A DISTRIBUTION NONSEPARABLE IN ENERGY AND 

PITCH ANGLE 

Up to now, we have assumed a very simple type of distribution. We 
wish now to discuss the stability with respect to more general types of 
distributions. One of the characteristics of the previous type of distri- 
bution is that it exhibits a single marginal stability frequency 
Moreover for higher frequencies, the gain is a very rapidly decreasing 
function of frequency. This has led several authors (Llemohn [1967]; 

Thorne [1968]) to hypothesize that whistler cutoffs in observed data are 
caused by hot plasma damping effects. Almost invariably, careful investi- 
gation has shown that cold plasma accessibility effects are a very 
plausible alternative explanation [Carpenter, 1968; Edgar, 1972]. Even 
in the amplification regime (A < A t values of the gain predicted 

by various workers look very high compared with observations. It is 
therefore interesting to investigate whether another type of distribution 
might exhibit not a single frequency marginal stability but a broad band 
marginal stability. 

We know fairly well (see Figure 3-4) the energy law of variation 
(«: p f however the pitch-angle law is unknown down to low energies (i.e., 
~ 100 eV) . Pitch-angle laws different from those considered in the last 
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section are needed to yield broadband stability ♦ A general distribution 
still separable and o: p can be expanded in a series of the variable 
sina* Therefore let us consider a distribution: 


^OH ^ 


-V 


q=0 


A sin a 
q 


For such a distribution; 


A 


KP 


q=0 



s 

q=0 


B 

vq 


(3.17) 


(3.18) 


This shows that A is independent of p , that is of A. Such a 
KP II R 

distribution would have exactly the same characteristics as that of Eq . 
(3.1) with respect to stability and we can state: 

Lemma ; a separable distribution in energy and pitch-angle with an 
energy dependence cc W ^ and an arbitrary pitch-angle dependence 
exhibits a single stability frequency . 

Identical conclusions are reached with a function 



. q 

sin a 


q ^ 0 . 


(3.19) 


Therefore, a fairly general class of distribution functions 

separable in W and a yields only one stability frequency. In order 

to look for a function which may yield marginal stability over a broad 

frequency band (such a function may be thought to have in the limit an 

infinite number of values A for which G =0), we now consider a 

St u 

nonseparable distribution in particle energy and pitch-angle. 

This is a natural choice. We have seen that an isotropic distribution 
yields damping for every frequency. The more the function is anisotropic 
(the higher the value of parameter q) the higher is the frequency A 
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Now high values of A correspond to low values of resonance energy, and 
vice-versa • It is therefore expected that a broadband marginal stability 
function is one which has a large anisotropy for low particle energy and 
low anisotropy for high particle energy* Therefore, such a function must 
have a pitch-angle anisotropy which is energy dependent* This can be 
seen equally well in a very simple manner from the interaction of a single 
particle with an electromagnetic wave propagating parallel to the static 
magnetic field* This interaction is described by the constant of the 
motion of the particle during the interaction (Eq. (2*19 a or b)) and its 
graphical interpretation (Figure 2*2)* For high resonance energy the 
interaction results predominantly in pitch-angle scattering (with virtually 
no energy exchange) whereas for low resonance energy both pitch-angle 
scattering and energy exchange occur* Therefore if waves and plasma 
constantly interact in the magnetosphere, as we have every reason to 
believe, the shape of the particle distribution may eventually be con- 
trolled by the process. This change of the particle distribution by the 
interaction with electromagnetic waves is described (within certain limits 
of validity) by the quasi-linear theory which shows practically pure 
pitch-angle diffusion for high particle energies, but for low particle 
energy, there is energy diffusion as well (e.g*, Kennel and Engelmann 
[1966])* A criterion for making the difference between pitch-angle 
diffusion only and both pitch angle and energy diffusion can be found in 
Gendrin [1968], for example. Pure pitch angle diffusion can be considered 
when: 

V « v (3 .20) 

P R 

In view of the preceding discussion, the hypothesis of separability 
in energy and pitch-angle may only be considered as valid for low 
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frequencies of interaction (as given by the criterion Eq . (3.20) and may 
be an oversimplification for higher frequencies. 

The simplest choice of a nonseparable function related to the previous 
calculations is to choose a superposition of functions oc p sin a with 
different v ^nd q: 


. q 

^ ^ a 

^OH “ ^ ^ V 

V q p 


(3.21) 


For mathematical convenience, we choose; 




. . V“3 


'^q 2 

exp[q{pj^/p) ^sin a/bl 


( 3 . 22 ) 


(By expanding the exponential in series, it is immediately seen that this 
function belongs to the class represented in Eq . (3.21). 

The parameters v, Q snd are chosen to match published particle 

data. The parameters v and q have the same physical meaning as given 
previously, i.e., v expresses the energy variation of the distribution 
function whereas q is an anisotropy factor. The new parameter 
expresses the coupling between energy and pitch-angle of the particle. 

Pi 

The normalizing condition is written with change of variable y = — 

1 «/2 ^ 

A ^ - f dy f expCqy *^sin^a3 y'^ ^sinada (3.23) 

^ •'o •'o 


and by expansion of the exponential in series; 


n 


2" n: 


V-3 


S ^-3 + nv (2n+l)I 
n=l <1 


(3.24) 


From Eq . (3.24) we obtain: 


1^) - 


cl Vo S ^ S B „»-(! ->,n3 

nr=0 q 


V+nv ^3 


(3.25) 
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Inserting Eq. (3 ,24) into Eq . (3*25) we get the gain rates for the 
distribution Eq* (3*22). 

Figure 3.8 represents G for fi = 10 (corresponding roughly to 

II u 

L = 4 at the equator) versus normalized frequency. The values of v 
and q are kept constant (4 and 5 respectively) but different values of 
are chosen to see the important effect of this newly introduced 
parameter . 

For the case = 0, which is the particular case of separability, 

the gain is very similar both qualitatively and quantitatively to that 

••VI q 

produced by the functions oc p sin However it can be seen that the 

gain decreases as v increases. For a value v =0.5, the maximum 

q q 

value of G (A) is already less than 1 db. Now a value v = 4 represents 

II u 

an upper limit in value of G and the gain is expected to be lower 

III ^ 

for higher values of v Figure 3.9 the function Eq. (3-22) has been 
represented for a value =0.5 in terms of differential fluxes versus 
energy and versus pitch angle with an appropriate normalization to closely 
fit Schield and Frank’s data [1970]. The experimental values of a were 
close to 90*^ (the field of view of the electron analyzer was directed 
earthward while the spacecraft was moving almost parallel to and near the 
magnetic equatorial plane) * The fit is better than a pure power law 
(«5c p ^) * In particular, the spectrum is softer for lower energies as 
indicated by the data. This may be a manifestation of increasing energy 
diffusion at low energies. Bogott and Mozer [1971] measured pitch-angle 
distributions on ATS-5 satellite at synchronous altitude during quiet time. 
The data show quite isotropic distributions in the measured energy range 
(40 keV and above) - Our proposed model of distribution also fits this 
observation. 
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FIGURE 3,8. PARAUJEX GAIN FOR A NONSEPARABLE DISTRIBUTION. 

Notice the Important influence of the parameter 
V of energy and pitch-angle coupling in the 
amplification regime. 6^,^ is the same parameter 
as in Figure 3.7. 
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FIGURE 3-9. ELECTRON DIFFERENTIAL FLUX FOR A NCWSEPARABLB 
DISTRIBUTION. Comparison Is made with the 
plasmasphere differential flux of Figure 3.4. 
The experimental values of pitch-angle a were 
close to 90°. 
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There may be some irrelevance in comparing Bogott and Mozer^s [1971] 
measurements which were made outside the plasmapause with measurements made 
inside the plasmapause. In favor of the comparison, from Schield and 
Frank [1970], high energy fluxes are fairly similar inside and outside 
the plasmapause. Secondly, the Kennel and Petschek [1966] mechanism 
predicts a small anisotropy for 40 keV electrons (A^p ^ ^ 0.3). 

This mechanism has recently been extended by Lyons et al [1972] to include 
arbitrary angle of propagation and provides a satisfactory explanation to 
both the ’^slot” between inner and outer Van Allen radiation belts and the 
shape of pitch angle distributions of high energy particles [D . Williams, 
private communication] . 

V 

E. NONPARALUEL PROPAGATION 

The purpose of this section is to compute the variation of the gain 

with the wave normal angle. The general expression of the gain is much 

more complex than for the parallel propagation case. For simplicity we 

“V . Q 

study first such simple separable distributions as p sin a. Then we 
study nonseparable distributions for small wave normal angle. In the 
latter case it is then possible to make an expansion around 9=0 of 
the expression (2.69), limiting the mathematical complexity. 

1 . Separable Distribution 

A Fortran program has been developed to compute the nonparallel 

gain rates for separable distributions with integer values of v and 

even integer values of q. The details are given in Appendix A, 

The gain rates have been computed for many combinations of the 

parameters v ^ A varying from 0.1 to 0.9 and 6 from 0^ 

to a fraction of a degree from the resonance cone angle 0 - The values 

R 
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of the ratio 


are displayed on Figures 3.10 and 3 >11 for a few 


G (9) 
u 

G (0) 
u 

values of v and q* The gain decreases with values of 9 because of 
increasing Landau damping. There is a distinction to be made here 
whether the gain was negative (damping) or positive (growth) for parallel 
propagation at a given value of A- As the absolute value of the gain 
often decreases for increasing wave-normal angle, it may well be that 
there is less attenuation, but there is never more amplification at 
0 0 than at 0=0 for A < 0.5. This casts a doubt upon Brinca's 

[1972] argument that the presence of minimal gain at 0 = 0 plays an 
important role in triggering emissions, at least with this kind of distri- 
bution- In fact the cases discussed by Brinca [1972] are cases for 
which there is only damping at 9=0. For A > 0.5, maximum growth at 
Q ^ 0 is possible but this effect is very small ( (G^(0 )/g^(O) always 
stays ^ 1). The frequency A = 0^5 is a transition frequency for several 
factors whose sense of variation with 0 is different for values of A 
above or below 0.5: and with 0 when 

A < 0.5. At the same time there are more Landau particles than fundamental 
cyclotron particles for A < 0.5. Therefore, the first cyclotron resonance 
part of the gain can only decrease (in absolute value) for 0^0 (and 
8 < 0^) and A < 0.5 . 

We have included in the program harmonics up to order |m|< 3. The 
conclusion is close to Brinca *s [1972] conclusions: harmonics of order 

m ^ 0,1 (Landau and fundamental cyclotron interaction) are almost always 

negligible. For A V^q ^ Pjj ^2 spectrum is hard enough so 

that the number of particles at Landau and second harmonics are not too 
different, the strength of the normal second harmonics interaction and the 
Landau interaction can be comparable. For A 0, p P » and the 
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FIGURE 3.11. A PLOT SIMILAR TO FIGURE 3.10 BUT FOR TWO 
DIFFERENT VALUES OF PITCH -ANGLE PARAMETER. 
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harmonic m = -1 can be non^negligible compared to the harmonic m 1 , 

It is seen from Figure 3*10 and Figure 3-11 that there is a range of 

values 6 < 9 (G (0 ) = 0) for which the distribution was unstable. 0 
s u s s 

is decreasing with increasing value of v» was already shown by Kennel 

[1966]- Therefore Kennel conclusions concerning an unstable case of 

angle < 0 are valid, though his approximation for Landau interaction is 
s 

incorrect (see Section 2D)- For A < 0-5, there are more Landau than 
cyclotron particles and the number of particles being weighted by the 
steepness of the energy spectrum, the ratio of available Landau particles 
to available cyclotron particles is increasing with v This ratio is 
increasing when A decreases- It is expected that 0^ is correspondingly 
decreasing. However, the relative number of particles available at 
resonance is weighted by the factor 0^- As pointed out in Section 2D, 

0Q decreases when A decreases . This has an opposite effect to the 
increasing number of Landau particles when A decreases. For the ’Vlassical'^ 
part of the Landau interaction, this just expresses the fact that the 
parallel component of the electric field decreases relatively to the 
perpendicular component when A decreases (see Figure 2.5). For a given 
value of \)f 6g is increasing with q. The increasing of anisotropy 

affects differently cyclotron and Landau interaction. In Eq. (2-70), the 

^0 2 1 

derivative ^ is weighted by a bigger factor (-sin a + t- ) for 

o s 1 } j\ 

2 

fundamental cyclotron interaction than for Landau interaction (-sin a) • 

Close to the resonance cone , there is a rapid change of topology in 

the refractive index around the Gendrin angle 0_ (see Figure 2,4) which 

G 

is reflected in the gain of Figures 3,10 and 3,11 where vertical arrows 
show the location of the Gendrin angle. The gain in absolute value even- 
tually becomes much bigger than for parallel propagation as 9 0 

R 
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because p -^0 and there is an increasing number of particles at 

resonance • With the values of v q chosen on Figures 3.10 and 

3.11, there is only damping as 0 0^^. 

For A « 1, p„. \ Pt. > anisotropic, 

RV“m) Rm O 


U 

m 


„ 1 » ^0 
p A d(sina) 


(3.26) 


Furthermore from Eq . (2.74), ® -* 0 when 0 -► 0„. Therefore the 

-m m R 

cyclotron harmonics tend to cancel each other. 

In the general case, using Eq. (2.74) for 6 0 and the functions 

R 

B™^(d) defined in Appendix A, we find that as 0-^0: 

vq R 


^ \^Z\ I 

m '■ 


roq ginm 
\h-2 , q A v+2 , q-2 


(3.27) 


_mm 


where the functions B are evaluated at d = (A-m)tan0^ 

vq m ‘ 


and the 

coefficient of proportionality is positive. For high values of A, we 


make a rough analysis treating d = Atan0 and d = (,!-[{) tanQ as 

0 R 1 R 


sma 


11 quantities and assuming that we can neglect harmonics of order ^ 0,1 > 


Growth will occur at the resonance cone if 

,2 

I ^ 

-(q+v) 


- — B 

Vf2,q' 2 v,q+2 


^1 ^ 

+ c r B.^-(q+v)B ^ (t^) > 0 (3.38) 


q 

A ^vq 


V,q+2 4 1-A 


(and damping in the opposite case). 

The coefficients B are given in Eq . (3.11) which shows that: 

vq 


V,q+2 q+v vq’ v+2,q \jq q+v 


(3.29) 


and Eq, (3.29) is rewritten as: 

(q+2) ^ -(v-2) + 


^ - (q+2) 


^ (^)'"> 0 
4 ^ 
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or 


q > 


" " 4 tv-2-d^] 

Jl 

1 - 1 + ^ 

A + 2 2 ^ A 


(3.30) 


For parallel propagation, the corresponding inequality was 

q > 2/(i “1) (see Eq. (3.9)). For example for \) - 5 and A = 0-6, 

A 

Eq . (3.30) requires q > ^ 6 for growth at the resonance cone whereas 
only q > 3 was necessary for growth at parallel propagation. 

This rough analysis indicates that distributions need to be signifi- 
cantly more anisotropic to yield growth for propagation close to the 
resonance cone than for parallel propagation. 

2 . Nonseparable Distribution 

By expanding Eq. (2.71) around 0=0, the variations of the 
gain for small angle can be found: 


2 cosa 




lpl / 

(2 

A (1 ^RO 




1-A Pj^ 

" 2 Po " H 


„ ^0 


- 1)] - 1- A" (y „“i 

1 1 Su 

l-A) J 




^ 2(1-A) J 






(3.31) 


5u 


where U = U( p„ ,P«)i and the quantity — is to be evaluated at 
II m II Km -L op 


*^ii II ^R1 ‘ 

Inserting Eq . (3.22) into Eq. (3.31), we find the variation of 
the gain for a nonseparable distribution and small angle. (Numerical 
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computations are made with a program developed in Appendix A-) This time 
the 9 range of instability is fairly small. The type of nonseparable 
distribution we have chosen gives only damping for the m = 0 resonance 
whereas it gives an almost marginally stable cyclotron gain over a wide 
range of frequencies- The cyclotron gain does not change much with 0 
small (its stays practically marginally stable for small g) and therefore 
the Landau interaction overtakes the cyclotron interaction very rapidly 
when 0 increases. This effect is more pronounced when the coefficient 
of coupling \) increases as shown on Figure 3.12 where a value of 

q 

V = 0.25 still gives growth at 0 = 10^ whereas for v - 0.50 there 
q q 

is already damping at 9 = 10^. 

F. CONCLUSION 

A new model of distribution function has been proposed. The new 
characteristic of the distribution is that it is nonseparable in particle 
energy and pitch angle, or in other terms, particles with different 
energies have different pitch-angle distributions. It has a very important 
effect on wave amplitude gains caused by wave-particle interaction, namely, 
this distribution produces a gain which is almost stable for a wide range 
of frequencies in contrast to previous separable distributions such as 
Liemohn's [1967] which yield large positive gains (amplification) for low 
frequencies and large negative gains (attenuation) for high frequencies. 
This offers an explanation for the accessibility/instability controversy 
about the high-frequency cutoff of nose whistlers presented in Chapter 1. 

As of now, electron energy spectrums have been measured over the range 
(100 eV - 40 keV) in which the bulk of interactions takes place between 
whistlers and electrons in the magnetosphere (see Schield and Frank^s 
[1970] data). On the other hand, pitch angle information is not yet 
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available in this energy range • Therefore immediate confirmation of our 
model is not possible* We feel it is reasonable though, because, 1) it 
predicts an increase in the slope of the energy spectrum in the low particle 
energy range, in agreement with Schield and Frank’s 11970] data; 2) it 
predicts almost isotropic fluxes in the high particle energy range, in 
agreement with Bogott and Mozer’s [1971] measurements and Kennel and 
Petschek’s [1966] predictions. 

A detailed study of the influence of the angle 0 between wave vector 

and static magnetic field upon the gain rates has been presented. For 

separable distributions, the study has been made for the complete range 

0 < 9 < 9 . We have emphasized the importance of the weighting factor 
R 

of the Landau interaction which weights the number of Landau particles. 

This factor and the number of resonant Landau particles vary in opposite 
phase with frequency. We have shown that significantly more anisotropic 
distributions are required for growth at the resonance cone than for 
parallel propgation. For nonseparable distributions, the gain is evaluated 
for small values of 0. It is found that the cyclotron part of the gain 
stays small and is rapidly overtaken by Landau damping when 0 increases. 
Our new model of distribution may be over simplified in the sense that it 
was derived by looking for a wave-plasma equilibrium over a large frequency 
band but only considering cyclotron interaction. A suggestion for a dis- 
tribution model giving equilibrium also in the presence of the Landau 
interaction will be given in the next chapter - 

This chapter was concerned with the evaluation of gain over a narrow 
region in space considered as homogeneous. To compare with actual data 
the gain must be integrated over a ray path. It is the object of the 
next chapter. 
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IV. INTEGRATED GAIN RATES AND CCMPARISON WITH DATA 


A . INTRODUCTION 

The variation of the wave amplitude gain per 1000 km has been 
studied in the previous chapter. The influence of the numerous parameters 
of concern have been investigated in detail. To compare with actual data, 
it remains to integrate the gain over a complete whistler wave ray path. 

To carry out the numerical computations we use the ray tracing program 
developed by Walter [1969] and refined by Angerami [1970] . Certain 
parameter values of the ray tracing program are used as input numbers to 
our computer program. This procedure is explained in more detail in 
Appendix A . 

In Section B we use our program, including arbitrary angle of propa- 
gation 9 for separable distributions and small values of 6 for non- 
separable distributions, to test the validity of the parallel propagation 
approximation (which assumes 0=0) for ducted whistlers. 

The observations of Angerami [1970] were very important in determining 
the nature of nose-whistler upper cutoffs. Taking the same values of ray 
tracing parameters as Angerami, we are able to test the validity of our 
model of distribution against models of separable distributions. This 
investigation is carried out in Section C. 

Recently, Dunckel and Helliwell [1973] observed signals with very 
high normalized frequency (above A = 0.9) and explained by ray tracing 
the accessibility characteristics of these signals. We investigate the 
instability characteristics of these signals in Section D. 

The magnetospherically reflected (MR) whistlers were first observed 
and explained in terras of accessibility by Smith and Angerami [1968]. 
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Edgar 11972] made an extensive study of this type of whistler. Using 
the results of some ray tracings of Edgar 11972], we study the instability 
characteristics of MR whistlers in Section E and compare our calculations 
to similar computations made by Thorne [1968]. 

In Section F, we siammarize the most important contributions of the 
chapter, and discuss their implications. In particular we note the 
relevance of our results to plasma injection experiments, and to plasma 
diagnostic techniques. 

B. DUCTED WHISTLERS 

The gain has been integrated along a realistic model of a whistler 
duct for a few distributions. The duct is described mathematically by 
the following expression [Angerami, 1970]; 

( 1 

n(L) = n ' ^ 1 + C exp 5— / (4.1) 

( 2-AL^ J 

where is the background density, C is the relative enhancement 

(depletion) of density at the center of the duct and AL represents 

the duct semi-thickness - 

We have chosen the following values: 

C = 0.2 



AL = 0*02 

n^(X=0) = 340 cm ^ 

n^dOOO km) = 3130 cm'^ (T = 1000 °K, 80% 0+ and 20% 

at 1000 km altitude)* 

The results are presented in Table 4.1 for a few separable distri- 
butions and a nonseparable one, where we give both exact and approximate 
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TABLE 4*1. GAIN FOR A DUCTED WHISTLER. GT and GP are 

respectively the total and parallel approximation 
gains (in decibels), integrated between conjugate 
points at 1000 km altitude along ray paths in a 
duct centered at L = 4 (see text for the other 
duct parameters)- The rays were started at vertical 
incidence at 1000 km altitude and at center of the 
duct, with the exception of the 6.65 kHz ray. 
is the normalized frequency at the equator. Three 
separable distributions are chosen and one non- 
separable distribution whose parameter values are 
determined in Section 4C * 


f (kHz) 


V 

q = 

5 

0 

5 

2 

5 

4 

■ 

5 

4.8 

V =0.4 

q 

1.36 

0.1 


-2.08 

+10.3 

+ 13.5 

1 

+0.56 



B 

-2.08 

+ 10.4 

+ 13.6 


+0.62 

2.72 

0.2 

GT = 

-8.94 

+ 17.1 

+24.4 


+ 1.32 



GP = 

-8.82 

+ 17.2 

+24.6 


+ 1.48 

4.08 

0.3 

GT 

-24.2 

+22.2 

+36.6 


+2 .52 



GP = 

-24.0 

+22.4 

+36.8 


+2.76 

5 .44 

0.4 

GT 

-54.6 

+22.8 

+47.4 


+3 . 92 



GP = 

-54.4 

+23.0 

+47.8 


+4.38 

6.65 

0.483 

GT = 

-115.4 

+ 14.5 

1 +57.0 


+2.48 



GP = 

-114.8 

+ 16.4 

+60 • 0 


+6 .62 


gain values. Our results confirm the validity of the parallel approxi- 
mation used by Liemohn [1967] in the case of ducted whistlers. In almost 
all cases the difference between the full treatment and the parallel 
approximation is negligible. There are two reasons for that. One reason 
is that the angle 9 stays small when a wave propagates in a duct. The 
second is that we deal with low frequencies (A < 0.5) and though there are 
more Landau particles than cyclotron particles, the parallel component 
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of the E field stays very small. The only big discrepancy is for the 
nonseparable distribution and the frequency f = 6.65 kHz, to which 
corresponds a normalized frequency A = 0.483. This frequency is very 
close to the limit frequency of duct trapping. The choice of input 
parameters is very critical in order to trap this wave and many attempts 
had to be made to find the correct values necessary for trapping, in 
contrast to lower frequencies for which a choice of vertical incidence 
at 1000 km at the middle of the duct was sufficient for trapping. It 
turned out that the variations of 0 for f =6.65 kHz were significantly 
higher than for lower frequencies. 

As expected from the results of Chapter 3, the gain is considerably 
lower for the nonseparable distributions (only a few db) than for the 
separable distributions. Therefore, unless pitch angle information is 
available, it is impossible to conclude that nose-whistler upper cutoffs 
are caused by hot plasma effects. Therefore the use of nose whistlers as 
a diagnostic tool for monitoring electrons may be limited. Our nonseparable 
distribution is constructed to imply a certain equilibrium between waves 
and plasma (neglecting sources and sinks) and in quiet times, nose 
whistlers may tell us only that indeed this situation is reached. In 
quiet times, it may be possible to correlate a noted change in electron 
distribution with a noted change in nose whistler amplitude characteristics, 
especially if there is a tendency towards isotropy. The observation on 
satellites of much higher cutoffs is very important because! 1) it is 
evidence of the accessibility explanation of nose whistler upper cutoff 
observed on the ground; 2) it sets up some limits on the values of the 
parameter's of the hot plasma distribution. An investigation of this high 
cutoff case is given in the next section. 


SEL 73-043 


84 



C. DUCT LEAKAGES 

The observations of Angerami [1970] provide an almost unique source 
of information for correlating the whistler amplitude spectrum with the 
distribution of energetic electrons within the plasmapause since his data 
and Schield and Frank’s [1970] data are from exactly the same period. In 
particular the spectra of Figure 3.4 were measured on June 23 , 1966. 

During the period between June 15 and June 23, the magnetic activity was 
quiet; index 1,2; the DST index was low and the plasmapause was 

located beyond 5.5. It may be therefore reasonable to assume there 
was no major change in the distribution function during that period. 

Briefly, the data recorded on OGO 3 on June 15, 1966, between L = 4.1 

and 4.7, showed whistler frequencies far above the cutoff frequencies 

observed on ground data. Sometimes these high frequencies were still 

ducted when they reached the satellite (see Figure 4.1a ). More often, 

they were already unducted ( Figure 4.1b ) . Careful ray tracing explained 

the propagation of these waves (see Figures 4.2 and 4.3 ) . For further 

details, the reader is referred to Angerami [1970]. The suggestions of 

/ 

Liemohn [1967] that the high frequency whistler cutoff on ground data at 
a normalized frequency of <^0.5 is caused by hot plasma effects is there- 
fore questionable, particularly because without ducting we expect increased 
damping with the class of distribution functions previously considered. 

To look at that question in detail, we chose the ray tracing 
parameters that Angerami used to explain the duct leakages . We integrated 
the gain rates on these ray paths and present the results in Table 4.2 
(see Table 4,2 caption for all details). None of the distributions are 
compatible with the data: either there will be too much damping or there 

will be too much amplitude difference between two different frequency 
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FIGURE 4.1. TIME FREQUENCY SPECTROGRAMS OF WHISTLERS RECEIVED ON OGO 3. Detailed 
explanation of the records is found in Angerami [1970]. (a) shows a 

whistler propagating in a duct up to 3/4 of the local gyrofrequency 
(see Figure 4.2). (b) high-frequency leakages from whistler ducts 

(see Angerami [1970] and Figure 4.3). 



FIGURE 4.2. RAY TRACING TO EXPLAIN WHISTLER UPPER CUTOFF 
FREQUENCY OF FIGURE 4.1a. The values of the 
ray tracing parameters at different points 
along the path are the following: 


A: 

L = 4.7407 

X = 19.48 

9 = 6.00 

A = 0.500 

W ^ - 404 
RO 

W = 404 

R1 

B: 

4.7467 

15 .42 

5.21 

0.601 

281 

124 

S: 

4.7406 

9.32 

-7 .45 

0.731 

158 

21.5 

C: 

4.7102 

7 .43 

-15.60 

0.748 

132 

15 .8 

D: 

4.4696 

2 .20 

-41.33 

0.884 

00 

16 .7 


L and A are the geomagnetic dipole coordinates. The angle 0 
(in degrees) is measured clockwise from the local magnetic field. 
A is the wave frequency normalized to the local magnetic field. 
W ^ and W ^ (in eV) are the minimum energies of resonance of 
Landau and fundamental cyclotron resonances respectively. 
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FIGURE 4.3. RAY TRACING TO EXPLAIN UPPER CUTOFF FREQUENCY 


OF LEAKAGE FROW DUCT 2, AS OBSERVED ON FIGURE 
4.1b (L component). 


A: 

L = 4.6110 

X = 24.66 

e =3.50 

A = 0.492 

*R0 - 

B: 

4.6011 

16.36 

-8.54 

0.751 

246 

C: 

4.3428 

10.98 

-33 .83 

0.766 

104 

S: 

4 .0747 

6.45 

-43 . 19 

0.703 

57.8 
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TABLE 4,2. GAIN FOR SEPARABLE DISTRIBUTIONS. The gain is 
expressed in decibels. Cases I and II refer 
respectively to Figures 4.1a, b. The frequencies 
shown are the following: Case I: f = 6.1 kHz 

is the highest observed frequency, and f - 4.17 
kHz is the frequency which value normalized to the 
local magnetic field is 0.5 at the satellite location; 
Case II: f = 8.7 kHz and f 7.7 kHz are respectively 

the highest and lowest observed frequency. The values 
of the parameters used in ray tracing of 6-1, 8.7 and 
7.7 kHz components are taken from Angerami [1970]. 

The variables v ^ energy and pitch angle 

parameters. The functions GT and GL are respectively 
the total and the Landau gains expressed in db, inte- 
grated between A and S (see Figures 4.2 and 4.3). 

For the unducted paths the values of the parallel 
approximation gain GP is also given. The normalization 
of the distribution function for the computations is 
5 = 4x10”^, to which corresponds n(W > 100 eV) 1cm 

is chosen lower than the lowest resonance energy 
a?ong the paths. 


Case 

f (kHz) 

v= 

4 

4 

4 

4 

6 

6 

6 




0 

2 

4 







GT= 

-3 94 

-117 







n 

GL=: 

-3 .2 

-1.5 

-2.34 




Hi 

i 

4,17 

— 

GTr: 

-117 

+ 15.4 

+65.2 

+86 . 6 

-58.6 

2.7 

+21.1 



GL= 

-0.75 

-0.37 

-0.5 

-0.6 

-1.2 

-1.8 

-1.2 

■ 

8.7 

GT= 

-1830 

-754 

-296 

-67.4 

-18000 

-5170 

-1450 

■ 


GL= 

-525 

-296 

-188 

-127 

-2110 

-851 

-432 

m 


GP= 

-2210 

-891 

-287 

+40.6 

-11000 

-3070 

-814 


mm 

GT= 

-871 

-300 

-70.8 

+32 .8 

-3190 

-800 

-141 



GL= 

-316 

-172 

-106 

-63.6 

-821 

-323 

-160 

■ 

■ 

GP=: 

-1100 

-348 

-817 

+171 

-2260 

-497 

-19.5 


components. (On the average, a difference of 20 db between two frequencies 
would be readily observable on the records.) 

Therefore we rule out separable distributions such as Liemohn's [1967] 
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as a possibility- It should be pointed out that for the unducted paths, 
computations with the precise values of 0 had to be done, in contrast to 
ducted propagation for which the parallel propagation approximation is 
valid- We have written a program which includes 7 harmonics (see Appendix 
A), however a reasonable approximation could be made by keeping only the 
m = 0 (Landau) and m = 1 (fundamental cyclotron) harmonics and neglecting 
all other cyclotron harmonics. 

We have repeated the duct leakage calculations using a nonseparable 
distribution of the class developed in the last chapter, and present the 
results in Table 4.3 . From the ducted path of Figure 4.2, we have 
determined, for two fixed values of v what value of the 

anisotropy parameter q which would yield small gain between points A 
and S and small differences in gains between any two frequencies. (As 
an example, we have chosen 6.1 kHz and 4.1 kHz which correspond respectively 
to the highest frequency observed in Figure 4.1a and the frequency to which 
corresponds A = 0.5 at the satellite.) A value of q between 4.8 and 
5.0 seems adequate. Note though the high sensitivity on the value of q. 
(There is a big difference in gain values with as little a change as q 
varying from 4.5 to 5.) The reason for the choice of parameter values is 
the following: v is chosen to match the measured energy spectrum and we 
have seen (Chapter 3) that ^ is a reasonable value- The parameter 

q is chosen so that the highest observed frequency is not damped. From 
Section 4 b, the parallel propagation approximation can be considered as 
approximately correct. From Figure 3.8, we see that the highest unstable 
frequency is less than the marginal stability frequency A of the distri- 

S X 

bution having the same values of v snd q but a value of v “0. In 
the case of v = Qnd v = can be obtained analytically; 
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TABIE 4,3. GAIN FOR NONSEPARABLE DISTRIBUTIONS. A table 

similar to 4.2 but for nonseparable distributions. 
V is the coefficient of coupling between energy 
aSd pitch angle. GC is the fundamental cyclotron 
part of the gain. 4x10"^, and = 20 eV. 


Case 

f (kHz) 


5 

4 .5 
0.4 

5 

4.8 

0.4 

5 

5 

0,4 

I 

6.1 

GT= 

GL= 

-18.9 

-1.97 

-1.87 

-1,80 

+8.0 

-1.7 

4.17 

GT= 

1 

1 GL= 

i 

+4.1 

-0.34 

+4,5 

-0.31 

+4.8 

1 -0.29 

II 

8^7 

oc= 

i 

i 

+ 12.3 

+ 116 

+ 181 

7.7 

S GC= 

+71.9 

+96.6 

+ 123 


^st = ' - 


1 -<l/b 

1-e ' 

q/b 


(4,2) 


For high values of q/b: 


St 


1 - b/q 


(4.3) 


We therefore expect a value 


q > 




(4.4) 


where = 0.73 corresponds to the highest frequency observed on 

Figure 4.1a, and b^ is the ratio B /b at the satellite position. 


Close to the equator: 


b ~ 1 + !• X^(rd) . 


(4.5) 


At the satellite, ~ 9*^ and b ~ 1.11, and we deduce from Eq. (4.4) 

s s 
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that q > 4 .1. 

The quantity is chosen so that there is small difference in 

gain between different frequencies- From Figure 3-8, such a difference 
diminishes when v increases, and we can set up a lower bound on v > 
which depends only weakly upon q- 

As it can be seen \yy looking at the values of k, from point to 

IS 

point along the ducted ray path, the instability characteristic of a 
nonseparable distribution is more complex than for a separable function- 
For a fixed wave frequency, a nonseparable distribution function can be 
unstable at equator and be stable off equator, along the same field line, 
in contrast to a separable distribution. Physically, this comes about 
due to the fact that the resonant energy is increasing going from equator 
poleward and for a nonseparable distribution, the anisotropy of the 
function at resonance is decreasing towards isotropy and correspondingly 
the distribution becomes less unstable. 

For the unducted path of Figure 4.3, only the fundamental cyclotron 

part of the gain could be determined because we developed a program only 

for small values of parameter d = !(A"-m)tan6l (see Appendix A)- Note 

m ‘ 

the high sensitivity of parameter q on the cyclotron contribution. As 
pointed out earlier, the instability behavior of a nonseparable distri- 
bution may be complex- Also, for the distributions we chose, there is 
both damping and growth along the path. Finally, because of unducting, 
the relative range of resonant energy is large and gain rates change 
markedly towards the resonance cone. 

The Landau contribution for the unducted path can be estimated to 
be of the same order for the class of separable and nonseparable distri- 
butions we have chosen. Both derivative terms with respect to energy 
SEL 73-043 92 



and pitch angle add up to contribute Landau damping {see Eq* (2 ,70)), in 
contrast to cyclotron interaction where they compete with each other. 

The essential difference between separable and nonseparable distributions 
was precisely that these derivative terms almost cancel each other for a 
large frequency band in the case of nonseparable distributions and do 
not in the other case. Therefore we expect from Tables 4.2 and 4,3 that 
the Landau interaction will dominate for our model. The nonseparable 
model we proposed is certainly not adequate for the unducted path for 
several reasons: 

1. We looked at wave-plasma equilibrium only with respect to 
cyclotron interaction. In the unducted case, a distribution at equilibrium 
with both Landau and cyclotron interactions should be looked for (this 
becomes especially important close to the resonance cone where the Landau 
interaction eventually dominates), for separable and nonseparable distri- 
butions . 

2. There is complete lack of information about particle distributions 

in the energy range below 100 eV and the whistler of Figure 4,1b inter- 
acts with electrons well below this energy. If used below 100 eV our 
model eventually violates the condition which is necessary for 

the validity of the gain expressions. 

We may suggest a type of distribution which is at equilibrium with 

respect to Landau interaction. Looking back at Figure 2.4, we see that 

N is fairly constant over a wide range of frequencies and angles 0 

(but a few degrees from 9 ) and so is the minimum energy of resonance 

R 

of Landau interaction from Eq. (2.12), 

A peak in the distribution function or at least a flattening in the 
vicinity of the average Landau energy of resonance would reduce Landau 
damping. To a minimum Landau resonance energy = 150 eV, there 

corresponds an average energy of interaction 

< w > = + < tan^a >] • (4.6) 


1 


/ 
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For V = 5, V =0.4 and q = 5, < W > = 750 eV. We see indeed such 

q 

a flattening around 750 eV on experimental spectra (see Figures 3.4a^b). 

The flattening of the spectrum being repeatable in two events separable 
by 20 days, it seems more than fortuitous and may well be related to 
Landau interaction. 

Several important results have been established in this section; 

1. Separable distributions such as Lieraohn^s [1967] cannot account 
for Angerami’s [1970] observations. 

2. The class of nonseparable distributions we proposed in Chapter 
3 is compatible with the high frequency cutoff at A ^ 0.75 of Figure 
4.1a. In turn the value of this cutoff determines a lower bound for one 
parameter (q). A lower bound on the coupling parameter v is determined 
by the condition that the difference in gain should be small for different 
frequencies. 

3. Our model is probably not adequate for very low energy electrons 
(a few tens of eV) because if it is accurate the unducted whistler 
components of Figure 4.1b should have been absorbed. This suggests that 
the observed flattening of the distribution in the 500 eV range reduces 
Landau damping . 

The nonseparable distribution we propose incorporates the idea of 
equilibrium between plasma and waves in the magnetosphere. This equilib- 
rium was investigated by Kennel and Petscheck [1966] and Lyons et al 
[1972] for high energy particles. We suggest here that waves may control 
efficiently the distribution of low energy electrons as well. Of course 
the picture of the mechanism should be completed by taking into account 
also sources and sinks of particles. 

Considerable attention has been focused recently on plasma Injection 
experiments [Brice, 1970; Cornwall, 1972]. By injecting cold plasma in 
the magnetosphere the resonance energy decreases and it is expected that 
wave-particle interactions become stronger as more particles will be 
available at resonance. The fact that, in certain (quiet) times, the 
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hot particle distribution is quite stable to electromagnetic waves casts 
a doubt on the efficiency of such an experiment. 

We have also demonstrated here that our method to study wave ampli- 
tude variations caused by wave— par tide interaction provides valuable 
information concerning distribution functions • Provided we know the 
flux of particles (which fixes the number of hot particles and one 
parameter, namely, the energy parameter \j) , our method determines very 
sensitively the anisotropy of the distribution function. As of now, the 
number of available data is scarce, but a controlled experiment looks 
very promising with the advent of ground based vlf transmitters such as 
the transmitter of the Stanford VLF group, located at Siple, Antarctica 
(75.55*^8, 83.55^W). Waves of frequency comparable to the minimum electron 
gyrofrequency of the Siple field line should be transmitted. The observa- 
tion of the amplitude characteristics of these waves onboard a satellite 
orbiting close to the equatorial plane would provide particle distribution 
information along the way described in this section. 

D. WHISTLERS OBSERVED AT FREQUENCIES NEAR THE LOCAL ELECTRON 

GYROFREQUENCY 

Dunckel and Helliwell [1973] observed signals at frequencies close 
to the local electron gyrofrequency, e.g., as close as ^ ^ 0.9 onboard 
OGO-1 satellite. They successfully interpreted with ray tracing how the 
wave could propagate. Since A is high, these signals must interact 
with very low energy electrons. The cyclotron resonance energy tends 
toward zero when h ^ I and secondly for A > 0.5, the wave becomes 
eventually unducted and propagates close to the resonance cone- This 
lowers the energy of resonance again. 
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Choosing the same values of ray tracing parameters they used, we 
investigated whether the distribution used in the preceding section is 
compatible with these observations or if we should have a different model. 

The caption of Figure 4.4 explains the data and shows the value of 
parameters and Table 4.4 shows the corresponding values of integrated 
gain from 1000 km altitude up to satellite location. 


TABLE 4.4. GAIN BETWEEN 1000 KM ALTITUDE AND SATELLITE. 

The gain is expressed in decibels, f = 80 kHz: 

The ray starts at 55° latitude, vertical incidence. 
In the ray tracing, the plasmapause is located at 
L = 3.9 (see Dunckel and Helliwell [1973]). The 
parameters of the ray tracing are; 



L 


q(0) 

A 

W 

RO(keV) 

w 

Rl(eV) 

m 

3 .6766 


-13.66 

0.883 

4^46 

78.1 

1 

3 .2894 


-18.42 

0.916 

1.65 

13.7 


6,:., = 4x10 ^ and n(W > 100 eV) ~ 1 cm"^; W = 20 eV. 
EX n. 


PATH 

■ 

5 

4.5 

0.4 

1 

5 

4.8 

0.4 

5 

5 

0.4 


5 

0 

5 

2 

5 

4 

AS 

GT= 

-7.61 

-6.12 

-5.28 


-128 

-11,3 

-1.17 


GL= 

-0.43 

-0.36 

-0.32 


-1.34 

o 

to 

-0.015 


GP= 

-6.27 

-5 .02 

-4.33 


-124 

-10.9 

-1.14 

AB 

GT= 

-72.6 

-59.4 

-51.9 


-1100 

-140 

-21.6 


GL= 

-4.62 

-3.93 

-3 .52 


-14.3 

-1.95 

-0.32 


GP= 

-44.4 

-36.0 

-31.3 


-916 

-126 

-17.7 


It is seen in Table 4*4 that the gain with a nonseparable distribu- 
tion or a separable distribution with an anisotropy of q = 2 would 
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FIGURE 4.4. RAY TRACING AT 80 KHZ TO SIMULATE OBSERVATION OF LF 
SIGNAL ON OGO 1. Ray commences at 1000 km altitude, 
vertical Incidence, and 55 magnetic latitude. 


S: 

L = 3.6766 X 

= 36.67° e 

= -13.66° A 

== 0.883 

RO 

4.46 keV W 

= 77.5 eV 

B: 

3 .2894 

32.27 

-18.42 

0.916 

1.62 

13.5 


Point B corresponds to the highest normalized frequency 
along the path. For further details see Dunckel and 
Helliwell [1973]. 
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yield small attenuation at the satellite location. The reason is because 

the satellite is far from the equator (X 37^) and there are few resonant 

b 

particles available at the satellite location because of two reasons: 1) 

Many particles have mirrored before the satellite location; and 2) the 
refractive index decreases with latitude which causes increasing energy 
of resonance and hence fewer particles at resonance. 

It is interesting though to notice that if the satellite had crossed 
the ray path a few degrees lower in latitude at \ ^ 32^, to which 
corresponds the highest value of A( 0.92) along the path, the damping 
would have been important for the nonseparable distribution. Such a 
satellite path would have been most interesting to test the validity of 
our model. 

For all distributions the Landau interaction is negligible, and 
the wave is interacting mostly with very low energy electrons in some 
portions of the path. 

This LF wave, though of much higher frequency than the waves con- 
sidered in the last section, has some identical characteristics. Namely, 
the wave propagates in the ducted mode up to a high normalized frequency 
and then rotates inward and propagates in the unducted mode. Both 
because the normalized frequency is high on the unducted part and because 
on the unducted part they eventually propagate very close to the resonance 
cone, these waves interact with very low energy electrons (i.e., less than 
100 eV) . Because there is a complete lack of data in that energy range 
and because satellite measurements become more difficult at low energy 
(spacecraft potential effect, etc.), whistlers may very well prove to be 
a unique tool for measuring very low energy particles. An experiment 
similar to the experiment described in the last section may be envisioned 


SEL 73-043 


98 



for that purpose. 


E. MR WHISTLERS 

Magnetospherically reflected (MR) whistlers were observed for the 
first time by Smith and Angerami [1968] and studied rather extensively 
by Edgar [1972]. An example of an MR whistler spectrogram is given in 
Figure 4.5. The explanation of such whistlers in terms of accessibility 
(see above references and Figure 4.6) is the following: These waves are 

generated at low L values and their frequency* normalized to the 
equatorial gyrofrequency of the field line of generation, is small and 
stays below the LHR during an important portion of their path. In the 
absence of ducts, the wave propagates with a large wave normal angle 
which generally increases with distance. Since the wave frequency is below 
the LHR, the refractive index surface is closed (see Chapter 2), there is 
no limit for the wave propagation angle and reflections can occur. 

According to the input latitude, the satellite can record different 
components which have been reflected a different number of times, labelled 
as MR 0+, 1-, 1+, etc., as sketched in Figure 4.6. 

Edgar [1972] explained the characteristics of the spectrogram of 
Figure 4.5 (and similar other spectrograms) through arguments invoking 
accessibility. We see several important features of MR whistlers in this 
spectrogram, namely, both upper and lower frequency cutoffs, and emission- 
like structures at the upper cutoffs. This last feature suggests that 
wave-particle interactions may play a role in MR whistlers and may account 
for both upper and lower frequency cutoffs. 

Thorne [1968] postulated that all but the Landau interaction can be 
neglected both because A « 1 and 0 is large along the major portion 
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FIGURE 4.5. FREQUENCY-TIME SPECTROGRAM OF A TYPICAL MR 

WHISTLER OBSERVED NEAR THE MAGNETIC EQUATOR. 
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FIGURE 4.6. SKETCHES OF RAY TRACINGS TO EXPLAIN THE FIRST 
THREE COMPONENTS OF WHISTLER SHOWN IN FIGURE 
4.5. The ray tracing parameters are given in 
Table 4.5. 
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of the path, and computed total amplitude gain over a few typical MR ray 
paths. He first used isotropic distribution <x — . He concluded that 
such distributions give too much damping and he resorted to distribution 
with a secondary peak around 10 keV. The peak causes Landau growth along 
part of the path, therefore reducing the damping. 

Edgar explained the upper frequency cutoff as well as "emissions" 
structures at the high frequencies purely on accessibility grounds. The 
high frequency distortions are caused by the presence of ducts, and density 
irregularities can account for the upper cutoffs. Edgar [1972] could not 
explain the low frequency cutoff but we can suggest a very simple explana- 
tion [R. A. Helliwell, private communication]; namely, the earth-ionosphere 
waveguide acts as a high frequency filter with cutoff frequency at ~ 1.5 
kHz during nighttime [Helliwell, 1965]. (Note that the spectrogram of 
Figure 4.5 was recorded in the early morning.) 

Therefore it seems of interest to re-examine Thorne's [1968] 
computations . 

We integrated the wave amplitude gains for separable distributions 
along ray paths whose parameters are defined in Edgar [1972] (see Table 
4^). We had to resort to separable distributions because of large values 
of 9 and we have developed a program for arbitrary 0 only for separable 
distributions. The choice though is not unreasonable. For cyclotron 
harmonics, the energy of resonance is very high and we have seen (Chapter 
3) that for high energy, the assumption of separability is reasonable. 

The Landau interaction is much less sensitive to nonseparability than 
the cyclotron interaction. 

We present in Figure 4.7 the results of our computations for the first 
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TABLE 4.5. RAY TRACING PARAMETERS TO EXPLAIN THE MR WHISTLER 
OF FIGURE 4.5. The parameters were chosen after 
Edgar [1972]. The rays begin at 500 km altitude, 
vertical Incidence, and latitudes between 20® and 40 
to intercept the satellite located at L ~ 2.4, 

X ~ 52 ‘^S. The diffusive equilibrium model is used 
with 50% IT^, 50% O'*" at 1000 km, T = 1600 ®K, and 
n(1000 km) = 1.15X10^, with an irregularity of the 


form n(r,L) = * 


expression and n = 1 


”de 

0.9 exp 


is the DE model 

■ 

The 


(L-4)- 

1.6 


sign of 9 
wise from 


is defined as positive measured clock- 


«0 


f ( kHz ) 

MR 

A 

g(o) 

WRo(keV) 

W^jCkaV) 

1.5 

1 - 

0.023 

-70.1 

0.314 

556 


1 + 


-102.3 

0.448 

821 


3- 


-84.1 

1.01 

1710 

3 

1 - 

0.046 

-69.5 

0.519 

217 


1 + 


- 101.9 

0.827 

328 


3- 


-83.5 

1.60 

750 

5 

1- 

0.076 

-68.7 

0.752 

119 


1+ 


-102.3 

1.36 

209 


3- 


t> 

CO 

00 

1 

1.48 

233 

7 

1 - 

0.105 

-67.5 

1.11 

1 — 1 

0 

00 


1 + 


-101.3 

1.63 

122 

10 

1- 

0.150 

-66 .6 

1.40 

45 .0 


three MR components . In this case, the parameters influence in a complex 
way the values of the gain and no simple discussion of the influence of 
each parameter individually is possible. We see already from Table 4.5 
that the ray path is fairly complicated. The ray path though is similar 
for different frequencies but the same MR component. The wave propagation 
angle is moving towards 90° as the number of the component is increasing. 
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When 6 is close to the Gendrin angle 6^ and O) is fairly close to the 

LHR, a few degrees difference in 9 may yield a big difference in resonance 

energy as illustrated by the difference in resonance energies for the 1+ and 

3- components of the 3 kHz frequency. It is important to compare the value 

of 0 at the satellite (which is close to the equator where the interaction 

is most important) with 9^. When 9 increases (as it does for increasing 

component number), two opposite situations develop depending upon whether 

9 is below or above 9«« Referring to Figure 2.4, when cu is slightly 

G 

above the LHR, N ^ decreases significantly when 0 increases while below 
0 , and increases very rapidly when 9 increases while above 9^. Corre- 
spend ingly, as the minimum energy of Landau resonance is inversely propor- 
tional to the strength of the interaction decreases (increases) when 

9 is below (above) 9 . This may be in favor of Landau damping for higher 
frequency and higher component number, assuming the ray path behavior is 
the same for different frequencies and the same component number- For 
instance, denoting by 9 the value of 8 satellite location, 

9^9 for f = 5 kHz and the 3- component, and therefore 9 would be 
S G ^ 

beyond 0^ for f = 7 kHz and the 3- component. 

The Landau interaction is almost always dominating, though neglecting 
the cyclotron interaction, as Thorne [1968] did, is Invalid in some cases. 

For example, for the case v = 4, q = 2 and the first MR component, cyclo- 
tron growth overtakes Landau damping for all frequencies. Apart from this 
case, the Landau damping dominates, but it would still be incorrect to 
neglect the cyclotron interaction, at least for v = 4. As pointed out in 
section 3E, though the cyclotron resonance energies are much higher than the 
Landau resonance energy for A « 1, the parallel component of the electric 
field, which excites the Landau resonance, is small. For increasing v, the 
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Landau interaction becomes more important relative to the cyclotron inter- 


action. For V = 5, the ratio of cyclotron gain to Landau gain is of the 

order of 10 to 207o and becomes negligible for v = 6 the order of 1 to 2 %) . 

The question arises of what is the importance of the gain around the 

turnaround point. In this case the approximation Eq. (2.71) is no longer 

valid because the ions were neglected, and we have to take the correct 

limit of Eq. (2.69). When 8 90^ and cd < cd , coscc 1, cos0 0 

LHR g 

and there is an apparent divergence of the order l/cos0 . But p ^ oo, 

Rm 

with the order (cos0)^ . At the same time d = |(A-m)tan0| w like 

m ' ' 

l/cos0 and the integral in Eq. (2.68) 0 like ^ ^ cos0 (see Appendix 

m 

A). It can be shown that H and F defined in Eqs . (2.68) and (2.69) 
remain finite. Therefore 


k. ^ (cos0)^ % =4 90^ 

ig 


(4.7) 


For V > 4, there is no interaction around the turnaround point. This 

result just expresses that when 0 90^, the resonant energy tends 

toward infinity, and as the number of particles must remain finite, the 

number of particles at infinite energy is null. (Boundedness of the 

number of particles means v > 3 and therefore from Eq. (4.7), k, 0 

ig 

as 0 5^ 90 . ) 

The damping is larger for v = 4 than v = 5 because of choice of 
normalization. The minimum energy of interaction is always bigger than 
100 eV and we keep the number of particles above 100 eV to be the same 
for every distribution. As already mentioned in Chapter 3, this means 
more resonant particles and hence stronger interaction for harder spectra. 

The influence of the anisotropy factor is more complex because of 
opposing effects (see Section 3C). Higher anisotropy means stronger 
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interaction because of increasing slope of the distribution with respect 
to pitch angle. It means also weaker interaction because of a decreasing 
number of particles at resonance (the average resonance energy increases 
and particles mirror closer to the equator) . Finally, the cyclotron 
interaction acts differently for isotropic and anisotropic distributions . 

The most important fact to notice is that for an anisotropy q = 2, 
the damping is never important enough to cause a sharp cutoff, in contrast 
to an isotropic distribution which could explain the higher frequency 
cutoff- Therefore, a distribution with a secondary peak in energy is not 
at all necessary to explain the observation of high frequency portions 
of MR whistler components as suggested by Thorne [1968]. The anisotropy 
of the distribution is sufficient for this purpose. Anisotropy must be 
present because of the loss cone (the more important, the lower the L 
shell) and isotropic distribution may not be likely in general, according 
to Sections 4 b and 4D. 

In fact, for v = S, even an isotropic distribution would yield 
reasonable attenuation. Moreover, let us be more explicit about the 
normalization we chose and express it in terms of differential flux at 

90^ pitch angle. At the satellite the total particle density we used 

-3 “3 

for ray tracing is 2500 cm . With our normalization (5 = 4x10 ), 

F 1 

there corresponds a differential flux at 10 keV (see Figure 3.5): 

4 -2 -1 -1 -1 

d$(10 keV) ~ 4x10 cm -sec -sr *eV , v=4, q=2; 

3 -2 -1 -1 

d$(10 keV) ^ 6x10 cm *sec .sr *eV , v=5, q-2. (4.8) 

Therefore, even for a value of \; = 4 , if we take a more reasonable 

4 “2 -^1 “1 -1 

value of the flux of the order of 10 cm -sec •sr *eV like Thorne 
[1968], we do not obtain much attenuation. 
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We deduce that the analysis of MR whistlers for monitoring the hot 
electrons in the magnetosphere is, at the present time, limited. Know- 
ledge of the particle flux is a prerequisite to drawing any conclusions. 

It may be interesting though in the future to test whether the high MR 
whistler frequencies may suffer high Landau damping relative to lower 
frequencies because high frequencies propagate sooner beyond the Gendrin 
angle. An experiment to test this effect should be done on satellites 
above ionosphere irregularities to eliminate the accessibility problem. 

An interesting experiment has been proposed by Ekigar [1972]. He showed 
that it would be possible for a wave excited close to the Gendrin angle 
to propagate in the MR mode along field lines, reflect and come back to 
the satellite (the "boomerang mode"). On the return path to the satellite, 
the wave would propagate beyond the Gendrin angle, favoring heavy Landau 
damping . 

F. CONCLUSION 

We have demonstrated that separable distributions proposed by Liemohm 
[1967] are incompatible with Angerarai's [1970] observations provided: 

1. The theory presented in Chapter 2 is adequate 

2. Particle fluxes at the time of Ar^erami's [1970] observations 
are comparable to fluxes observed by Schield and Frank [1970]. 

We have shown that the model of nonseparable distribution we have 
proposed in chapter 3 is consistent with both these observations and 
particle observations. We summarize in this section the implications 
of this new model, and the application of our analysis to a diagnostic 
tool for monitoring the energetic electrons in the magnetosphere. 

Our model of distribution involves a certain equilibrium between 
low energetic particles and electromagnetic waves in the magnetosphere 
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(neglecting sources and sinks of particles)- The fact that the hot 
particle distribution may be quite stable to electromagnetic waves casts 
a doubt upon the efficiency of plasma injection experiments. 

With our analysis, we have obtained valuable information concerning 
the hot plasma in the magnetosphere. Knowing the energy parameter v 
from Schield and Frank’s [1970] data, we assessed a lower bound on the 
anisotropy factor q of our model from the value of the upper cutoff 
frequency of the whistler shown in Figure 4.1a. From the condition that 
the difference in gain should be small between this upper frequency and 
any lower frequency, we assessed a lower bound on the coupling factor 
V of our model. This method could be used in a systematic way in a 

q 

controlled experiment which would involve a ground based transmitter and 
a satellite orbiting near the equatorial plane. Data similar to Angerami’s 
[1970] or Dunckel and Helliwell [1973] should be obtained with such an 
experiment. An experiment using MR whistlers and a satellite seems also 
desirable to obtain information about hot electrons. The boomerang mode 
looks most suitable for that purpose. 
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V. GENERAL CONCLUSION AND SUGGESTIONS FOR FURTHER WORK 


The purpose of the present work was the consideration of two prob- 
lems. The first problem involves the study of wave amplitude variations 
for arbitrary wave frequency below the electron gyrof requency and arbitrary 
direction of propagation. The second problem involves the determination 
of a realistic model of energetic electrons of the magnetosphere. 

The instability characteristics of a few whistler observations, 
whose accessibility characteristics were already thoroughly established, 
have been investigated, and limits of distribution function models have 
been set up. 

The theory we used to compute the wave amplitude variations caused 
by wave-particle interaction (i.e., whistler -electron interaction) was 
based upon several assumptions. In Section A, we briefly discuss where 
some of these assumptions break down. Section B will present a summary 
of our study. Finally, in Section C we give some suggestions for future 
work . 

A. LIMITS OF THE THEORY 

The theory we developed was based on the assumptions that; 

A) The expression of the imaginary part k^ of the wave vector k 
is obtained from the coupled set of Maxwell equations and Vlasov equation 
after linearization . 

l^il « |k| • (5.1) 

Our purpose here is to discuss succinctly when these assumptions may 
break down. From Chapter 2, Assumption A is no longer valid if particle 
trapping by the wave becomes important. This occurs if wave signal 
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duration is of the order of trapping time, or more* We look briefly at 
this question in Appendix B, where we take into consideration: 

1. Natural whistlers* In this case, though the duration of the 
signal is shorter than a trapping time, the question arises whether the 
electron can be trapped because of wave dispersion and space inhomogeneity 

2. Ducted whistlers. For unducted whistlers, the wave moves across 
field lines and interacts during a shorter time with electrons which are 
compelled to drift along magnetic field lines than for ducted whistlers. 

3. The region around the equator, because the important part of 
the interaction takes place there (Chapter 3). 

We used the criterion developed by Sudan and Ott [1971] and Dysthe 
[1971] and conclude that in those conditions, trapping can only take place 
for waves much more intense than average intensity waves [Burtis, 1969]. 

Assumption B insures that; 

1* The real part of wave refractive index need not include hot 
plasma corrections - 

2. The ray propagates along the ray path predicted by cold plasma 
theory and the group velocity keeps its full meaning. A ’’signal 
velocity” concept [Brillouin, 1960] need not be introduced. 

In Appendix C, we look at question (1) for parallel propagating waves 
We conclude that nose whistler dispersion cannot be used to test (1), and 
that hot plasma corrections should be included at the resonance cone. 

B. SUMMARY OF RESULTS 

Based upon the above mentioned theory, a general expression has been 
derived for the wave-amplitude variations of an electromagnetic wave in a 
inagnetoplasma * This expression *is algebraically simpler than previously 
derived expressions [Kennel, 1966; Brinca, 1972] and, as such enables 
easier computations and allows some conclusions by inspection; 
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1- The instability character of each harmonic of the interaction 


is straightforward • 

2* The validity of the low-frequency approximation of Kennel [1966] 

is tested* It is valid for 0 < 9 but not for 9 -♦ 0 ; Kennel [1966] 

G R 

approximation is incorrect for the Landau interaction because it neglects 

an important term, even for 9 < 9 * However this does not invalidate 

G 

his conclusion that at 0 / 0, there is growth for a cone of angle 
9 < 9 , though his value of 0 is incorrect* 

s s 

3* Significantly more anisotropic distributions are required for 
growth when 9 -♦ 0„ than when Q =♦ 0. 

A new model of distribution function has been proposed. This new 
class of distribution is nonseparable in particle energy and pitch angle, 
or in other terms, particles of different energy have different pitch- 
angle distributions * 

Our model is expressed mathematically by 

V 2 

exp[q(v /v) ^sin a] 

f(v,a) “ (5.2) 

V 

V 

where v^ is a normalization parameter* The parameters v> q> and 
have a simple physical meaning: v and q express the energy and 
pitch-angle dependence respectively, and v expresses the coupling 
between energy and pitch angle* 

This function can be almost marginally stable over a large frequency 
band for parallel and almost parallel wave propagation, in contrast to 
previous separable distributions such as Liemohn's [1967]. This model is 
more realistic than separable models. It predicts an increasing slope in 
spectrum with decreasing energy, in agreement with Schield and Frank’s 
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tl970] data. It predicts near isotropy at high energy, in agreement with 
Bogott and Mozer's [1971] data. There is indirect confirmation of our 
model, in the high energy range, by Kennel and Petschek’s [1966] mechanism 
extended to arbitrary angle of propagation by Lyons et al [1972]. 

Fortran programs have been developed to compute total wave amplitude 
variations caused by whistler-electron interaction, integrated over a 
ray path. The following conclusions were reached: 

1. For a whistler trapped in a duct, the parallel propagation 
approximation is valid for all frequencies (that is, below the frequency 
of untrapping A ^ 0.5) and is valid for nonseparable distributions for 
all frequencies except frequencies close to A ^ 0.5. 

2. From almost simultaneous measurements of electrons [Schield and 

Frank, 1970] and whistler observations [Angerami, 1970] we conclude that 
separable distributions previously proposed by Liemohn [1967] are in- 
compatible with observation- For the class of nonseparable distributions 
we have proposed, bounds on the values of the three parameters, v> q> 
and assessed, v 5 is determined by the experimental spectrum 

[Schield and Frank, 1970]. q ^ 5 is determined by the highest observed 
frequency, and 0*5 is determined by the condition that the 

difference in gain between two frequencies is small. The sensitivity 

of the parameter q is important. There is a strong suggestion that 

the flattening of the distribution around 750 eV in Schield and Frank 

[1970] is caused by Landau interaction, because a flattening in the 

distribution reduces the Landau interaction which eventually dominates 

for moderately anisotropic distributions when 0 9 . There is a 

R 

suggestion also that a peak in energy should exist around a few tens of 
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eV, otherwise the gain rates would be too large* 


3* It is possible to observe whistler waves in the magnetosphere 
at a frequency close to the local gyrofrequency , provided the observer 
location in latitude is far from the equator, simply because there are 
few resonant particles at high latitude* If the distribution is aniso- 
tropic, many particles mirror before reaching high latitude* Furthermore 
inside the plasmasphere , the refractive index increases with latitude, 
leading to an increase in the energy of resonance, and fewer particles 
resonate with the wave* 

4* Computations of total wave amplitude gain have been done over 
some MR whistler ray paths, similar to Thorne^s [1968] computations* 
However in Thorne’s paper the cyclotron interactions were neglected on 
the grounds that the cyclotron resonance energies are much higher than 
the Landau resonance energy, for A « 1* This assumption is incorrect 
if the spectrum is hard enough, because, though there are many more Landau 
particles than cyclotron particles, the Landau interaction is weighted 
by a factor which decreases with frequency- However the Landau inter- 
action always dominates after the first reflection. Anisotropic distri- 
butions greatly reduce the damping so that it is possible to explain MR 
whistler attenuation only by anisotropy without resorting to peaked 
energy distributions as Thorne [1968] did. 

Our model has some implications: 

1. This model involves a certain equilibrium between plasma and 
the electromagnetic waves that propagate in the magnetosphere. 
This situation was predicted by Kennel and Petschek [1966], and 
Lyons et al [1972] for high energy particle fluxes. Our analysis 
suggests that waves efficiently control the low energy particles 
also . 

2. Plasma injection experiments rely on the fact that increasing the 
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cold plasma density lowers the particle resonance energy and 
therefore increases the number of resonant particles. Wave- 
particle interaction effects should then become more important, 
and noticeable effects may ensue. However, at times when the 
distribution is quasi-stable (to electromagnetic waves) the 
efficiency of the experiment is doubtful. 

Extension of our analysis to a diagnostic method for monitoring the 
distribution of energetic electrons in the magnetosphere may be envisioned. 
However a selection in the type of whistler data is necessary: 

1. Ground recorded whistlers are of some value for the purpose of 
hot plasma diagnostics. They may be useful to confirm a known change 
in distribution or to show that the distribution is not isotropic. How- 
ever they cannot be used to obtain the precise shape of the distribution. 

2* MR whistlers recorded on satellites and generated close to the 
earth are of limited value. The particle fluxes necessary to yield 
noticeable effects seem high, and a precise knowledge of the actual flux 
is prerequisite to any further conclusions. High MR whistler frequencies 
propagate sooner beyond the Gendrin angle than low frequencies as the 
number of reflections increases. As a result, the Landau damping increases 
with frequency and it will be interesting to test whether the high frequency 
components suffer heavy Landau absorption after the first or second reflec- 
tions . To eliminate the problem of accessibility a satellite-to-sa tellite 
experiment is most desirable. The "boomerang mode" [Edgar, 1972] has the 
merit of requiring only one satellite, transmission and reception being done 
on the same satellite. It has the other advantage of using propagation at 
the Gendrin angle and beyond, therefore increasing Landau attenuation, 
and also, optimizing the transmitting power [Wang, 1970]. A strong 
recommendation is made for a detailed study of this experiment. 

3. A controlled experiment using a ground-based transmitter and a 
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satellite would be most valuable for the purpose of hot plasma diagnostics, 
provided that the frequency of the transmitted wave is comparable to, or 
higher than, the minimum electron gyrofrequency of the field line of 
propagation and provided that the satellite orbits close to the equatorial 
plane* In these conditions the energies of interaction can be very low, 
for two reasons. The wave propagates a certain time in the ducted mode 
with a high value of A (low energy of cyclotron resonance) after which 
it propagates in the unducted mode and approachs the resonance cone angle 
(where both Landau and cyclotron interactions become important) . Because 
the lower the particle energy, the more difficult the satellite particle 
measurements, this method may very well prove to be unique in measuring 
very low energy particles (of the order of a few tens of eV). 

4. Because of the scarcity of available relevant data, we could not 

make a systematic test of how best to conduct the process of measurement 

for the purpose of diagnostics. Two orders of sophistication in the 

process are foreseen. The first method would eliminate an accessibility 

study and would require only the measurements of upper and lower cutoff 

frequencies . From the measurement of the higher frequency a lower bound 

on the parameter q is determined and from the condition that the 

difference in gain for different observed frequencies should be minimized, 

a lower bound on the parameter is estimated (along the way described 

in Section 4C). For ducted propagation, the knowledge of the particle 

flux seems a prerequisite (determination of parameter v) • On the contrary, 

for unducted propagation, knowledge of v may not be necessary, as the 

conditions from one ray path to another may be quite different and from 

the observed differences, v itself might be determined (or, it might 

be deduced that a different energy law variation is needed to represent 
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the energy spectrum)* The second method would require the meas tirement of 
dG/dO), This method is more sophisticated than the previous one because 
it would first require a study of dc/do) due to accessibility (a problem 
only recently approached [Dantas, 1972]), to clearly assess the relative 
importance of accessibility and instability effects. 

C, SUGGESTIONS FOR FUTURE WORK 

Besides the experiments we have mentioned, several suggestions for 
future work can be made: 

1* Similar computations should be done outside the plasmapause. 
Going from inside to outside the plasmapause, the cold plasma 
density drops abruptly* There is a change also in hot plasma 
distribution, but apparently it is not the same change and 
probably the instability characteristics are different inside 
and outside the plasmapause. 

2. A program should be developed to include a general distribution 
function for arbitrary direction of wave propagation. In 
particular, this will allow a study of very low energy electrons, 
when a wave becomes unducted and eventually propagates very 
close to the resonance cone angle. At the same time, hot plasma 
corrections to the real part of the refractive index should be 
looked for. 

3. The stochastic analysis of Kennel and Petschek [1966], and Lyons 
et al [1972], which takes into account possible sources and 
sinks of particles, should be extended to higher frequency waves 
and low energy particles. 
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APPENDIX A. PROGRAM TO COMPUTE NONPARALLEL GAINS 


Separable distributions ; 

Inserting given by Eq. (3.6) into Eq. (2.71) yields: 
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(A.l) 


^iL ’^ic the contribution of the Landau and the 

cyclotron interactions respectively to the overall whistler mode wave 
interaction. 

In the last expression we have to compute integrals of the form: 
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(A. 2) 
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Inserting the coefficients B ’ given by Eqs . (A. 2) into (A.l) 

v.q 


yields ; 
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(A. 3) 


From Eq. {A. 2), there is an identity between distributions of different 
order in v and q: 


B „ = B - B „ 

V»q+2 v>q v+2,q 


(A .4) 


Similarly, there exists a couple of recurrence relationships between 
sucessive resonance harmonics. Using the identity 


t 2 . . 2 2m d 12] 


B ’ (d ) = B + [qB o“(q+v)B^,_^^ „] (A. 5) 

v,q m v,q v+2,q-2 v+2,q 


and using the identity 
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(A. 6) 


. . 1,2 2 . 2m^ J2. 

Vl = ■ 2 ■" ^2 \ ’ 

v,q m 2 v,q V,q ^2 v+2,q-2 

m 

These last two relations permit us to deduce all the coefficients 

j mn 4 ^ 00,11, from the values of and ♦ 

V»o[ V>Q V#q 

When V is an integer, the coefficients can be integrated 

analytically by means of relations involving different kinds of Bessel 
functions derived in Watson [1958], for example: 


CD 

/ 


a + d 


(A. 7) 


where and are the modified Bessel functions of first and second 

kind respectively - 

By derivation with respect to d, we can obtain the values of the 


coefficients B™ , m 0.1, and v even: 

mm d 5 ^ mm., 

®40 “ - 5 3a '"so' - »*<=• 


(A. 8) 


From the asymptotic expressions of I (x) and K (x) for 

n n 

|x| » n, 1 (see Watson fl958] for example), we deduce: 
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(A-9) 


Keeping the first term in Eq. (A. 9), we see that 


„mm 

^40 
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(A. 10) 


mm 


For higher values of v, the coefficient B ^ tends also 

vO 

asymptotically towards a l/d variation, and from Eqs . (A. 4), (A, 5), and 
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(A. 6) all the coefficients have this asymptotic behavior. 

For odd integer values of v? can also obtain the coefficients 

in closed form [T. F. Bell, private communication], though not by 

vq 

straightforward derivation. Even for even integer values of \j but for 

high values of v ^^d q, a closed form for B requires calculations 

VQ 

of high order derivatives- This may be tedious and it is easier to integrate 

numerically the coefficients B for 0 < d < d , and use the asymp- 

VQ ® 

totic form Eq . (A. 10) for d > d . 

a 

A listing of the program is given in Table A.l , and an example of 
the computer output is given in Table A. 2 - 
* Nonseparable distributions ; 

A program is developed for small values of 9, for a distribution 
function given by Eq . (3-22). From Eq- (3.25), is given by; 
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(v+nv +2n) B „ . „ 

q V-2+nVq,4+2ii 


V-2+nVq» 2+2n 


(A. 11) 


A listing of the program is given in Table A. 3 , and an example of 
the computer output is given in Table A. 4. 
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table A.l. 


LISTING OF THE PROGRAM TO COMPUTE THE GAIN INTEGRATED ALONG 
A RAY PATH FOR SEPARABLE DISTRIBUTIONS. 


1 

2 


SWATp [ V 

C THIS PROGRAM CDMPDTFS THC HQOF V)AVE AMPltTunP TAINS HM 

C JNTrGPATED OVER A PAY PATH WHDSE » AR AhET PS An F f,IVFM "V TMf' STANFORn 
C VLF GHnUP PAY TRACING PR'IGRAM, POSITIVE INnCATlVTl GAIN MFANS WAVE 
C GPOMTHIOAMPINCI . THP ENtHGETlC eLfCTRaN niSTlPUTinN IS PKOPfluT I riNAl TO 
C S INI alpha >«ttO/P + *NU» where ALPHA AND P ARE srSPFCTIVFLY THC ELfTTPON 
C PITCH-ANCLE AND MOMENTUM, NU GAN TAKE JNTFCER VALUES FPOH 4 TO 10* 

C AND a CAN TAKE EVFN INTEGER VALUES FROM 0 TO I D-?*1NT tUU/2 > . M IS THE 
C ORDER OF CVCLOTHCN HARHONICS. THIS PROGRAM INCLUDES VALUES FROM M=0 TO 
C A8SIMJ = 3 IM=0 FUR LANDAU HARHflNlC, M=l FOR FUNDAMENTAL CYCLOTRON 
C HARMONIC! 

C 


C INPUT: 

C 

C 

C 

C 

C 

C fPUTPJT 

C 

C 

C 

C 

C 

C 

c 
c 
c 
c 
c 
c 
c 
c 

C 

c 
c 
c 
c 
c 
c 
r. 
c 


FREQUENCY OF TUG HAvF(KH7.I, INITIAL ANO FINAL GROUP T(niES 
I L INFS ^6 ANU 47 J . 

L VALUE, LATI TunFIDFG.) , NORHALUED FRFOUrNCY, WAVE PROPAGATION 
ANGLEIDEC-I, dens ITV(CM-3I tGROUP TIMEISFC) I L I NE S 57 AND Sfl ) . 

FREQUENCY OF THE WAVFIKHZI, INITIAL AND FINAL GROUP TIMFS 
L VALUE, L AT! TUDF IIIFG. » . NORMALIZFO FREQUENCY, WAVE PROPAGATION 
angle (DEC. 1, nFNSITY(CM-3l , PLASMA FRFOUENCY NORMAL TZFO TO 
GYRUFREOUE^lCY t PARALLEL CDHPONFNT OF RFFRACTIVE iNOO^t MINIMUM 
RESONANCE ENEoCy FQP LANOAll AND FUSOAmCNTAI. CYCLDTRUM INTER- 
ftCIIONS IKFVI. 

GAINS iNTEGRATEn FROM INITIAL TIKE, INCLUDING R ES PECT I V CL Y: 

ALL HAP MONICS I , LANDAU HAR M(iN I C » F UNDAMF NT A L HARMUN I C, HIGHER 
harmonics, PARALLEL PROPAGATION APPROXIMATION, NU, AND 0. 

high PASMA FRFOUENCY APPPOXIMATION FDR REFRACTIVE INDEX, 8UT 
INCLUDE IONS, RELATIVITY NOT 1 NCUi nFfl *T 0 PRINT DIFFERENT VALUES 
OF NU AND 0 I BUT IN THF PREVIOUSLY DtFlNEP RANGF) , CHANCF 
accordingly lines 126 TO 14A, 191 TO 211, AND 254 TO 274, 

TO INCLUDE DESIRED NUMRER OF HAHKCMCS, CHOSE APPPOPPIATF VALUE 
OF M IN LINE 244. THE INTERVAL OF INTEGRATION IS nTATr>=3,02 SEC 
OR MORE AND CAN RF CHANGFO APPROPRIATELY IN LlNCS 63 AND 64, 

THE DISTRIBUTION FUNCTION IS NORMAL I ZFU TO OF <=0.00 4, WHERE 
DE=PAT10 HOT plasma DENSITY FOR ENFPGY>100FV TO COLO PLASMA 

density; chose appropriate value df of in line 246. 


FUNCTION factorial 
FUNCTION FFFINJ 
IF IN.EQ.O) GO in 2 


3 

4 

5 

6 
7 
6 
9 


FFF=1 

DO 1 I=i,N 

1 fff=fff*i 
GO TO 3 

2 FFF=1 

3 RETURN 
END 


10 
1 I 
12 

13 

14 

15 

16 
17 


FUhCTIDN AN0(NU,NQ> 

ANO= I 

J0M = NQ/2«-l 
DO II JQ=1,JCM 
II ANQ= (FLO ATI JO 1-0.5 >«ANg 
ANQ=AN0*2«(NI>-3) /FFr(NO/? » 
RETURN 
RNO 


16 

19 

20 
21 
22 

23 

24 

25 


FUr 1C T I n 4N 0 ( NU , NO ) 

RNO- I 

JO«=NQ/2+^l 
DO 10 J0=1,JCM 

10 BNO = GNOy ( NU/2.-2 , + Fl.OAT UQ) J 
ftNy=FFE ( N-JZ? )*BN0/ Z 
RETURN 
END 
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26 

27 

2d 

2y 


3 0 


3 2 

33 

34 
3f> 

36 

37 
?8 

39 

40 

41 

42 

4 3 
44 


46 

46 

47 

48 

49 
S3 

5 I 

52 

53 

54 

55 

56 

57 

SB 

50 

6 0 
6| 
62 

63 

64 

65 

66 
67 
6d 
6 9 
70 


DIMENSION S120) ,SC<20>,SU20),SC1 120 » 
niKENS I1‘>M C.US 120 I • Qvisr ( 20J «0MSL < 201 f GWSS 120 J 
ojfiFNSmN Of. VIS 1201 
DIKFNSI ON 00 (11 1 2 1 5 1 1 

C «\ 0 ’ 1 NTEGRAL rjJOM 0 TO Pl /2 OP S 1 NX*CDSJC**( NU -3 1 * 00 . UP^TANX )**2 . 
r 0 VARIES > ro 5 BY STEP O.U 

G 60 l«» 0 tll- 010=01 ♦ BOl*r*t?» = R 1 0 = . 1 ) ,FTC... 

C I IsNMVAP 1 AT IDN, VA»/I 4 linN» I 3 =D VARIATION 

OIMENSI^•'^ BNmt 6 , 2 ) 

C n FU9 NO VAtt*V NU-3 12 0 VAR, - Q/2 + 1 13= R V AP. I AT1 1 0:il 

C p\(llrb«ll=B( 4 T 0 t 2 r 2 >t liH ( 1 0 r 6 , 1 > = B t 5 » 0 1 2 t 2 ) » E TC . . * 
r AKii 1 I - S.PWm 4 . 2 lO, ? 1 . 3 N 1 10 « 5 1 2 I I 5 r 2 1 0 , 2 > t € TC, ♦ . 


INTEGER A:J,AC 
C AU=NUiAy*Q. 

PEAL rtiK ,KIKLAN*kI KCY l 

PEAL Lft,KlG 
pPftl. INVL 
REAL LAT 
LDGICAI Tt'ur 

nr ?00 L=l ,2 


nn 200 NU=ltU 
REAP <5,241 «0lNLr,1.,l) 
24 fORylM (FQ,6> 


or 200 K = 1, 10 ^ ^ 

HFAlXS, ?31 1 R01NU,L , J*l K-l »*5 » , J=2, 6) 

23 FnOMATI 5E16. 7J 
200 MnTINUE 

C PEAil 1NPUI DATA FRDH COL > “1 ASMA RAY TRACING 
10 00 1 PHAD <bi46Jl FKC. , TGNEy.TGLAST 
460 F()PMAT(P|0.4tFlOw6,Fn*6» 

l*RITn6,4S50l E«C , TONE Vi ,101 AST 
4550 FPR.^AH * IFKC. = ' , FID, 4, * TG D=* ,510-6, * 
C FKr. IS FRfOUrNCY SIGNAL IN K^7 
TPUC = -fALSE. 

nit 113 l = l» 2 C 
rows (1 1 = 0. 

GWSI I 1 = 0, 

GwSLM) = 0. 

GWSCin = D. 

GWSSUl = 0. 

118 CGNMNUE 

IDODO PEACH 5,490, r NO*90Ol 

% INVL ,lAT,t A, THETA, igPNS,TG 

490 rriRMATlFlO-5 ,510.2, FID.4,F10.2, I 10 
T»UC = .NDT.TPUC 
IF (TRUC) Tr,CLIi = TG 
J F 1 4 NOT . T i>U C > TGri E W=T G 
OTATG = ARSC TGNCw-TGnLn ) 

IF ( DTATG-LT ,0.02> TRUC = .NOT* TRlJC 

jr ( ptatg.lt ,0.021 on to ioodo 
PEGHA n = 3.V4l502/l{n 
IF (1 Hhr\.FQ -00. 1 th^ta- 39,99 
THFTAR = THP TA^nCr.PAO 

PFTA r LA*y. 904*SOPT(fLI]AT( IDENSl J /FKC 


TGI AST = » I riO.M 


nci = n-LA) *TAN<THFTAP J 
plan = LA’^TAMTHFTARI 

c mCkS IS hFNSITY in rapt ICLLS/CUf»lC CENTIMFTFF’S 
c ^FTAsPLAStn FJFOU‘=NCY/GYTa'=PEqutNf.y 

C LAcWAVE FkrQUFNrY/CYRnFRr-juENrY 

C TH^TA=ANGLF flF PP9RAGArrHN WITH PFSPEFT TO THE •’AGNETir FIELD 
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71 

72 

73 
7^1 

75 

76 

77 

78 

79 

80 
81 
82 
83 
8^ 

85 


86 

87 


Q6 

R9 

90 

91 

92 


93 

96 

95 

96 

97 
9{i 

99 

100 

101 

102 

103 


106 


105 

106 


137 

LOd 

109 

no 

m 


costa - AOS( COS( THFTAP.) J 
SINTA - ABS{ SINI THETAP ) 1 
OFLTA ni-l/( A3->LA)«‘‘2 
CnSTAD ■= CnSTA-LA*D£LTA 
TALf^HA^blMTA/ (Z*lCn$TAO M 
C THETA HFP.E is FXPPESSEO IN OTCRFfS 
r ALPHAS l/SOR 7(1* TALPHA**21 
P.IND = PFTA/SC»T UCnSTA-LAOOELTA J*LA) 

C P.INO = PEFftACT 1 VE INDEX 
P10=.02 

PPTNn=:=ftfcTA/SCRTt (1-LA»*LA1 
C PP. INO-PA»AL< FI PCFRACTive INOFX 
epAP.* R IND^COSTA 
upft>aRrNn«S INTA 
M = COSTAo 1 1 /LA-I.A I 
SO = f < t A« I l-LA'Coan 

SPa.<6FTA/LA 1**2 
r SP = STI X P PAKAHFTFF 

PP1*AHSU A-1 |/IPRlMn*LA) 

C SO* STIX n PARAKETtP 
rO=2*A8Sl SPl 

PBIWAT = U-LAJ/(Rf>AR*LA) 

C PklWAT IS THE lONP. MOMENTUM OF RCSONANCF 
C PR l=MnMENTUM FpR PARA(.LFL PfI OPAdA T 1 ON 
WO = 512. 

^;r 1 - WD*( PP 1WAT1**2/? 

PPOWAT 9 l/PPAR 

WLAN - W0*( PR0WAT>**2/2 

C PROWAT IS THF PFSC'NANCt MCfAfNTUM FOR lANCAU iNTFPACTinN 
PWCIHL = wn/WLAN 
C VI0“M0*C**2 

C WCl and WLAr; apf PFSPFCTIVFIV the minimum fnfrgifs of PFSCNEAMCC fop 
C FUWAHFNTAl AND LANDAU PHSONANCCS (IN KFVl 
WP|7E(6,680) 

680 FORMATt'OS* L *-AT LAMRO* THFTA BETA 

S OFNS Tf; MUCOS WLAN WC 1 M 

WRITE(6»6V3J INVLtLA I,L/' , THE T A , BFT A * I t-F N$ r TO , R PAR ,WL AN , WC I 
6 70 FflRMAT(FlO^ 5,F13,2tFl0.6,fl3.2iFlO.?.UO*Fl0.6,Fl0.2, 2C11.3> 
WRITfct6<660l 

46 0 format (» OUT GML GWC GWS PGW 

$ m fj* i 
K-0 

ir ('ucu.oi on to 63 

157 CONTINUE 

CMl =(( RFTA/l.fl 1*62*1 l+CDSTAl + BRER ♦♦?«L A/ M I /CO 
CM2 = I ( HFT.\/LA )♦*?♦( 1-C[ISTA 1 ♦ ftPE R* * Z*L A/ « J /C 0 
t C Ml cCOEFF ir, I FNT OF HESS ELF FUNCTION OF ORDER (M-U IN KIG ANO 
C CM2=CCEFF. of Sf<SEL FUNCTION OF ORDER (M*l) 

^3 pvt* ^ ACS ILA-H)/I RPAR*LA» 

C PRM = r^QPMALlZFL' MfJHFNTUM GF MFSONNANCF I TO P39H0*C) 
n -AFISI (LA-MJ *'SI NTA/CnSTA) 
ir lO.PQ.O) CD TO 159 
C STATEMENT lb9 FOR PAPALLFLE PROPAGATION 
C BEGIN OFFINITION OF AND Jl*»2 TERMS 

C INTERPOL AT I ON FC)P HIO-OI 
55 RJ-10*0 
J=INT(PJ > 
on 93 L=l*2 
DO 95 MU-1 » I I' 

IF ID. IT. 0. 1) RNINUil.LJ “BOlNU^LfH * 0+6261 OD* ( HOI NUf L f 7.1 


125 


SEL 73-043 



11 ^ 


111 

llA 


115 
i t6 
117 
I 10 
IIV 
120 
121 
122 
123 

12^ 

125 

126 
127 
12C 
!?•> 

130 

131 

132 

133 
13A 
1 35 
U6 
137 
1 38 

139 
l^tU 

141 

142 

143 

144 
1^5 
14i. 
a47 

140 
14>^^ 

150 

151 

153 


1 54 
155 
1 56 
157 
l5o 
15'* 
160 
161 


f -n0tNU,L,lJ J 

IF I (O.Gr.O.ll .4WD, jr)*1T.5.J 1 

$RrCNU»l,l) == (3J-J) (nO(NU,L , J»2) - R01NU,LiJ»m ♦ 6 D(NUtL,J+IJ 
C ASYrtPTDTIC FXPRFSSION WHEN D>5 

rF(D.GF.5-. » BN(M'J, I ,L)-P01 NU»Lt 51 1 
95 CONTI NUF 
C END INTERPOLATION 
C COMPUTE 6N0 tOQtUll f=OR 

no 12 L=l»2 
PTRUF=L-l 
DO 12 J0 = l,5 
JNM»11 -2«J0 
N(jTRUF=2*JQ 
on 10 jr;u = l*JNrt 

PN{ JNUr JO^l » Ll = RN(JNU,JO*L» - PN tJNU+2 . J <3 , L 1 

10 continue 

12 CONTINUE 

C END OF DEFINITION UP J0*P? AN" Jl««2 TFPHS 
IF r^.fvE.O) GO TO 163 
C BEGIN LANDAU TEE F 


no 

5 100 1 

I=l» 

10 

IF 

i 1 

.FQ, 

.1 ) 

AU-4 

IF 

{ I 

.FO. 

.2) 

AU=4 

IF 

4 I 

.EU. 

.31 

AU-4 

IF 

4 [ 

,F0. 

.4) 

AU=4 

IF 

4 4 

0, 

,5) 

AU-5 

IF 

1 I 

,EU. 

»6 1 

CU = 5 

IF 

I I 

.FP.7) 

AiJ = 5 

IF 

( I 

,Fe. 

-6> 

A4/-6 

IF 

( 1 

*( 0. 

.v> 

nU=6 

IF 

i [ 

.FO. 

.10 4 

AU = 6 

IF 

4 I 

.F(i. 

. 1 } 

AQ^O 

IF 

4 I 

.FO, 

.2 4 

A0 = 2 

IF 

1 [ 

• EH. 

.3 1 

A0 = 4 

IF 

4 1 

.EO. 

.4 » 

Ag-c 

IF 

4 1 

♦ FO . 

.5 1 

Ag-j 

If 

4 I 


.5 1 

Ag»o 

IF 

t I 

.EC , 

.6 I 

AQ=2 

IK 

4 I 

.EQ, 

.7 J 

AO-4 

ir 

( I 


,U ) 

A0 = 0 

IF 

4 I 

.FQ. 

.9 ) 

AD=? 

V 

4 I 

.EU. 

.101 

AQ = 4 


NLI^AU 

NC-AQ 

JNU=NU-3 

JQ=N0/2^1 

Tl - IRPAP^PPCR I **2«RN1JNU*2 rJQ* I J/C0**2 

IF ID.LT. 0.2 I T2= D* ?1NC1 1 NU» NO* 2 1 *RPAP* RPfR « S D* V/CO * *2 

\h t0.r,[.3.2» 

IT? = UND*2 1 J-;, 1 1-INU*N0*2 1 ♦FN UNU+2t JO+It 1 1 ) 

* p AR CP P FP * SO*- H /■ I n*r. O'* «'21 

T3 - ( Sn*H»#*2*5NUNU,J0 + l, 2 I /CO* *2 
SUNLAN =-I NII*Ng» "I Tl-T2«-T31»4 

SUML AN^SU^ILAN*! PIO/P^M **(NU-b I * I COSTA-LA I /It 
If n>,GE,5.» SU>LAN - SU>^LAN*=5/D 
SUW=SIJ'*LAN 
sn I = SU‘3 
SL ( I J i S UM. AN 
5100 CONTINUF 
C FNO LANDAU TE^N 
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16 ? 

163 

164 


16^. 
I isU 
167 

16V 
I/O 
IV I 
I7i 
173 
1 fU 

175 

176 

177 

17b 

170 

lao 

i;n 

U)3 
I »;4 

1 .\t. 
137 
1^33 
139 
l')0 


1’->1 

m 

194 
1S5 
196 
14 7 
196 
199 
?30 
201 
202 
29^ 
2 04 
295 
206 

20 7 

2on 

209 

210 

21 L 
212 
213 
?!<► 


^=l 

Cn TO 157 
163 COWTINUE 

r E«D COMPUTINtJ nif CfJP4FKl6NT 8( NU , 0 i M ,M M 01 
C BFOIN TERMS J2**2 AND 30*3 2 AMO SUPERIOR CPDERS 
C rjJNPUTE PtMU,0*2»21 

14 NUM-9 
mtruE=m+i 
PM^KTRUE-2 

IF ( NIH.LT.O 1 GP Ta looi 

nn 2 8 

ONI 12-NU, 6, 1 I = (fM(NOrl*l> - 2 * ( NU» 3 ) ( NU*2 , 1 1 2 1 /O* *2 

8Nll2-NUr6, ? HI 

28 continue 

MjM=NUM-2 

IF INUM.LT.O) GU rn 996 
on 29 NU>-l,NLW 

PM12-NU,5,11 - 2-f!Ur6, U - RN ( I 0 -NiJ , 6 , U 

RM 12-NUi:: f2 I - -HN(NII,2, 1 >+3N< 12-NU|5, 1 1 I /?+'*'! ( N(l*? , I t2 > *7 
4 *1 1/DI **2 

29 CONTiN'jr 

69 V CONTI NOE 

C COMPUTE : A:\I9 AND F\(NU»0- = C A\D 2,0 r2l 

0-^ 15 i = lt2 
DO 13 311 = 1 »3 
NDTRUF-2*30+2 

go rn von 

Di; 1 1 3NU = 1,JNK 

nMll?-jNlJ,5- JO,L 1 = AN{ l?-3NII ,6-JQ,L ) - PN I 1D-3NU , 6-3 0 , L I 

I I cuntinue 

13 CCMTINUF 
99B CONTIN'JE 

15 COi<;tini.ie 

C TOP TPkNS J?*=«2 AND J3*3 2 
C RFGIN ri^ST CYCLCTKDN 


no 

51 J1 

= 1, 

19 

IF 

( 1 

• Eg. 

1 1 

AIJ=4 


( I 

,*^0, 

?) 

A 11^4 

IF 

I 1 


3 1 

AU = 4 

IF 

1 I 

• f:j. 

41 

AU=4 

K 

1 I 

.(■Q* 

5 1 

AU=5 

1 F 

i 1 

.FO. 

6 1 

At/=5 

IF 

1 I 

.FO* 

7) 

AU=5 

IF 

{ I 

,F(0. 

8) 

AU=6 

IF 

1 I 

.FO. 

9) 

Atj=S 

JF 

I 1 

.to. 

131 

AU = f, 

IF 

( I 

.ro. 

1 1 

AD = 0 

IF 

1 I 

*t3. 

? 1 

AD = ? 

JF 

I ( 

.tin. 

3 1 

•*<J=4 

IF 

( 1 

.ED. 

4 1 

40=6 

IF 

1 I 

. = Q. 

5 1 

A0=0 

IF 

( r 

.rf}. 

6 J 

A9=2 

|F 

( I 

-FQ, 

7 ) 

A-g=4 

ir 

1 [ 

-Fg. 

8 ) 

AD=0 

IF 

r 1 

• CD. 

o J 

A-j=2 

IF 

1 1 

.FQ . 

lOJ 

AD = 4 

Nlr- 

.AO 





:A0 





NO/?+l 
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Zlb 

Zio 

217 

?AC 

219 

220 
221 
222 
223 

Z?A> 

226 

227 

221 ) 

229 


23 0 
231 
2J2 
233 
23^ 

235 

236 

237 

238 

239 


2^0 

241 

242 

243 

244 
24'i 

246 

247 

245 
2 49 
2Si0 


251 

2f>2 

253 

254 

255 
?S6 
26/ 
25 0 
259 
2 60 
261 


jM=Nl/-3 

12 = 24C»'Kf V2 * ( ^n| 12'JN, 7-JQ,2 - ('lO+NUI * 

% PN< 12-Jt- ,<>-J0,2 J > 

IF lO.LT.O.?) 12=0 
IF \ l^^S\ V.) , GT. i J T2-0 
IF 1,1,2 

1 rrmiNUE 

ir ( I .GT . I ) G'i n: ? 

rPi2=CMi 
FM1 = CM2 
CM2=Crtl2 

2 Tl - C'n*>*2 ♦ ( P'nIIJN, JO, 1 >*‘^*F'Q/LA - I Nt}+*’UM>»ri < J\. JO* I , I M 

IF TlwO 

IF IM.FJ.Z) 

t.Tl - r.^U<'*2 > I A’ JO, A - (Ng*NUl<'f»M| J*J, JO* 1,2 ) I 

13 = CF.2*'<'2 * 1 1?’JN, 7-JO, 1 J A - (N0*9(,/) * 

% PN( 12-J^‘,4'J0, 1 I ) 

?F (rt.FO.3) 

m - Cl I ■'“2 * ( 12-JU, 7-J3, 1) «W*N0/LA - (MO+MilJ * 

4 3N( 12-jr 1 \ I 

IF ( (l),l T .0,2 I »)^,F0-3 I J 11=0 

IF (D.LI,?.?) 13 = 0 
IF ( lABStM I .Cl.l 1 11=0 
WKCVC - T1+T2+13 

SM«0YC = 'iU*'rvC*< PlO/P'- '+J *»< NU-3l*Aft«i( inSTA-LAl/tM-LAl 1 
ir <Ll.r.F,5.l ?i;*1CYC = SU'1CYC#5/f) 

<;r<i) = r.LJv.cYr 

su i = su I ♦ Sf. n I 

IF ( M.FO. n sen r I =su^<cvc 

5 101 OGhjllNUE 
C FNO FIrfSl OYCLOlAnri 
C lETijPN TO XlI.liFR HARMONIC 

ir (.31 154, I «f,, 153 

158 W--M 

GO TO 15 7 
156 M = i 1 

If I IA5S ( 'U , GT.3 > GO T(' nOl 
GO TO 157 

C CND FflUPM rn HtGHFH HAPMOMC 

loot CONTINIJF 

159 CONTTNUF 

c subprogram fmp padali.fl gain 

nE = .0J4 

plat = LAT«Drf;yAn 
X = SIM PLAT I 
C INVI. = INVARI AMT LATITU'TI; 

C OF=PATIO hot PMSMA O'^MSITY fop EN6RGV>IOOEV TO COLO PLA$«A IICnSITY 
C LAT=l ATITUOE 

C fU-MOBNAU iFJ) “lACI'FTIC FJFLD TO FOUATORIAL VAUlF 
WT=S0Pn l*3ftX**2 in l-A*^*^l**3 
Plft-PiO/PPl 
00 5102 1-1,10 
IF ( i .fg.l 1 AO-4 
IF (T.F0.21 AU=4 
IF (I.E0.3> AU=4 
IF U.llO.AJ AU = 4 
IF II.eO.M AU=5 
IF ( I *F0.61 AU-5 
IF I1.F0.7J AU=5 
IF U.FQ.8I AU=:=6 
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w. 

?h3 

26 5 
Z6u 
26 7 
26 » 
26-? 
270 
771 

272 

273 
27^, 

275 

276 

277 
27 H 

27'V 

2ft0 

2BL 


2P? 

ZftiV 

?ft 5 

266 

2 H 7 

268 

7B9 

290 

291 
2'»2 
29 3 
29^ 

295 

296 

29 7 

298 

299 

300 

301 

302 

303 

304 

305 

306 

30 7 

30 B 
309 
3 J 0 

31 1 

312 

313 


314 
31 5 

316 

317 


IF (I. to . 9 I AU=r, 

IF U. to . 101 AU = 6 
[F tl.FQ.l I 3 Q =0 
IF ( 1 .FQ .2 1 A"j = 2 
IF U.fg .3 I = 

IF n .tO .4 J AQ = 6 
IF \ A 0=0 

IF n.FO .6 ) 30 = 0 
ir n.ea .6 ) AQ = 2 
IF ( I ,E 0,7 J A 0=4 
IF IT.tQ.fl J AQ -0 
IF I I .to , 9 > A 0 = 2 
IF U.rO.XOl A 0^4 
iVU = A'J 
^ 0 -Ag 

COEfKC='-AflC( KU,Mv;j*I 3 .l 4 l 592 /B)*De*(l/ 6 T»**<N 0 / 2 » 

PSU.HsBMOl W»N 0 J »(2 + < l-l/LA) + NO) 

% *Pl«'»*<NU-3 ) 

PKIK = -ODHFKrfPSU*^ 

AnVGC=SOPTUA)* S 09 .T( l-LA»** 3 * 2 /DF,TA 

PW 1 W*<-PKI K+APVlir*H iNi: 

C AV&C AnO APC&C = > 4 DDIJLUS DF GP’flOP VFLOCI T 1 ES (GENF'? At AND P AP Al tELCASFS I 
C NnftrtALlZED TO C 
C .rflw*s AftE JMl I 3 MEGAI /Rb( 0 , 1 «= 0 A) 

C PKIK=KIK FOP PAP ALL FL PPlJI'AGAUON 

PC 0 l:Fr» = -l 600 C 3 *BTOLA*PP IND/ I«^VL »*3 
PG 1 0 00 =P COE FG*P A I -C . 

C tNO SUBPROGRAM FOP PARALLEL GAIN 
CnCFKl- COFFKC 
4 *CALPHA/( CrSTM <<*2 
K IKL A?g = COFFK I*SL I I ) 

KIK-Crif FK 1 *S ID 
KIKf.V l=COFFK I»SC U I J 

AVGC = SOPT (LAafCOSTAD I*( 4 *(CnSTAD I *« 2 f $ INTA*« 2 M / 

I. ( PFTA«C 0 STAt 2 *LA*l 1 -m LTAJJ 
CUFWI W-AVGC^RINO 
VJIWLAN--KlKLANt'CnfWIiJ 
WCYl = -K IKCYl^C OFWIM 
WiW=-KlK«COFWlW 

CnEFG=-l 6000 0 «L A*«T*Q 1 ^J 0 / T NVL* + 3 
G 1000 =»GDEFG*KIK 

gioool-coffg*kiklan 

C 1000 C=C 0 EFG*.kIKCYI 
GIOOOS = GIOCO-OIOOOL-GIOOOC 

IF I AftS (PCIOOO I .GT. 1 . C -30 \ BATIOG = GIOOO/PCIQOD 

IF (AHS(Pt;i 300 l.Lr.l.F -30 > RATIGG ^ 123456 

IF ID.EQ.Ol RATIOG^I 

COEFGVJ = 86 R 6 *FKG *6 .233184 

PGW = CtiEFGW^nTATG^PWlW 

GW = CDFFGW*D 7 ATG*W 1 VJ 

GWL - rorFGW#OTATCf‘HlWLAN 

GWC - CD£FGW*DT 4 TGFWIWCYI 

PGws < n - pghsi 1 1 ♦ «'Gw 
GwS ni = GWS C I J * GW 
GWSLU 1 = GWSUr I ♦ GWL 

Gwsci n = Gwsrn i + r.wc 
cwssin = Gwsn j-gwslh i-Gwscn ) 

WRtTf:< 6»4301 GWSI n »GWSLl I \ tCWSCtl >»CWSS< 1 1 » PCWSM J , AU^ AQ 

430 F 0 PHAT( 5 E 12 . 3 » 214 I 
5102 CONTINUE 


IF ( TG.EQ.TGLAST ) GO TO lOOOl 
GO TO 10000 
909 STOP 
FNO 
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TABLE A. 2. AN EXAMPLE OF THE COMPUTER OUTPUT OF THE PROGRAM LISTED IN TABLE A.l 


FKC» U.7DO0 TGO=: O.QOODOO TGl.A ST= 2.011502 


L 

lat 



THETA 

BFTA OEMS 

TG 

Miicns 

WLAN 


.4G1 


4.61105 

24.11 0 

.50-11 


1.2 7 9. 

76 345 

0. 

0204 

1 9.51 

0.672F 

00 

0.625E 

□0 

GWT 

Owl 

owe 


CWS 

PGM 

MH 

0 






-O.f>06e 01 

-0.190F-01 

-0 .603 F 

01 

-3,fl90P-OZ 

-0.6C6F 01 

4 

0 






-0. 757F-01 

-O.802E-O2 

-o.7ort 

-01 

■3.376F-02 

-0.767F-01 

4 

2 






0.784E 00 

-Q.507E-52 

O.702F 

00 

0. 709F-02 

0.7846 00 

4 

4 






C.66 0E 0 0 

-0.304F-02 

0 .657F 

00 

D. 6236-02 

0.66 OF 00 

4 

6 






-0.32ef 01 

-0,124E-01 

-0..326F 

01 

-0.707F-03 

-0.3ZTE 01 

5 

0 






-0.349E-01 

-0.321F-J2 

-0.319E-01 

0. Z046-Q3 

-0.1326-01 

5 

2 






0.269E 00 

-0.115e-02 

0.290F 

00 

0.391F-03 

0.290r 00 

5 

4 






-0.J49E 01 

-0.730E-02 

- 3.149F 

01 

-0.6016-04 

-0. 149E 01 

6 

0 






-0»l37E-0l 

-0.145F-02 

-0.123F- 

-01 

0.135E-04 

-C, 126F-01 

6 

2 






•J.958E-01 

-0.37VE-03 

0.962F- 

-OJ 

0. 2456-04 

0.965E- 31 

6 

4 






L 

LAT 

L AMhOA 


theta I 

BETA 0FM5 

TG 


MlfCDS 

WLAN 


WCl 


4.AU25 

23.51 0 

.52 74 

2 

.94 10, i 

07 34? 

0. 

0439 

20.15 

0.610E 

DO 

0.506E 

00 

fJWT 

GVrL 

owe 


GW5 

PCM 

NH 

0 






-0.135E 02 

-0.375E-01 

-0.1356 

02 

-0.173E-01 

-0.135F 02 

4 

0 






-0.?59P 00 

-0.175E-01 

-0.349T 

30 

0.6956-02 

-0.36 IF 00 

4 

2 






0.16 01 

-0.102E-01 

0.167F 

01 

0. 137E-01 

0.167F 01 

4 

4 






0.147F 01 

-0.616E-02 

0-l46r 

01 

'>.12^6-01 

0.147E 01 

4 

6 






-0.776E 01 

-0.252E-01 

-0.773F 

01 

-0.1416-02 

-0.776P 01 

5 

0 






-0.173E 00 

-0.661E-02 

-0-l67b 

03 

0. 3806-03 

-C.1706 00 

5 

2 






0.h52£ 00 

-0.239F-02 

0»65^r 

00 

0.701F-O3 

0.654F 00 

5 

4 






-0.376E 01 

-0.151F-01 

-0.375E 

01 

-0.1276-03 

-0.376F 01 

6 

0 






-0.725e-0l 

-0.3036-02 

- 0-695 F- 

‘01 

0. 2716-04 

-0.703F-DI 

6 

2 






0.22VS 00 

-0.912F-03 

0.230C 

00 

0.511F-04 

0.231F 00 

6 

4 






L 

LAT 

1 ampoa 


THETA ftFTA 0£\S 

TG 


HUGOS 

WLAN 


MCI 


4.61 1 hS 

22.62 0, 

► 54 90 

2 

.52 10.44 339 

0. 

0733 

2 0,96 

O.503E 

00 

0.3936 

00 

n.jT 

GWL 

GwC 


CWS 

PGM 

MO 

0 






-0.236E 02 

-0.557E-01 

-0.236E 

02 

-0.2486-01 

-0.236F 02 

4 

0 






-0.105E 01 

-0.261E-01 

-0.103E 

01 

0.952E-02 

-0.105F 01 

4 

2 






0.272E 01 

-0.162E-DI 

0,27?'= 

01 

0.1966-01 

0.2 72 F 01 

4 

4 






0.254F 01 

-0.936£-Q2 

0.253F 

01 

o.ieoF-oi 

a.253r 01 

4 

6 






-0.1 A6e 02 

-0.3H5E-01 

-0.146F 

02 

-0.212E-02 

-0.K6F 02 

5 

0 






-0.550E 00 

-0.1036-01 

-0.540F 

00 

0.5466-03 

-0.545F 00 

5 

2 






0.114E 01 

-0.374E-02 

0.114F 

01 

0.1176-02 

0.115E 01 

5 

4 






-0.772E 01 

-0.2366-01 

-0.769F 

01 

-13.196F-03 

-0.772F 01 

6 

0 






-0.253E 00 

-0.4866-02 

-0.248F 

00 

0.397E-04 

-0.250E 00 

6 

2 






0,435e 00 

-0.132F-02 

0.4366 

00 

0.600E-04 

0.437F 00 

6 

4 






1 

LAT 

LAM80A 

ThFTA 8FTA OEMS 

TG 


HUCDS 

WLAN 


wci 


4.61193 

22.38 0. 

,5631 

2 

.24 10.66 336 

0. 

093 7 

21.47 

0.555F 

00 

0.334E 

00 

r.WT 

Gwl 

GWe 


GWS 

PGM 

NU 

0 






-0.3l0e 02 

-0. 6616-01 

-0.3D9E 

02 

-0.2686-01 

-0.3I0E 02 

4 

0 






-0.169E 01 

-0.3l0e-i)l 

-0.I67F 

01 

0.1066-01 

-0. I70F 01 

4 

2 






0.342C 01 

-0.182F-01 

0.341 F 

01 

0.228E-01 

0.341E 01 

4 

4 






0.330E 01 

-0. 112E-01 

0.329F 

01 

0.2126-01 

0.329F 01 

4 

6 






-0.201E 02 

-0.463E-01 

-0.200 E 

02 

-0.2526-02 

-0.2016 02 

5 

0 






-0.932E 00 

-0.1246-01 

-0.920F 

oa 

O.627E-03 

-0.926E 00 

5 

2 






0.149E 01 

-0.456E-02 

0.149E 

01 

0.139E-02 

0. 1506 01 

5 

4 






-O.lllE 02 

-0.2076-01 

-0.1116 

02 

-0.Z37E-03 

-O.IIIF 02 

6 

0 






-0.452E 00 

-0.6006-02 

-0.446E 

00 

0.4646-04 

-0.4486 00 

6 

2 






0.593E 00 

-0.1656-02 

0.595E 

00 

0.9756-04 

0.596E 00 

6 

4 






L 

LAT 

LAMBDA 

THFTA beta OENS 

TG 


HUGOS 

MLAN 


MCI 


4.61234 

21.94 0. 

5770 

i. 

.97 10.69 334 

0.1148 

22-02 

0.526F 

00 

0.284E 

00 



1 

2 

4 

5 

6 

7 

8 

g 


c THIS PROGftArt CGMPlilES TrtE VlHlSfLEtt HUDE WAVE AMPLtTljnE GAINS t IN CB.I, 
INTEGRATED OVER A R AV PATH v<H05F PARAMETERS AWE GIVEN 6Y ThF STAAPCRO 
VLF GROUP RAY TRACING PRUGRA*^. POSITIVE INEGATIVE) GAIN NEANS WAVE 
GROWTHinAMPING) . THE ENERGETIC ELECTRON DISTIBUTION IS EQUAL TO PE* 

I Al/A*Pn*H( P-PH)*EXP( Q^IP l/PK»NUO*SlNt AL PHA I «*2 ^3 T ) / P**NU ,WH.FP E 
alpha ANO P are respectively the electron PITCH-ANGLE ANO HOMEWTUH, 

H IS THF STEP EUNCTtONt Al AND Plt=lOOEVl ARE NORMALIZATION CONSTANTS, 
RT IS THE magnetic EltLO NORMALIZED TO THE EQUATORIAL VALJEt NU CAA 
take any value >4, Q AND NOO CAN TAKE ARPITRARV VALUES. 


INPUT; 


INITIAL ANO final GROUP TIMES 


NUtCiNUO tLINES 46 ANO 47) 

FREQUENCY OF THE WAVElKHZIr 
ILINFS 58 AND 59) 

L VALUE t LAT ITUOEI DEG ,) f NORMALIZED FREQUENCY, WAVE PROPAGATION 
ANGLEIDEG.), DEN5I TV ( CM-3) , GROUP TIMEISEC) I L INFS 69 AND 70). 


FREQUENCY OF THE WAVEIKHZ). INITIAL AND FINAL GROUP TIMES 
L VALUE, LATITUDEiDEG, » , NORMALIZED FREQUENCY, WAVE PROPAGATION 
ANGLE ( DEG. 1 . OE NS 1 TY I CM-3 ) , PLASMA FREQUENCY NORMALIZED TO 
ANGLEIQPG, I, DF NS 1 T Y< CM-3 I .GROUP TIMEISEC) UlNES 69 AND 73). 
GYROFRFOUENCY. PARALLFL COMPONENT OF REFRACTIVE INDEX, MINIMUM 
RFSONANCE energy FOP LANDAU AND FUNDAMENTAL CYCLOTRON INTER- 
ACTIONS IKEVI, CO, Dl. AND D2. ^ . 

GAINS integrated FROM INITIAL TIME, INCLUDING RESPECT [ VFL Y; 
HARMONICS M-0.1 ,2i LANDAU hARMQNI C ( M= 0 ) ; FUNDAMENTAL CYCLOTRON 
HARMON I CIM=U i 2ND CYCLOTRON H ARMON 1 C( M= 2 I ; PARALLEL PROPAGA- 
TION APPROXIMAT ICN, 

LIMITS FRFOUFNCY APPROXIMATION FOR REFRACTIVE INOEXi BUT 

INCLUDE IONS, RELATIVITY NOT INCLUDED. 

THE INTERVAL OF INTEGRATION IS OTATG=0.02 SEC 
Oft MORE ANO CAN BE CHANGED APPROPRIATELY IN LINES 75 AND 76. 
THF DISTRISUTION FUNCTION IS NORMALIZED TO DF=0.004, WHERE 
OE-RATIO HOT PLASMA DENSITY FOR PNERGVPlOOFV TO COLO PLASMA 
DENSITY) CHOSE APPROPRIATE VALUE OF DE IN LtNF 90. 

THIS PROGRAM IS APPROXIMATE AND ONLY VALID FOR NEAR 
PARALLEL PROP.\G AT ION, HARMONIC QF ORDER N CAN BE CONSIDERED 
AS APPROXIMATELY CORRECT WHEN DM<0,3 WHERE nM=D0.Dl.D2 !F 
M = 0 ,1 , 2 . 

CHOSE appropriate value DF PH IN LINE 137. 

fUNCTtCN FACTORIAL 
FUNCTTCN FFFINI 
IF IN,E0.0> GO TO 2 
FEF=l 

DO 1 1 = 1. N 

1 

GO TO 3 

2 FFF=I 

3 RETURN 
ENO 


10 

n 

12 

13 

14 

15 

16 
17 


FUNCTION ANG(NU,MO» 

ANC-l 

J0M=NQ/2*1 
on 11 JO^l.JOM 
11 ANO»IFl!nAn J0)-0,5 )*AN0 
ANC=ANQ*2* INU-3I/FFF1NQ/2 ) 
RETURN 
FNO 


10 function 3nqiun.no ) 

19 BNO-1 

20 J0M=N0^2FI 

21 OD 10 JO-l.JQM 

22 10 RNQ=BN0/(UN/2.-2. fFLOATI U0> 1 

23 BN0=FFF(NQ/2J *0NOy 2 

24 RETURN 

25 END 


TABLE A. 3. LISTING OF THE PROGRAM TO COMPUTE THE GAIN INTEGRATED ALONG 
A RAY PATH FOR NONSEPARABLE DISTRIBUTICNS AND SMALL WAVE 
NORMAL ANGLE. 
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26 FUNCTION TNOK UN^NO.PI 

27 BNOT - IUN-2J /UJN-2 > 

20 IF {NO.FO.OI GO TO & 

2q NCH ^ NO/2 

30 00 9 K=I,NOM 

31 11 BNOT - 0NOT ♦! -I FF M NO/2 ) 1/ (FFF « M l + FFF (N0/2-H H 

$ ♦ P* + (UN-2*2«H)/(UN-2<-2*MJ 
22 9 CCNTINUe 

33 8 CONTINUF 

3« RFTURN 

35 FND 


36 REAL KIKS 

37 REAL NC 

3 P RFal KIKLU 

39 RFAL K IK IKLANtKI KCYl 

AO RFAt 1A,K1G 

A1 INVL 

42 REAL LAT 

43 real K iKlfK IKC.MkCU 

44 R«=AL NU 

45 LOGICAL TRUC 

46 ftfAri5rl99l NU,NO»ON 

47 199 fnPMAT (3F10*6) 

48 kRITF(6,198J NU.NL'tUN 

49 198 FDPMATI * ^JU= • , F 10 . 3 r ' 0= * , F 1 0 .3 , ' NUO= ' , F 1(1 . 3 1 

50 P1PI=3,141592 

51 A»1,/<NU-3) 

52 6 = 0 

53 2000 N=N*^1 

54 ANFW ^ NO *FfF|fn / f F F F I 2*N + 1 ) « I NU-3 »N»OWU 

55 A = A ♦ AN»EW 

56 IF ( AA?( ANFUW AR9( A J •GT. .001) GO TO 2C00 

5T 2001 A r i./A 

C RFAP INPUT OATA FUCK COLO PLASMA RAY TRACING 

58 lODOl RFAO t 5 , 460 tFNO^t 9g j FKC » TGNEVi , TGL A$T 

59 460 format (F 10.4,F1D.6»F10.6 ) 

60 WBrTE(6*4550) FKC t TG NR N « TGL AST 

61 4550 FORMAT r 1FKC=‘ , F10.’4» * TGO^ • , F 1 0 »6 , • TGL A $ T- • , F 1 0 . 6 J 

C FKC IS FPF0U5NCY SIGNAL IN KH2 

62 TRUC = .e^ALSR. 

63 PG6SUP = 0. 

64 GWT$UM=0. 

65 GHL5UM=0. 

66 GWCSUK^C. 

67 GViLUSP-Or 

68 GHSSUM = n. 

69 GtiCUSM = 0, 

70 10000 MPADI5*490) INVL *LAT«L 1 tTKCTA, IOGNS^TG 

71 490 FCRMATIFIO ,5.F 10.2 *f i0.4,F10.Z, I 10 ,F10.6) 

72 TRUC ,= .NOT. TRUC 

73 IF ITRUC) TGOLO=rG 

14 I.FOr.TRUCI TuN£m=TG 

75 TTATG = AOS 1 TGNFW-TGniOI 

76 IF I OTATG.LT .0.02 I TA'.JC - .NOT. TrtUC 

77 IF I OT ATG.L T.O .02) 00 TP 10000 

70 OrCRAD = 3,141592/180 

75 TH'=TAP = ThFTA*DFGRAO 

9C BFTA - 1 A#8.9844SyftTIFL0ATnnENS) >/FKC 

81 CCl = Il-L A) #rAN( THFTAR ) 

82 CLAN = LA*TANITH=T AR) 

C lOFNS IS OFNSITY IN P A k T I C L t S /CU 0 1 C CfNTlHFTCRS 
C PFTA-PLASMA P'lEOUFNCY/GVftnFfiFagf nicy 
C LAtVAVE fRFOurNCY/GY^f^FFtECU''NCY 

C THFTA = ANGL*: PF PROPAGATIPN 4lTH RFSPFCT 1 I TUf MAGWETIC FirtO 
C END OF mn PLASMA PA^AMEICR CfCLARATION 
03 CnSTA = A3SICI.1S I THfcTAR n 

H4 SINTA - A!^S (SI M THETA'- n 

85 CFLTA =1-1 /(i, 3*LA) **=2 

66 C^^STAn = CCSTA-la*DELTA 

87 TALPHA-SIVTA/ ( Z’XC.OSfAl )) 

C THFTA HROF |S FXPRFSS'f^ nfpppps 
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an CALPHA-l/SOH T( l+TALPH4*+Z) 

B9 rind - F}i:t (CDSTA-LA*DeLTAI*LA) 

C R TNO-RFFRACTIVE INDEX 

90 Pl0=.O2 

91 PRtNf^=OeTA/$ORTt 1 l-LA)+L A) 

C PR|Nn=PARALLFL REFRACTIVE INDEX 
gz RPAft- RINH^COSTA 

93 rPFR=PIND*SINTA 

9A He: COSTA* (1/LA-LAJ 

95 Sr = BFta**Z / tLAAU~LA**ZM 

96 SP--(BFTA/LA)**Z 
C SP e: STIX R parameter 

97 PRlVtAT = II-LA)/< PPAR*LA1 

C PKIWAT IS the fund. NIIMENTUN OF PFSDNANCE 

c subfrocpam for parallel gain 

90 DF=.004 

99 PLAT LAr*DFGRAD 

100 3t T 5|N(RLAT» 

C lNVL=INVflP I ANT LATITODF 

f. DF-kATID HOT PLASMA DENSITY FOR eNERGV>lOOFV Tfl COLO PLASMA DFKSITY 
C iat=latitudf 

C BT=N0RMAL I MAGNETIC FIFL TO ECUATOPIAi VALUE 
I 01 BT=S0RTI l*^*X*A2) / C l-X«*2)**3 

IC2 PRP.A6SII.A-1 1/ (PRINC*LA) 

103 PPP = PIO/PRP 

ICA CnCFKC = -A*tPIPl/8l*0F 

105 PKIK = PPR**(NU-31 * ( P I P 1 /4 1 *DE *!1N0 1 NL,0) 

106 N-0 

1C7 1000 N-N+l 

ICP PRIKNl, = I 1/8T )**N«I I +N J • PP R *♦ I NU-3+N*ON 1 ♦ 

* I PlPt/4 I-sDE^lNO^+N/FFPINn » 6NC ( NU ♦N*ON . Z *N J 

109 PKIK = PKIK + PKIKNH 

no IF ( ABSI PRIXNWI/ AOS IPKTR J .GT..001I 

$ GO TO 1000 

111 K=1 

112 AKIK = I 1 /6T )**N*I-N/LA J * PPR** (NlJ-3*N*0N J ♦ 

t ( PIPl /A ) +nE*I N0*AN/FFF( M I ♦ BNQ ( NU»N*ON r Z*N 1 

113 lOCl N = N«-1 

114 AKIXNW . I l/RT)**N*(-N/LA J * PPP** I NU-3^-N*0N }♦ 

i I PIPl/A J*OF*(N0**N/FFF< N) J ♦ BNQ I NU+N*QN , Z*N 1 

11* AKIK = AKIK ♦ AKIKNW 

116 IF I ABSI AKIKNWJ/ ABSI AKIK) ,0T,,00lI 

i GO TO lOOl 

117 PKIK = A * (PKIK+AKIK) 

110 APVGC=SOPT(LA )* SORTH-LA) **3«2^6FTA 

119 AVGC-SORT a A* ( CDSTAO I*l4*ICnSTAn J**ZfS INTA**2) > / 

i ( 3FTA*CCSTA*2«LA* U-CFLTAJJ 
izo PWH«=-PKrK*APVGC*R INO 

C AVGC AND APCGC = MODULUS OF GROUP VELOCl TI CStGFNFPAL AND P ARAU FLCASES ) 
C NOPPALIZED TO C 
C WU'S APF IMIGMECA)/RF IOMEGA! 

C PK1K=R1K FOP OAPALtPL PPaPAGATION 

121 PCUEFG*- 1400CO*BT*LA*oR 1 NO /1NVL**3 

122 PG1000=PCnFFG*PK IK 

123 PPIMAT = I 1-LA)/IRPAP*LA » 

C Pftlv^AT IS THF FUND. MOMENTUM OF PESONANCF 
C PRl^eHOHFNTUM FPR PARALtFL PPCPAGATION 
12^ WO - 512. 

125 wCl = W0*1 PRIWAT J**2/Z 

126 PROWAT i L/RPAR 


/ 


133 


SEL 73-043 



127 

128 


12^ 

130 

131 

132 

133 
13^ 

135 

136 

137 
IBS 
13S 
lAO 
lAl 

1A3 

lAA 

IA5 

IA6 

1A7 

148 

149 
15C 

151 

1*2 

153 

154 


155 

156 

157 

158 
15S 
160 

161 

162 

163 

U4 

165 

166 


167 

168 

169 

170 


WLAN - WO+(PROWAT» ♦♦Z/? 

C PROWAT IS THF RESONANCE Hn»'F^TUM FOR LANDAU INTERACTION 
RWCIWL = WCl/HLAN 
C W0=P*0«'C«*2 

C WCl AND WLAN ARE RESPECTlVPLV ThE ENERGIES OF RFSONFANCF FCR 

C FUNRAMFNTAL ANO LANOAU RESONANCES UN KEVl 
DO ^ LA^SINTA/COSTA 

01 ^ I 1-LA I *SINTA/CUSTA 

02 = I 2-LA|*SINTA/CUSTA 
OO'ABS ID0> 

Dl^APSlOll 

D2=ABS ID2) 

WRITFI 6 r251 

25 FORMAT! *U L LAT LA ThFTA BETA Dtf 

»NS TC HUGOS WLAN WCl 00 01 02' > 

WRl TF 16,26 J |NVL,LAT,LA, THET A,BFT A , lOFNS ,TC, RPAR tWLAN, WCl ,DC, 01,02 

26 FORMAT {5Fl0.4«UO.Fl0.4,Fia.2«2ElK3f 3F5.2J 
CO = 2*APSISPi 

CIL !RPAR*P.PCR/C0J**2 
C2L - RPAR*RP6R*M#Sf:)7C0*«2 
C3l. ^ 1S0*H/C0>*'*2 

Cil = nBFTA/LAI*A2*u*CCSTA)frRPFft**2*LA/lJ/CD 
C21 = necTA/LA)**2A(l*CCSTAJ»RPFR**2ALA/Zl/C0 
PP0=1/RPAR 

PRI-^U-LAJ/ILA^RPARJ 

POP = PIO/PRO 

PIR = PlO/PRl 

PR2 = (2-L A>/UA«RPAR > 

P2R - P10/PR2 
C BEGIN LANDAU TERM 

SL = -NU*ICPSTA/LA-l >*IP0R»**INU-3J*4 
(C1L*8NC!NU+2,0J-C2L*6NCI NU,Z»*no+ PNO < NU-2 , 4 J ♦nO**2 »C 3L /4 J 
N=0 

300C N^N+l 

SLNW = -1 WBT) ♦♦N« INO**N/FFF(NM * ( NU+2»N+N «CN ) ♦ 1 PO R J # I NU-3+N«CN 1 
% «4« (COSTA/LA-l ) * I CU*BN0INU*'N*0N*2,2«N| - C2L*RN0<NU^N*CN , 

l2«N*2l*n0t C3CAFlMi5lNU»N4CN-2,4^2#Nl*D0«*2/4J 
SL - SL+SLNW 

IF I ABSI SLNW) 7A9S1 SL) ,CiT..OOn GO TO 3000 
C END lANOAU TERM 

SLU= -NU4K0ST A/LA-n*l PCU**lNU-3)44 ♦ 1 C IL»BN C 1 NU*2 , 0 H 

N-0 

2999 N-N41 

St UNW= -n /3T ) ♦*N*(NCl**N7f FFINI I ♦ I NU + 2*N+N4QNI ♦ ( P OR ) ♦ I N U-3 + N*CN t 
% *44(C0STA/LA-l ) * t CU«RNUNU + N*aN + 2,2»N) ) 

SLU - SlU+SLUNW 

IF 1A6SI SLUNVO/AHS ISL'J) .GT,,a01J GO TO 2999 
C BEGIN 1ST CYCLOTRON HARMONIC 

SC =r -pip**iisiu-ii « hcosta-lai/i 1 -la n ♦ cii**z 
% «(NU*0N0(NU,2) 50*01<'*2*NU*BNOINU-Z,4» I 

N=0 

3001 N=N^1 

SCNW = -PlR*«lNU+N*0N-3)*nC0STA-LA I/I 1-LAI )*C I l**2*l l/BT »**N 
( >♦ 

% 1 INU^Z+NtN^CNl^BNOl NU + N*ON,2*2*N)-4b04Dl**2>^lNU4^2«N^N*t}Nl 
i » FJNOt NU-2*N«'CN,4+2*N) I 
SC = SC^SCNW 

IF IAR51SCNW)/ABSI SCI .GT..OOD GO TO 3001 
N=l 

ASC = PlR**(NU+N*CN-3)* MCDSTA-LA) /I l-LAU*Cll*#2AU/flT)*«N 
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171 

172 


173 

174 


175 

176 
1 77 
178 


179 

180 
181 
162 


163 

184 


1 ES 
186 

187 

183 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 
2C5 
206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 


i * ( 1 1 * i NL)tN’*‘0N*2*N ) /LA-C. 5 0*D 1»'»2 

* *l2<'^i*6H0(NU+N*(5N-2« 2*M ♦^21/14 » ) 

3002 N-N+1 

A SCNk= PIS** tNU^N^QN- 3 I <» ( < COSTA-LA I / ( l*LA | ) *C1 1 1 / BT 1 ♦♦N 

% * (MO»*M/PFF(N n * \ NU*N+ONt 2«N ) /LA-C.5 0*01**? 

5 ♦(2*K*3M5<Ng + N*0^-2, Z+K+2)/LA) ) 

ASC = ASC+ASCNW 

IF 1 AftS(A9CNV<) /A 3S1 ASC > .Gr..OOU CH TO 3002 
C FNO 1ST CYCifTHON HA^IMQNIC 
C 0FGIN 7N0 CYCLDTRCN HARMflMrC 

SC2= -P2R**lNU-3> ♦ nCnSTA-LA)/(2-LAH * C2l**2 
i * NU*6!ViC(NU-2»4J*D2+*2/4 
N = D 

3D04 K'^^N+l 

SC2NW- -P2P**<NU+M*0N-31 *( (COST A-LM / I 2-L a j ) *C21**2*1 l/PT )**K‘ 
% ♦(NO*«N/FFF(N M « 

i 0 2*^2*(NU + 2*N*N*QN > fRNO (NU-2 + V*0N t4+2*N) 

SC2 = SC2 + SC2NW 

IF (A6S( SC2h»vn/A95i ISC2 ) ,GT,.00l J GO TO 3004 
N = 1 

ASCP= P29**(NU+iN*0N-3 (CnSTA-LA> / (2-LAJ » «C2 I **2< ( I / ftT 1 
$ *(NO»*?J/FFF(Nn* r02**Z/4)* 
i <,*K*riK5(N0»N*CN-2 »2«N + ?)/LA 
3005 N=N+1 

ASC2Nw=P 2B**(NU4-M*0N-3»*( (COSTA-LA J / (2-LAJ »*C2l**2*l 1/8TI**M 
$ *(N0**N,^FFF(NJ >♦ (D2042/41* 

% 4+N*BNCr^U+N*tfN-Z »2*N + 2 l/LA 
ASC2 = ASC2 * ASC2NK 

IF I AP.S(A5C2NWJ/ABS( ASC2 I ,GT, *0011 CD TD 3005 
C FNT ZfiD CYCLCTBON HARHOMlC 

CDFFKI = C0FFKC*CALPHA/CCSTA**2 

KUL » CnFFKl*SL 

KTKC - Ca£FKl*(SC+ASC» 

PlPHIN-2. 24 

PlPMI^) = PlP^n^7PlP 

IF (PI PMIM.GT.l. ) GO TD 293 

CN - BNOT ( NU.O.PIPHIN) 

M=0 

20C5 

CNNW = NO«6M*PlPW!^»•iM*CNJ*BNOT(NO♦M•ON^2♦MtPlP^<^M|/FFFCH» 

CN = CN+CNNM 

IF (ABSICNNWJ /ABSICM ,67,0.001) GO TO 2005 
P1PM|N=1, 

CD - BNOTINUrOfPlPMfN) 

M=0 

20C6 

CDNW - N0A*M«PIP«I N**(M«ONl48NOT(NO+M*ONi2*«fPlFMTN)/FFFl)«) 

CD = CD+CQNW 

IF (ABS (CONWl/ABS ICOI ,GT, 0*001) 60 TO 2006 
C06F = CN/CO 
KTKC = KIKC4C08F 
293 CONTTNUF 

KIK=K1KL+K IKC 

COEFG = -160000*LA#BT*RTNC/INVL+43 

GIOOO = COEFG » Kin 

KIKS = CDEFK[*(SC2+A$C2J 

61000L» COFFC * KI KL 

61C00C= C06FG * kiKC 

GIOOOS - CDEFG*KmS 

KIKLU - COEFKI*SLU 
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217 GLU = COFFG*KtKtU 

218 GCU - PGl 000*1 tCOSTA-LA»/H-LA) l*(CALPKA/COSTA*«2J*<ft[WO/PRlNDI 

* «Cll»*2*(PPft/PlP 

21^ KIKCU - GCU / COtFG 

220 APVGC=S0RT(LA)* SORT H-LA J+*3*2/ 8FTA 

221 PWIW=^-PKIK*APVGC*PRIM0 

222 Wl«T -KIK *AVGC*R1N0 

223 WEHt =- KIKL*AVGC*RIND 

224 XIWC = -KIKC«AVGC*RINn 

225 H[VILU = -KIKLU*AvGC»ft|ND 

226 WtViCU =» -K JKCU*AVGC*R1ND 

227 WIVIS » -K!KS*AVGC«R1N0 

C AVGC AND APCGC = MODULUS OF GROUP VELOCIT I FS (GENERAL AND P AR ALLELC ASES I 
C NORMALIZED TO C 
C WI«»S ARF INCCHFGAJ/REICMFGA J 
C PKIK-KIK FCR PARALLEL PROPAGATION 
22B PCOeFG=-lGOOOO*BT*LA*PRINO /INVL*»3 

229 PGIOOO-PCOEPG*PKIK 

23C COEFGk.- 06H6*FKC*6. 2831 64 

231 PGM = COEFGM+DTATG+PH IW 

232 CCC = COFFGk*DTATG 

233 GkT €CC*MIW1 

234 C'At ^ CCC + HIWL 

235 GWC = CCC*WlViC 

236 CWLU ’ CCC4MIVILU 

237 GWCU CCC+HIWCU 

23B GWS = CCC*N|WS 

239 PGHSUM - PGWSUM * PGW 

2A0 GVISSUM - GkSSUH ♦ CkS 

241 CHTSUM - GhISUM ♦ GWT 

242 GHLSUH = GWLSUM ♦ GWL 

243 GHCSUM - CUCStjM ♦ GhC 

244 GWLU5H = GkLUSM ♦ GWLU 

245 GkCUSw - CkCUSM ♦ GWCU 

246 PRINT 703 

247 703 FORMAT! » GHT GHL GWCl GWC2 PGWJ 

248 WRITE (6«430! GWTSUM rGWLSUM,GWCSUM«GWSSUM«PGHSUM 

249 43D F0BMAT(5Fll*3> 

250 IF ! TC.FO.TGLASTI GD TO lOOOl 

251 GO TO lOOOO 

252 999 STOP 

253 END 
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TABLE A. 4. AN EXAMPLE OF THE COMPUTER OUTPUT OF THE PROGRAM LISTED IN TABLE A. 3 


fKC 


- 0 . 


- 0 . 


- 0 . 


- 0 , 


- 0 . 


- 0 . 


-c. 


-D, 


-0 


-0. 


- 0 . 


= 80.OCOO TGO- 

0.CO0700 TGLAST- 0.524752 



L 

L6T 

LA 

THETA (3FTA 


PiPNS 

?.'3623 

47.9900 

0.3297 

-13.3700 1.0544 


611 

GkT 

f.Wl 

GW Cl 

GWC 2 PGW 



703C-C2 

-C.351F-C2 

-0. 352E- 

02 -0.373F-C4 -0.385F- 

0 2 


L 

IftT 

LA 

THFTA 86TA 


DFNS 

3*9t20 

43. 14CO 

0.5015 

-4.6200 1-5754 


562 


GWL 

GHCl 

GhC2 PGW 



<►860-01 

-0,52HF-D? 

-0,4330- 

Qt -Q.603F-Q4 -0.439F- 

■01 


L 

LAT 

LA 

THFTA *3FTA 


nFNS 

3.9150 

40.4900 

J.7404 

-6,6900 1.6979 


521 

GhT 

GWi 

GHCl 

GHC2 PGW 



2730 CO 

-C,1520-CL 

-C.257P 

Cu -0.8flOF-C4 -0.2586 

CO 


L 

i AT 

LA 

THETA 86TA 


DFNS 

3.U 1 

18.4ZC0 

0. 8364 

-10.9370 2.0098 


495 

GkT 

CUL 

GWCl 

GwrZ PGW 



lOF Cl 

-C.3B0F-01 

-0.134F 

01 -0.211F-C3 -0, 1296 

01 


L 

LftT 

LA 

theta RfTA 


OFNS 

3.6766 

36.6700 

0.f1832 

-13.6600 2.1865 


486 

Gk'T 

GkL 

GHCl 

GHC2 PGW 



5280 01 

-0.324E CO 

-0.496F 

01 -0.4CRF-C3 -0.433E 

01 


L 

1 AT 

LA 

thfta beta 


OEMS 

3.5797 

35.5500 

0.9CC0 

-15.0700 2. 2258 


465 

GhT 

gkp. 

GHCl 

GVIC2 PGW 



ICfeF 02 

-Q.633F 00 

-0.992F 

01 -0.5R5F-03 -0.7996 

01 


L 

LAT 

LA 

'H6TA bPTA 


OF NS 

3. 5549 

34. 15C0 

0.9116 

-lh.6070 2,2568 


436 

GWT 

CWL 

Gh!Cl 

G,<C2 PGW 



2<»4F CZ 

-C.142F 01 

-C. 230C 

02 -J.n92'^-C3 -0.l58f 

02 


L 

LAT 

LA 

THPTA 8FTA 


OENS 

3.37B4 

3 3. 29CiO 

0.9150 

-17.4600 2.2699 


486 

GViT 

GWL 

GHCl 

GWC 2 PGW 



• 357F C2 

-C.2170 Cl 

-0.336E 

02 -0.1146-02 -0.2216 

02 


t 

LAT 

LA 

tk^ta get a 


DENS 

3.2894 

32.2700 

0.9164 

-ia.42£)0 2.2604 


491 

GWT 

GhL 

GHCL 

CHC2 PGW 



,5190 02 

-0.352E 01 

-0.484E 

02 -0. 15 16-C2 -0.3136 

02 


L 

LAT 

LA 

THKTA RFTA 


DFNS 

3. 1840 

31.0200 

0.9154 

-19.5400 2.2871 


495 

GkT 

CKL 

GWCl 

GWC 2 PGW 



,7620 C2 

-C.614E 01 

-0, 7CIF 

02 -0.216F-02 -0.4446 

02 


L 

LAT 

LA 

thcta beta 


DFNS 

3.1245 

30,2900 

0.913U 

-20.1700 2.2901 


49 U 

GWT 

GWi 

GHCl 

GWC2 PGW 



.929F C2 

-C.016F 01 

-0.8476 

02 -0.263F-C2 -0.527F 

02 


L 

LAT 

LA 

thfta oeta 


UENS 

3.0579 

29.4500 

0.9112 

-2 0.9C00 2.2927 


502 

GWT 

GWL 

GWC 1. 

GWC2 PGW 




7G 

»<ucns 

WLAN 


HCL 

DO 

01 

02 

0.0215 

2.22 

0.51TE 

02 

0.Z14E 03 

0.08 

0*16 

0.40 

TG 

PUCCS 

WLAN 


WCI 

00 

01 

02 

0.C429 

3. 19 

0.251E 

02 

0.1306 02 

0.05 

0.03 

O.ll 

TG 

MUCCS 

WLAN 


WCI 

00 

Cl 

02 

0.0652 

4.35 

0.135E 

02 

0.166E 01 

0.09 

0.03 

0.15 

TG 

MUCCS 

WLAN 


WCI 

00 

ni 

02 

0.093G 

5.87 

0.743E 

01 

0.284E OO 

0.16 

0.03 

0.22 

TG 

MUCCS 

WLAN 


WC 1 

00 

01 

02 

n.l270 

7.57 

0.446E 

01 

0.7816-01 

0,21 

0.03 

0*27 

TG 

MUCCS 

WL AN 


MCI 

DO 

D1 

D2 

0.1529 

8.60 

C,33QE 

01 

0.408E-01 

0.24 

0.03 

0.30 

TC 

Mucrs 

WLAN 


WCI 

00 

Dl 

02 

0. 1907 

10.41 

0.2366 

01 

0.222E-01 

0.27 

0*03 

0.32 

1 G 

MUCOS 

WLAN 


HCl 

no 

Dl 

02 

0,2172 

11.39 

0.197E 

01 

0.17C6-01 

0.29 

Q.03 

0.34 

TG 

MUCCS 

WLAN 


WCI 

no 

01 

02 

0.2522 

12.45 

0.L65E 

OL 

0.1376-01 

0.31 

0*03 

0.36 

TG 

MCCOS 

WLAN 


WCI 

CO 

Dl 

02 

0.3004 

13.56 

0. 139E 

01 

0.119E-01 

0*32 

0.03 

0.36 

TG 

MUCPS 

WLAN 


WCl 

DO 

ni 

D2 

0.3315 

14.09 

0-129E 

01 

0. 1156-01 

0,34 

0.03 

0.40 

TG 

MUCfS 

WLAN 


MCI 

DO 

Dl 

02 

0.3704 

14.61 

0.L20E 

01 

O.l 146-01 

0,35 

0.03 

0.42 



APPENDIX B. POSSIBLE TRAPPING BY A NATURAL WHISTLER 


When a whistler is emitted by a lightning discharge of duration much 
shorter than the trapping time, the question arises whether the wave can 
change the trajectory of the particle enough to trap it. We make an 
analysis here on the assumption that the medium is one dimensional (axis 
Oz along . The interaction is only important around the equator, 

and the propagation is parallel (for nonparallel propagation the wave 
propagates across L shells and interacts during a shorter time with 
particles compelled to drift along magnetic field lines than for parallel 
propagation) . 

Cyclotron trapping ; 

In this case we assume parallel propagation. Let us suppose that 
at some point 

V|i = (B.l) 


where 


V 

R 


CO -CO 
c 


k 


(B-2) 


= parallel velocity of an electron • 

Along the lines of Dysthe [1971] and Sudan and Ott [1971], we will 
consider the particle as trapped (see Section 2B) if the condition 


d 

dt 


(v -V ) = 0 
II R 


(B-3) 


can be realized during a trapping time = 2:t (co^^kv^^) or corre- 

spondingly alor^ a length L = Vjj as the particle moves. 

We can write Eq. (B.3) as: 


Av. = Av 


R 


(B.4) 
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as the particle 


where Av,| and represent the variations of 


travels Az 


AV|| is the sum of two terms depending upon: 

1. the interaction of the wave on the particle (depends on the 

angle \|/ between and B^); 

2. the variation of the static magnetic field, i*e*i particle 


mirrorxng , 


( V v^ ^ dOD^ \ 

— sinilf - 7^ ^ j 


(B.5) 


dv is the sum of two terms depending upon: 

R 

1. the variation of electron gyrofrequency with distance (because 
of the medium inhomogeneity); 

2. the variation of wave frequency with distance (because of wave 
d is per s ion) 


J/3 

R I 2 m - 0 ) O) dCD / dz Ikv / dz R 

\ c pc' ^ g c / 


Av„ = 


^ V AZ (B.6) 


doi -1 5cjo 5a> 

^ “ V,| 3t ^ 


(B.7) 


(electron and wave are counters treaming ) . 
The whistler time delay is given by 


) r dz 

(Z.CD) = J — 


"O ® 


SO that : 


on g 


We see from Eqs . (B.5) and (B.6) that Eq. (B.4) is fulfilled if, 


l3 2 l\ 

'"'Ai i 'r (2 * -2^^ 53T ; d. *\v^ 


\ ‘'“o ri 
I Vo 
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Around the equator: 



(B.ll) 


and we assume u> ~ constant. 

P 

According to Eq. (B.ll), the variation of co is a second order 

c 

quantity in (z/r ) around the equator. The variation of cD is a first- 

ill 

order quantity but for o> -♦ 0; therefore, we can neglect the first term 
in the RHS of Eq* (B.IO) and express at the equator ^ 0 ) , 

We use the hyperbola approximation for the whistler dispersion [Bernard, 


1973 ] : 


E 


Om 

2 


03 

cE 

/o)(CD ^-co) 
cE 


D 


Om 


VpE 



(B.12) 


and for the D*E* model the values of the dimensionless parameters A and 
are [Bernard, 1973]: 

A ^ .25 , 

cIq - -7 . (B.13) 


Using Eq. (B.13) we get: 



03 ^(303-03 ^)-A03(0>fa3 ) 

cE cE cE 

3/2, T2 

U)/ (cO^g-03) 


(B.14) 


which can be inserted into Eq . (B.IO) (neglecting dtu^/dz) to yield: 


ktana CD , > 
cl — 



. 2. .1/2 1/2 

4C0 (CD^^-CD) CD^^ 


(B.15) 


Whistlers trigger emissions primarily at co ~ cd /2 and also at 

cE 

lower frequencies at the bottom of whistler traces. 
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As an example, we choose CC = 30^, L - 3, and 

X = 2jt/k - 2 km. Then Eq • (B-15) reads 

HI * 


O) 


O) 


cl -4 

Ss - 3.5X10 

cE 


(B.16) 


for 


--- 1.2X10 \b/m^, B, S 4xlo“^%b/m^ . 
OE '1 ' 


(B.17) 


The minimum field required to have trapping is excessively high 
compared to typical values (B^ 10 -10 Wb/m ). For a frequency 

where is comparable to, or higher than, its value at 

trapping would need the same order of magnitude for B^ • 

Low frequencies are not more favorable either because here 

CD 

k — 7 = Jl^ decreases with frequency and v_ increases. This makes it 
qJW ^ R 

c 

more difficult for inequality (B-10) to be fulfilled. 

The RHS of Eq. (B.IO) may vanish at a point off the equator. In 
this case approximations (B.ll) to (B.12) are no longer valid, and a 
more refined analysis should be made. 

Landau trapping : 

The Gendrin mode is the most favorable for trapping because the 
group velocity and resonant velocity are equal and ray direction and 
particle drift directions are the same. As we are concerned here with 
ducted whistlers, we will suppose cos6 ^ 1. 

Computations similar to those above yield: 


CD 

'^R " k 


(B.18) 


Av = 


eE 

1 1^ 

M V.. 


sin(kz-CDt) 


V dCD 

1 _1 C 

v,^ 2CD dz 
If c 


Az , 


(B.19) 
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(B.20) 



Vi 

1 

1 

dCD 

p 

dCD 

\ — 

/J_ 

1 ' 

j dCD 


[(" 

CO -CO 

c 

CD 

p 

dCD J 
c ' 

1 dz 

1 kv 
\ g 

’ ^ ; 

) dz 


A condition similar to Eq. (B.IO) is found for the minimum 
field amplitude required for trapping. 


eE^k 

"m 




/ n dCD \ 

dCD 

/ \ 

/, 2 CO <0 p \ 

y ^ 2co ” CO dco / 

\ c p c / 

c 

1 K T \ 1 

dz 

< 1 
OQ 1 

1 

y 

< 


(B.21) 


This time, v ^ v and the wave dispersion can be neglected, and 
R g 

the space inhomogeneity is more effective in preventing trapping. 
Neglecting the wave dispersion, Eq • (B.21) can be fulfilled up to a 


certain distance z from the equator given by: 

-1 


— < (O) V ) ^ 

^ CE 
E 


2 oj 
tan a — 
2o) 


dCD 

0^ p 

CD dCD 

P C 


2neE^ 

9\ MR^ 
w E 


(B.22) 


Let us take the same example as previously = 2km, cd = 
a = 45*^, and L = 3 . 


E_ = E ^ V B, sin0 (B.23) 

1 z pi 


Eq. (B.22) reads 


z 10 

— < 1.6X10 B sin9 
— 1 
E 


(B.24) 


There will be trapping if this value is bigger than a trapping length. 
For a ducted whistler, 9 < 30° and trapping would require a minimum 
field of 

~ 4xl0"^\b/m^ (B.25) 
This field is still high compared to typical values. 
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APPENDIX C. HOT PLASMA CORRECTION TO THE REAL PART 


OF THE REFRACTIVE INDEX 


The refractive index is given by the expression: 

2 oo r 


= 1 




2 2 
k V 

1 


2(a>kv) 


(kV|. “0>f(X ) 


2 ky, -cjyhCD 

i 




For small v ,, we treat 


kv. 


CD -CD 
C 


as a first order quantity and we arrive 


at : 


2 2 2 
_ 03^ CD a) ^ ^ CD 

= 1 + -T— ^ ^ ^ , k < V > 2 


2 

— k^ < !i > 

2 2 ^2 ^ ’ 


(C.2) 


where : 


< V 


■ > =f 

>1 Jq J- 


We rewrite Eq . (C.2) as: 


N 


N (1 + e„) , 
c n 


(C.3) 


where : 


^ ^ 


O-A)^ ' ' 2c/ 


and where N is the cold plasma refractive index expression, 
c 

2 1 2 

For isotropic distributions, < '^ii > - o ^ ^ reduces 

•I 2 J, H 


to 


V 2 
II ^ 


-. 0 A ^ 

^ 3 ^ ^ * 

“ 2(1-A)^ 

This result was already derived by Guthart [1964] where he chose for 
a Maxwellian distribution for which < v ^ 2 KT/m. 

Equation (C.4) is also given by the full adiabatic macroscopic 
theory [Quemada, 1968]. 


(C.4) 
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We choose now Eq, (3.6) for 

- ^3q^yq -q/2 \ 

®H ■ (v-5)(q+3) 


Then 


1 - 


(1-A)(q+2)1 


“J ^EH 


m 


(C.5) 


As an example, we chose \j - 6 and q = 6 (the distribution is thought 

to be quite anisotropic towards low energies) and evaluate roughly the 

corrective term e at the equator where the contribution to whistler 
HE 

dispersion is most important: 

With W„ = 10 eV, at L ~ 4, g ^ 10, ~ 10“^ (for v = 6) 

H EH 

and 

®HE « 1) ~ - 4X10"‘^, 

%E ~ " 4X10"^. (C.7) 

Such a difference in the dispersion would be unnoticeable . 

Unfortunately, a similar computation for 0^0 becomes fairly 

complex as a cubic equation has to be solved. This was done by Wang 

[1970] for the scalar pressure theory [Denisse and Delacroix, 1963]. He 

found corrections became important close to the resonance cone, as 

expected. A solution for Q 0 will contain terms such as Eq. (C.6) but 

where 1-A is replaced by cos0-A snd e increases very rapidly for 

H 

9 =* Wang's [1970] computations could be improved using the full 

adiabatic theory [Quemada, 1968]. 
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